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The Specific Heat of Liquid Helium tT 


By R. W. Hit and O. V. Lounasmaa uh 
Clarendon Laboratory, Oxford 


[Received November 27, 1956] 


ABSTRACT 


New measurements of the specific heat of liquid helium between 1-80 
and 5-05°K are presented, and some of their consequences for the helium 
vapour pressure scale of temperature discussed. 


§ 1. IyrRopvuctIon 


TE work to be reported here forms part of an investigation of the entropy 
diagram of helium in the range 2 to 20°K and 0 to 100 atmospheres pressure, 
It is presented in advance of the main work because of its bearing on the 
helium vapour pressure scale of temperature, a subject of considerable 
current interest. In determining the diagram, entropy differences only are 
measured, and before absolute values can be given it is necessary to know 
the value at any one point. For our purposes, the most convenient 
reference point would be the liquid in equilibrium with the vapour at a 
chosen temperature, but the entropy is not known with sufficient accuracy. 
It was for this reason that new measurements of the specific heat were 
undertaken. 

Below the A-point, the most recent data are those of Kramers et al. (1952) 
between 0-20 and 1-85°xk, and those of Hercus and Wilks (1954) from 1-10 
to 2:05°K. In the immediate vicinity of the \-point the best data are 
those of Keesom and Keesom (1935), while at higher temperatures there 
exist no more recent measurements than those of Keesom and Clusius 
(1932). Two difficulties arise: first, while the values of Kramers et al. 
and Keesom and Keesom fit smoothly together, those of Hercus and 
Wilks conflict with both; second, the results of Keesom and Clusius 
show a considerable scatter, and are said to be uncertain to 3%. 

A new determination of the specific heat above the A-point was evi- 
dently desirable, and could conveniently be undertaken in our apparatus. 
Measurements below the A-point were not foreseen when the apparatus 
was designed, and presented some difficulty. However, the nature of the 
disagreement between the results of Kramers ef a. and those of Hercus and 
Wilks strongly suggests the presence of an undetected systematic error in 
one of them. <A few points in a limited temperature range should then 
suffice to determine where the error lies. This was, in fact, as much as 
could be done without redesigning the apparatus. 

— + Communicated by Dr. K. Mendelssohn, F.R.S. 
t On leave from Turku University, Finland. 
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§ 2, EXPERIMENTAL 


The colorimeter was made of copper and was filled as far as possible 
with fine copper wire, leaving a volume of 21-3 cm* for the helium. A 
filling tube of german silver, 0-3 mm internal diameter, was used ; it ran 
directly from the calorimeter to a valve at room temperature, thereby 
keeping the dead volume of the system very small. The calorimeter 
was suspended inside a vacuum jacket immersed in liquid helium at its 
normal boiling point. Exchange gas was used to cool the calorimeter 
to liquid hydrogen temperatures and was then removed, and further 
cooling obtained by condensing liquid helium into a separate compart- 
ment of the calorimeter itself. By pumping this helium away, tempera- 
tures lower than that of the bath could be achieved down to a limit 
imposed by the size of the vessel. With large helium samples this limit 
was a little above the A-point, in accordance with the use for which the 
system was designed. When it was later decided to carry out measure- 
ments of the specific heat below the A-point, it was found that 1-8°K 
could be reached with a reasonable size of sample, and this was considered 
to be adequate for our purpose. 

The particular virtue of this method of cooling is that exchange gas is 
used only at comparatively high temperatures where it is easily removed. 
Good thermal insulation is therefore obtained, and the temperature drifts 
encountered were only a few millidegrees per minute at worst. The small 
correction for these drifts is easily made as thermal equilibrium is estab- 
lished in about 10 seconds. As the critical point is approached the 
equilibrium time increases, but does not exceed one minute. 

The thermometer used was a carbon composition radio resistor (47 ohm, 
4 watt, L.A.B.). For use below 4:2°K it was calibrated during every 
experiment against the helium vapour pressure on the scale designated 
by 755, (Clement et al. 1955). No change of calibration could be detected 
during an experiment, and changes between experiments were small. 
This constancy was particularly valuable at high temperatures, where 
the thermometer was calibrated by using the calorimeter as the bulb of 
a constant volume gas thermometer. Calibrations of this type could not 
be carried out during a specific heat experiment, and adequate correction 
for changes is possible only if they are sufficiently small. 


§ 3. RESULTS 


The saturation specific heat C, cannot be measured with a calorimeter 
completely full of liquid, and an important correction must be applied for 
the heat supplied for evaporating liquid and heating the vapour. This 
correction can be obtained by elementary thermodynamics, yet the 
derivation does not appear to have been given in the literature of liquid 
helium. 

Let S, be the entropy of n moles of liquid and vapour in equilibrium, 
then S,=n,S,+n,S, where n, and n, are the number of moles of liquid 
and vapour respectively, and S, and 8, are the molar entropies. Also 
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let V, and V,, be the molar volumes and V the actual calorimeter volume. 
Then S,=n,S 1+(n—n,)(S,+L/T) where L is the latent heat at tempera- 
ture 7’. Using the relations (n —n,) V +n V {= V and L/T=(V,— V)dp/dT 
we have S,=nS,+(V—nV,)dp/dT, so that 

C,=Td8 aT =[C,—Td/aT {(V —nV,)dp[aT }\/n. 
In this expression the total specific heat C', must, of course, be corrected 
for the heat capacity of the calorimeter. 

Below the A-point the vapour correction is generally small, but it 
becomes important at higher temperatures where it can be positive or 
negative, depending on the value of n in relation to V. Some examples 
of this are given in table 1 ; here, as elsewhere, gram units are preferred. 


Table 1 


A=Td/dT {(V —nV,)dp/dT} (A and C, in joules deg ; 
C,, in joules deg! gm-*) 


to molar ones so as to conform to the general usage. The experimental 
results are shown in fig. 1 ; no distinction is made between points obtained 
with different values of x, but the scatter does not exceed 19%. The values 
of Keesom and Clusius, recalculated to the 7';;, temperature scale, are 
included for comparison, and may be seen to scatter fairly evenly about 
the present ones. It is evident, however, that they are not accurate to 
3% near the normal boiling point. Values obtained below the A-point are 
also shown in fig. 2, where the larger scale affords easy comparison with 
those of Kramers ef al. and Keesom and Keesom, represented by the 
solid line, and those of Hercus and Wilks, represented by the dotted line. 
All results are in terms of 7';5,. It is obvious that our few points confirm: 
the Leiden values ; therefore, although it is still not clear how the cae 
crepancy between the two sets of results has arisen, we have taken the 
i to be the correct ones. 

ies aan of the results above the A-point is limited ipearaien! 4 
the accuracy with which the sensitivity of the eo daca co e 
determined ; this is about 0-5°%, and the specific heat va. oe an a 
sidered to be accurate to 1% or better. Below the es ae oe 
insulation is less satisfactory, due to the helium film in x e ig oe 
This is probably the reason for the greater scatter in the experimen 


results. 
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Fig. 1 
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In calculating the entropy values, the value of 0-535 joules deg-! gm—1 at 
1-80°K (Kramers e¢ al.) has been taken as a starting point. Thereafter, 
the present results were used exclusively. Entropy and smoothed 
specific heat values are given in table 2. 


Table 2. The Specific Heat and Entropy of Liquid Helium under its 
Saturation Vapour Pressure 
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The unit for C, and S is joules deg-! gm-! 


§ 4. Discussion 


At the present time, physicists are faced with two helium vapour 
pressure temperature scales, one of which, proposed by van Dijk and 
Durieux (1955), is based on the thermodynamic equation relating the 
vapour pressure to the thermal properties of the liquid and Meigie 
phases. This scale has recently been criticised by Keller (1956), w ~ 
article contains a full account of the methods by which the two scales 
have been formulated. The basis of Keller’s criticism is that the bh 
Dijk scale is uncertain within rather wide limits because of hile ies 
in the relevant thermal properties. In particular, the pono ¥ Mas. 
the Kramers and the Hercus and Wilks specific heats leads ie core 
about the value of Lo, the latent heat at absolute Zero. 4 he ‘ 
possible choice of values for the liquid specific heat above the mete 5 
within the range of the Keesom and Clusius data, leads to uncertain 
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which are much greater than the supposed difference between the scales. 
Keller therefore concludes that the apparent close agreement of the two. 
scales is to some extent fortuitous. 

The present data should help to reduce these uncertainties in two ways. 
Firstly, the confirmation of the specific heats of Kramers et al. would 
appear to justify van Dijk’s neglect of the Hercus and Wilks data when 
estimating L,. Secondly, the increased accuracy in the specific heat 
above the \-point materially reduces the uncertainty in the 3—4°K range. 
The full consequences of our data for this part of the van Dijk scale will 
not be evident until a full recalculation can be carried out, but, since the 
uncertainty here is almost directly proportional to the uncertainty in 
the liquid entropy (Keller), it seems that at least a threefold increase in 
precision should be obtained. 
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“Li Emission from Heavy Element Disintegrations Produced 
by 5-7 BeV Protons + 


By 8. J. Goupsack, W. O. Lock and B. A. Munir 
Department of Physics, Univetsity of Birmingham 


[Received August 24, 1956] 


SUMMARY 


A stack of Ilford G5 emulsions exposed to the 5-7 Bev proton beam of 
the Berkeley Bevatron has been searched for disintegrations of silver and 
bromine nuclei with ejected *Li fragments, and 91 examples found. The’ 
experimental results can be explained qualitatively by assuming that 
most of the fragments from stars with 12 or more black prongs are ejected’ 
as part of the evaporation process but a few of the fragments from the 
smaller stars (less than 12 black prongs) are ejected during the initial 
stages of the collision. 


§ 1. IyTRopDUCTION 


A PREVIOUS paper from this laboratory (Munir 1956) has given data on 
the ejection of *Li fragments from nuclear disintegrations produced in G5 
emulsions by 950 Mev protons. Through the kindness of Professor E. J. 
Lofgren at Berkeley, California, a stack of G5 emulsions was exposed for 
us to the 5-7 Bev proton beam of the Bevatron. These emulsions were 
searched for stars with ejected *Li fragments as in the work at 950 Mev. 
The high background in the emulsions, and the difficulty of establishing 
a reliable criterion to distinguish between light and heavy element dis- 
integrations, led us to confine our attention to definite heavy element 
disintegrations. 


§ 2. EXPERIMENTAL PROCEDURE AND RESULTS 


Ilford G5 emulsions, 400 » and 600 » thick were used. The primary 
proton beam was in the plane of the emulsion and parallel to one edge of 
the plates. The emulsions were searched under a x 10 objective and 
15 eyepieces. Each star with more than 5 black prongs and with a 
minimum track in the beam direction was examined under a x 45 oil 
objective. All black prongs were followed either to their ends in the 
emulsion or to the point at which they left the plates. The precautions 
were in general the same as those observed previously (Munir 1956). All 
the ‘ hammer tracks ’ found were assumed to be due to *Li nuclei stopping 


in the emulsion. 
0S ee ee 


+ Communicated by P. B. Moon, F.R.S. _ 
+ Now at-the American University of Beirut, Lebanon. 
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In all, 7700 stars with 5 or more black prongs were examined and 91 
examples of *Li emission found. Of the parent stars, 90 had 6 or more 
black prongs and must therefore correspond to heavy element disinte- 
grations ; the remaining parent star had 5 black prongs but a total of 14 
prongs and must represent a heavy element disintegration. A few of the 
7700 stars examined had only a small number of associated grey prongs 
in addition to 5 or 6 black prongs and may, of course, correspond to light 
element disintegrations. Thus the figures for the frequency of *Li emission 
are probably lower limits. 

It should also be pointed out that in this paper we are dealing with a 
slightly biased selection of heavy element disintegrations, for our method 
of searching the emulsions would miss a heavy element disintegration 
consisting of a very small number of prongs. Such interactions could 
only be found by seaching along the tracks of the primary protons. On 
the other hand, at a primary energy of 5-7 Bev the number of small heavy 
element stars will clearly be rather low. Thus the bias is probably not 
serious. 


Fig. 1 
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Energy distribution of the ®Li fragments. 


78 


The energy distribution of SLi fragments is shown in fig. 1. It has a 
cut-off at 6 mev. This value may be taken to be the height of the potential 
barrier facing a *Li nucleus in a silver or bromine nucleus when excited by 
the impact of a 5-7 Bev proton. The most energetic fragment recorded 
had an energy of 78 Mev. ‘The distribution has been corrected in the usual 
way for loss of fragments leaving the emulsion before coming to rest. 

The angular distribution of the fragments is plotted in intervals of equal 
solid angle in fig. 2. There are 55 fragments (observed number) in the 
forward direction (angle to primary direction less than 90°) and 36 in the 
backward direction. 
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The disintegrations with ejected ®Li fragments have been divided into 
groups with different numbers of prongs and the frequency of emission of 
*Li fragments calculated for each group. These frequencies together 
with other relevant data are given in the table ; they have been corrected 
for loss due to fragments leaving the emulsion before coming to rest. The 
frequencies are plotted as a function of the mean excitation energy in fig. 3; 
the full line is that calculated by LeCouteur (private communication) on 
the basis of evaporation theory. 


Fig. 2 
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§ 3. DISCUSSION 


Figure 1 shows that most of the fragments have energies between 6 
and (say) 40 mev but that there is a pronounced high energy tail. Suppose 
now that we separate the data into two groups depending on the number 
of black prongs in the parent star. Figures 4(a) and 4(6) are plots of the 
energy distributions, of the fragments from stars with less than 12 black 
prongs (‘small ’ stars) and with 12 or more black prongs (‘ large ° stars) 
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Frequency of ejection of the *Li fragments plotted as a function of the 


excitation energy of the parent stars; the solid line is that calculated 
on the basis of evaporation theory (LeCouteur, private communication). 


respectively. The mean excitation energy for each of these groups is 
300. mev and 550 mev respectively ; there are 52 observed events in the 
first group and 39 in the second. 
If the *Li fragments are ejected during the ‘ evaporation process’ then 
their energy spectrum N(#)dH should be of the well-known form 
N(£)dE= eM ly) 


exp 


Le 4b 


dH 
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where 7’ is the nuclear temperature and V the potential barrier. In the 
present instance we may take V—6 mev (see § 2) then 7’'=4-3 mev for 
U=300 mev and 7-8 mev for U=550 mev (see Lock et al. 1955). The energy 
spectra for these two excitation energies have been calculated using the 
above formula and plotted as full lines in figs. 4(a) and 4(b), normalized 
to the number of *Li fragments in each plot. The curve in fig. 4(b) 


Fig. 4 
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(a) Angular distribution of the *Li fragments ejected from stars with less than 
12 black prongs ; the shaded area corresponds to those fragments with 
more than 36 Mev kinetic energy. (6) Angular distribution of the SLi 
fragments ejected from stars with 12 or more black prongs, 
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the fragments from the large stars are ejected as part of the evaporation 
process ; this is, in fact, borne out by their angular distribution, shown in 
fig. 5(b), which is essentially isotropic, as it should be for an evaporation 
process. On the other hand, a substantial fraction of the fragments from 
the small stars have energies in excess of those predicted by the evapo- 
ration theory. Further, their angular distribution (fig. 5(a)) is markedly 
anisotropic. This strongly suggests that fragments are also being ejected 
during the initial nucleon—meson cascade in the struck nucleus ; such 
fragments would be ejected in the forward direction and would have 
relatively high energies. To be more specific we suppose that the 12 
fragments (observed number) in fig. 4(a) with energy more than 36 Mev 
arise from the cascade process. These 12 events are shown shaded in 
fig. 5(a); their forward to backward ratio is 9:3. The forward to back- 
ward ratio of the remaining 40 events is 25 : 15=1-7-0-6 compared with 
34 : 18=1-9-L0-6 for all 52 events. The value of 1-7--0-6 is not in dis- 
agreement with a value of unity for this ratio; we have made no correction 
for centre of mass motion which may be appreciable in some cases. 


Frequency of Ejection of Li Fragments 
(a) At 5-7 Bev 

Figure 2 shows that the observed frequencies of emission at 5-7 Bev 
(points c and d) of the *Li fragments are higher by a factor of about two 
than those predicted by LeCouteur on the evaporation theory. This 
discrepancy remains if one considers only the fragments emitted from the 
‘large’ stars which, as shown above, are otherwise consistent with the 
evaporation hypothesis. 

On the other hand, we may compare the relative frequencies of emission 
as predicted by LeCouteur for 300 and 550 Mev excitation energy (‘small ’ 
and ‘large’ stars). If we suppose that all the 42-7 fragments from the 
1950 ‘ large ’ stars of 550 Mev excitation energy are evaporation particles, 
then from the 5750 ‘small’ stars of 300 Mev mean excitation energy one 
should expect 49-4 ‘evaporated ’ *Li fragments. The corrected number 
of observed fragments with energy less-than 36 mev is 42-1. Thus the 
remaining 17-2 events of higher energy probably arise from some other 
cause such as the ‘ cascade ’ hypothesis which we have discussed above. 

We may perhaps regard the ‘ small ’ stars as those in which the primary 
proton makes an edge collision with the struck nucleus of silver or bromine. 
Thus the resulting disintegration is rather small and any aggregate of 
nucleons, such as a ®Li fragment, which by some means may acquire some 
tens of ev, has a good chance of escaping from the nucleus without making 
a collision. On the other hand, the ‘large’ stars may represent those 
cases in which the primary proton makes a head-on collision with a nucleus. 
Any nuclear fragment accelerated in the initial nucleon—meson cascade 
now has a much smaller chance of escaping intact from the disintegrating 
nucleus. We may expect “Li fragments to arise largely. or perhaps 
entirely, from the ‘evaporation’ stage of the disintegration. The 
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discrepancy between the observed and predicted frequencies of emission of 
the fragments may be overcome in a further refinement of the theory or 
it may arise through an underestimate of the excitation energy involved. 


(b) Comparison with Frequencies of Ejection at Other Knergies 


Point a on fig. 3 is the frequency of ejection of ®Li fragments at 340 Mev. 
This was calculated from the results of Wright (1950) for the cross section 
of ®Li ejection from silver and bromine nuclei combined with the inelastic 
scattering cross sections for these nuclei at 305 mev determined by Millburn 
et al. (1954). Point b corresponds to the results of Munir (1956) at 950 Mev. 

Both points a and b are well below the theoretical curve given by 
LeCouteur, which is in contrast to the situation at 5-7 Bev which we have 
discussed in § 2. The origin of the discrepancy between the experimental 
data and the predictions of LeCouteur, based on the evaporation theory, 
is not clear. 


§ 4. CONCLUSIONS 

SLi fragments are ejected from the nuclear disintegrations of silver and 
bromine produced by 5-7 Bev protons with a much greater frequency than 
that expected on the basis of the usual evaporation theory. The results 
are consistent (expect for the absolute frequency of emission) with the 
assumption that the fragments from the stars with 12 or more black 
prongs arise almost entirely from the evaporation process, while those 
from the stars with less than 12 black prongs are partly ‘ evaporation 
particles ’ and partly ‘ cascade’ particles. The exact mechanism which 
gives rise to the high energy (36-70 Mev) forward moving fragments 
remains to be elucidated. 
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ABSTRACT 


An experiment has been carried out at Saskatoon, Canada, at a nominal 
geomagnetic latitude of 60-5°N, which shows the form of the low energy 
end of the cosmic ray alpha particle spectrum independent of the earth’s 
magnetic field. The differential spectrum shows a maximum at about 
300 Mev per nucleon. The number of particles on the low energy side of 
the maximum is too great to be accounted for by secondary production, 
and it is believed that this spectrum represents the differential spectrum 
of primary alpha particles outside the earth’s magnetic field. <A re- 
analysis of previous experiments at a lower geomagnetic latitude, 55°N, 
in the American continent, gives a spectrum very similar to the one 
obtained in Saskatoon. As a consequence of this it appears that the cut- 
off energies experienced by cosmic ray particles over the American 
continent are probably lower than those predicted by geomagnetic 
theory. 


PART I 
§ 1. INTRODUCTION 


Because of the uncertainty introduced in attempts to infer the primary 
spectrum of cosmic rays from latitude measurements, it was felt that a 
direct experiment to determine the energy spectrum of the various 
components should be undertaken. Consequently, a series of balloon 
flights was carried out at Saskatoon, Canada, in the summer of 1954. 
This period was a time of minimum solar activity. The flight from 
which all of the information has been obtained so far was flown on the 
18th June 1954, during the daylight hours. The experiment consisted of 
exposing a stack of 83, 3}in. by 93 in. 500-micron thick emulsions to 
cosmic rays at high geomagnetic latitude. The latitude of Saskatoon, 
Canada, was chosen in the following way. On the basis of then-current 
balloon techniques, and with the payload determined for this flight, the 
minimum pressure which was felt to be practical was 10 g/cm’. This 
pressure determines a range cut-off imposed by the air above for all 
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components of the primaries, and if one is to study energy spectra in the 
absence of geomagnetic effects, the air cut-off for the lowest Z component 
to be studied must be equal to or greater than the energy cut-off imposed 
by the magnetic field at the latitude at which the experiment is to be 
carried out. The residual pressure which could be reached in the balloon 
experiment thereby specified an equivalent geomagnetic latitude for 
alpha-particles of 60:-5°N. In view of the results which follow, it appears 
that Saskatoon may be at an even higher latitude than 60-5°N. 

The detecting emulsion block was mounted in a frame which could be 
rotated by a small motor activated by a pressure device. The emulsion 
stack was enclosed in a pressurized spherical aluminium gondola. The 
instrumentation for the flight included an aneroid device which activated 
the motor which turned the photographic emulsion stack through 180° 
about a horizontal axis after the balloon had reached ceilmg. The same 
signal which turned over the emulsion stack also opened the balloon 
appendix in order to insure as high an altitude in the flight as possible. 
In addition to this turnover switch, there was a ballasting servo system 
which ‘followed up’ as the balloon rose and set itself to ballast if the 
balloon descended as much as one millibar below its maximum ceiling 
altitude. There was also a ‘ down-camera ’ for photographing the ground 
and determining the trajectory, as well as a telemetering system utilizing 
a temperature controlled bellows which transmitted the pressure measured 
during flight. In addition, a camera photographed a Wallace and 
Tierman pressure gauge to give an independent record of the altitude of 
the flight. 

All of the methods of pressure determination in the flight agreed to 
within 0-5 of a millibar, and the mean pressure at which the plates were 
exposed in their upright position at altitude was 10-8-+-0-6 g/cm?. 
Because the plates were in one position during ascent and were turned 
over at ceiling, one can separate those particles for which the direction of 
motion is known into a group that came in at ceiling altitude, and a group 
which came in during ascent. 

Figure 1 shows the time—altitude curve for this flight. It can be seen 
from this curve that the period during which the plates were in their 
upright position represents a very constant level region of the flight. 

Because of the importance of the determination of the charge of particles 
passing through the stack, a photometric device was developed to measure 
the ionization over a wide range of energy loss. The profile of the track 
is displayed on an oscilloscope and measurements of the track can be made 
with various slit widths. This measurement gives the equivalent 
obscuration of the track as compared with that which would have been 
produced by a wire of the width indicated in microns. Figure 2 shows 
the way in which this densitometer measurement varies with dE/dx, and 
for convenience the energy loss for charged particles at different speeds 
is also indicated in this graph. This figure shows the general characteristics 
of the G—-5 emulsion in that the obscuration between about 20 and 100 
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times minimum ionization varies very little. It is desirable, then, to 
make ionization measurement. with the densitometer in regions of the 
track where this plateau is not important. This condition is satisfied 
for the Saskatoon alpha-particles. 


§ 2. ALPHA-PARTICLES AT SASKATOON 

The plates were scanned under x 360 for particles which satisfied the 

following criteria : 
1. Ionization greater than 3-5 minimum. 
2. Projected length in each emulsion greater than 8 mm. 
3. Zenith angle less than 45°. 

Two regions of the stack were scanned. One entire vertical side was 
scanned for all particles which satisfied the above criteria and which 
entered the stack from the outside. The other region scanned was a 
line parallel to and 1 cm below the top surface of the stack. All particles 


crossing this line which satisfied the above criteria and which entered 
from the top of the stack were recorded. The total scan yielded 190 


Fig. 3 
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alpha-particles, plus four fast heavier nuclei, and approximately 1500 
contamination particles, mostly slow protons. 

Figure 3 shows a densitometer calibration appropriate for the alpha- 
particle scan. Measurements were made under a linear magnification of 
100 times. All the tracks found were scattered by the Fowler coordinate 
method with 500 micron cells, and 20 of the slow alphas were followed 
until they ended. 

Figure 4 shows the resolution between protons and alpha-particles 
obtained by plotting the densitometer measurements versus the scattering 
parameter obtained without noise elimination (‘raw ’ scattering). The 
lines are the theoretical lines based on a scattering constant of 26 and the 
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known densitometer calibration. There are approximately 15 second 
differences on each track and an equal number of densitometer readings 
were taken throughout the emulsion, from the air surface to the glass 
surface, to average out the depth variation of this measurement. 

To investigate the errors in the ionization waenoasnet Aree 
ments were made on the track of a fast particle which passed thr ous a 
emulsions. Using the dependence of width upon ionization the ae ec 
standard deviation of the energy determinations of alpha-particles from 


20% 1 interval ~ 
the ionization measurements was found to be 20'% in the energy 


162 P. H. Fowler et al. on the Low Energy End of the 


100 to 400 Mev per nucleon, and of the order of 30°%, between 400 and 
500 Mev per nucleon. Particles of energy less than 100 Mev per nucleon 
almost invariably ended in the stack and their energy was determined 
accurately by their range. Table 1 shows the pertinent values with 


Table 1. Comparison of Predicted and Actual Energies for Alphas 


Energy pre- | 7 ‘ 

Energy from | dicted from Teied Genie °/ Error in | % Error in 
ens range ionization scattering energy of energy of 
number (mev/nucleon) eerie (ev/nucleon) column 3 column 4 


respect to those followed particles which ended in the emulsions and for 
which the energy could be determined precisely. The mean error on the 
energy determined from the raw scattering measurements on these slow 
particles was about 33%. This value was approximately that expected 
from the number of available cells. The mean error in energy determined 
from the ionization measurements for these particles was only 11°, some- 
what better than the estimate based on the fluctuations observed on the 
fast tracks. The estimate derived from the ionization was used for the 
determination of the energy in all cases in which the particle did not end. 
The scattering measurements were merely used in conjunction with those 
of ionization in order to allow identification of the particles (see fig. 4). 
The flux values obtained must be corrected for the absorption of the 
alpha-particles in the air above the plates and for secondary alpha-particles 
produced intheair. The method of correction is described in the Appendix. 
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Table 2 gives the differential and integral fluxes together with the 
magnitude of the correction for secondary alpha-particles produced in the 
air above the plates. The fluxes were measured under an average of 
12 g/cm? of air. The energy listed in the table is the energy of the alpha- 
particles at the top of the atmosphere. 
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The differential and integral energy spectra at the top of the atmosphere 
observed over Saskatoon and Minnesota. 


(a) The differential energy spectrum over Saskatoon. The crosses show the 
effects of the correction for the secondary production of alpha-particles. 

(b) The integral energy spectrum over Saskatoon. Also shown are the results 
of Waddington (1956 b) over Southern England (1) and those of Linsley 
(1955) over America (2). 

(c) The differential energy spectrum over Minnesota. The open circles show 
the results obtained by Waddington (1954) in this stack, while the closed 
circles are those obtained in the present analysis. The crosses have the 
same significance as on (@). 

(d) The integral energy spectrum over Minnesota from the data of Waddington 
(1954) and the present analysis. 


Figure 5(a) shows the differential energy spectrum determined from the 
densitometer energy measurements with the energy of the alpha-particles 
corrected to the top of the atmosphere. The flux corrected for secondary 
production of alphas is also shown. This spectrum shows a definite 
maximum at 300 Mev per nucleon, but also indicates the presence of 
primary alpha-particles with energies down to the average air cut-off 
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value of 124mev per nucleon. The flux of alpha-particles between 
100-300 Mev per nucleon is of the order of 10 times too high to be explained 
by fragmentation of heavier particles, and cannot be ascribed to inelastic 
collisions of primary alpha-particles, or to alpha-particles from other 
nuclear disintegrations. The presence of this maximum would give rise 
to a ‘knee’ in the integral spectrum at an equivalent geomagnetic 
latitude of 57°-58°, but, as noted above, the ‘ knee ’ is not due to a sharp 
cut-off but is due instead to a maximum in the alpha-particle differential 
energy spectrum. Figure 5 (b) is an integral energy spectrum of the alpha- 
particles. The data plotted are those determinations in which the energy 
was determined directly. 


PARTE LE 
§ 3. INTRODUCTION 


Recent experiments have shown that the cut-off energies of cosmic ray 
particles cannot be calculated in the conventional manner from geo- 
magnetic theory. In particular the experiments of Simpson et al. (1956) 
and Waddington (1956 b), first reported at the Mexico Conference, have 
shown respectively that the geomagnetic equator does not coincide with 
the positions of minimum cosmic ray intensity, and that the cut-off energies 
observed over England are considerably higher than those predicted by 
geomagnetic theory. 

As a result it is clearly desirable that direct experimental determination 
of the true cut-off energies should be made at a number of points on the 
earth’s surface. Such determination are particularly important over the 
American continent since the majority of cosmic ray investigations have 
been made there. 

For this reason we have re-examined the experimental data on the 
multiply charged particles obtained by Dainton, Fowler and Kent (1952), 
in astack of glass-backed emulsions flown on October 5, 1950 at Minneapolis ; 
nominal geomagnetic latitude 55°N. This stack was flown at a mean 
altitude of 96 000 feet, and the mean amount of matter above the emulsions 
was 14 g/em?. Asaresult the alpha-particle air cut-off energy was 130 Mev 
per nucleon. The time latitude curve is shown in fig. 6. 

At this geomagnetic latitude and longitude the cut-off energy for 
multiply charged particles incident vertically should, assuming a central 
dipole field, be 300 mev per nucleon. If instead an eccentric dipole field 
displaced towards 154°H is assumed, Jory (1956), then the cut-off energy 
should be 350 mev per nucleon. At these high latitudes the dependence 
of the cut-off energy on the zenith angle should not result in an appreciable 
spread in the experimentally observed cut-off energy. . 

In the above experiment it was noted that an appreciable proportion. of 
the observed particles appeared to have entered the atmosphere with 
energies less than the expected cut-off value. The significance of the 


presence of the slow particles was not fully appreciated, and they were not 


analysed in detail. However, for the reasons outlined above, we have 
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now made a more careful and critical analysis of the low energy alpha- 
particles observed in this experiment, and find evidence that suggests that 
the cut-off energy over Minnesota was equal to or less than the air cut-off 
energy of 130 mev per nucleon. 
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§ 4, Re-ANALYSIS 


In the investigation of Dainton e¢ al. the emulsions were scanned for 
tracks which were more than 6 mm long and which had a 6-ray density 
(NV) greater than that of a relativistic Lithium nucleus (1 per 100). The 
variation of the scattering parameter, «, with the 5-ray density for the 
relevant tracks is shown in fig. 7. There was an appreciable scanning 
bias against tracks having 6-ray densities near the minimum value, and 
as a result we have only considered the slower alpha-particles, with higher 
d-ray densities. The particles accepted were those which when entering 
the emulsions had an estimated residual range of less than 32 mm. 
These residual ranges were estimated (for those particles that did not come 
to rest) by determining the residual range corresponding to the measured 
values of « and NV, by referring to the position of the corresponding point 
on fig. 7, and then allowing for the distance between the centre of the 
measured segment and the point of entry. This procedure represents a 
convenient approximation to a more rigorous estimate of the range. 

The scattering parameters were obtained from a noise elimination 
procedure between 400 and 200 micron cells, and were all greater than 
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about 0:070°/100u. It seems reasonable to assume that these values 
cannot have been seriously affected by any forms of distortion. That the 
above estimates of the residual ranges are not seriously affected by 
systematic errors, is shown, both by the close agreement between values 
calculated in this manner with those observed directly on particles ending 
in, and by the fact that 34 particles were observed to end as compared 
with the 31 that were expected to do so. 
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Altogether 196 alpha-particles with estimated residual ranges of less 
than 3-2 cm of emulsion were found in this stack. In order to obtain 
values on the differential energy spectrum from this data, we defined 
certain groups of tracks in terms of their geometry, and thus in terms of 
the mean amount of matter they had traversed before being detected. 
Such a division enabled the fluxes to be calculated of those particles 
which entered the top of the atmosphere within a certain defined range of 
energy. Inall four such groups were defined containing 162 oF the alpha- 
particles. The details of the four groups are shown in table 3. In this 
table the successive columns show :— 

(i) the group number ; 

(ii) the number of alpha-particles in te group ; | 

(iii) the effective collecting area in m° steradian sec ; 
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(iv) the mean amount of matter through which the particles have passed 
before entry into the emulsion g ; 

(v) the mean residual of the individual amounts of overlying matter dg ; 

(vi) the mean energy at the top of the atmosphere H, calculated from g 

(vil) the spread in energy AZ, due to the spread in the residual ranges 
accepted, 6-32 mm ; 

(viii) the flux at the top of the atmosphere corrected only for absorption, 
in particles/m?/ster/sec/Mev ; 

(ix) and (x) the percentage corrections due to fragmentation of heavier 
particles and due to the stripping of alpha-particles, Appendix (a) 
and (bd) ; 

(xi) the fully corrected flux at the top of the atmosphere. 

In order to extrapolate the observed fluxes to the top of the atmosphere 
we have taken a mean free path of alpha-particles for absorption in air of 
45g/cm?, Bradt and Peters (1950), Webber and McDonald (1955). The effect 
of the secondary production of alpha-particles has been considered 
separately, and the derivations of the corrections are shown in the Appendix. 
The corrected differential flux values derived in this analysis are shown in 
fig. 5 (c) together with those derived by Waddington (1954, 1956 b) in the 
same emulsions. Also shown in this fig. 5(d), is the integral spectrum 
obtained in these emulsions from this analysis and from the data of 
Waddington, for comparison with that obtained at Saskatoon. 
Waddington’s data was obtained from a scan made for alpha-particles of 
all energies and is directly comparable with the scan made for alpha- 
particles at Saskatoon. 

It should be noted that although the differential spectrum obtained 
over Minnesota does not appear to rise to a pronounced maximum similar 
to that observed over Saskatoon, the two distributions are not significantly 
different. 


§ 5. CONCLUSIONS 


From these results we conclude that : 

(1) The differential energy spectra of alpha-particles at Minneapolis in 
1950 and at Saskatoon in 1954 were closely similar down to the lowest 
observable energies, in spite of the fact that they were obtained at very 
different geomagnetic latitude. In both experiments the spectra show a 
maximum at about 300 mev per nucleon and then fall gradually to about 
130 mev per nucleon. Therefore, since Saskatoon is considerably nearer 
the north magnetic pole than Minneapolis, the cut-off energy over 
Minneapolis must have been less than 130 Mev per nucleon. 

(2) Because of the similarity of these two spectra, it appears that the 
spectra must represent the distribution in energy of the cosmic ray alpha- 
particles before they arrive in the vicinity of the earth, as suggested by 
Ney at Mexico. 

(3) This reduced cut-off energy implies that the latitude from which 
the cut-off energy should be calculated, the * cosmic ray é latitude, was at 
Minnesota, at least three degrees greater than the geomagnetic latitude, 
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(4) Any rise in the total ionization count greater than a few percent 
observed north of the cosmic ray latitude of Minnesota, must imply that 
the cosmic ray protons have a different rigidity spectrum to that of the 
alpha-particles. 

(5) The present results are in disagreement with those of Ellis et al. 
(1954) who conclude that their data indicates a sharp cut-off for alpha- 
particles at an energy of 280 mev per nucleon. 

(6) It seems unsatisfactory to postulate that the spectrum in the low 
energy region is cut-off sharply by a magnetic field in space, and then the 
cut-off is ‘ smeared ’ by traversal of matter. Since calculations show that 
the residual atmosphere and packing above the emulsions in these experi- 
ments was not able to produce any appreciable number of low energy 
particles, it would be necessary to invoke a large amount of absorber in 
space to smear out a sharp cut-off sufficiently to give the observed spectra. 
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Corrections for Helium Nucler Produced by Secondary Processes 


There appear to be three processes by which an appreciable number of 
secondary helium nuclei could be produced. 

(a) The production of secondary alpha-particles by the fragmentation 
of multiply charged particles in the matter above the emulsions. 

In order to estimate the number of alpha-particles produced in this 
manner we have two assumptions. Firstly, that in an interaction which 
results in the fragmentation of an incident multiply charged particle the 
velocity of the fragments will be the same as that of the primary particle. 
Secondly, that the mean free paths and fragmentation probabilities are 
independent of energy. Both these assumptions appear to be reasonable. 

If we consider a detector placed y g/cm? from the top of the atmosphere, 
then the number of secondary alpha-particles J, arriving at that detector 
from a layer of the atmosphere of thickness dx and distance a from the top 
of the atmosphere, due to the interactions of a parallel beam of multiply 
charged particles with charge Z, is given by : 


dl=J(Z, ee exp (—a/A,). exp —(y—a)/A,.da . . . (1) 
Z 


where P,, is the mean number of alpha-particles produced by a particle 
of charge Z in an interaction ; \, and , are the absorption mean free paths 
of the heavy particle and the alpha-particle respectively ; A, is the inter- 
action mean free path of the heavy particle ; and J (Z, R) is the number 
of particles of charge Z with a residual range, at the top of the atmosphere, 
betweeen Rk and R+ AR. Her Rand R+ AR are residual ranges such that 
the alpha-particles which result from interactions at a depth x arrive at 
the emulsions (depth y) with residual ranges of 6 and 32 mm respectively, 
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AR is therefore independent of both « and y for a given charge, and is 
equal to (32—6) . 2/Z mm of emulsion. 

At high energies, where R is effectively independent of x, eqn. (1) may 
be integrated over all x to give the usual relation : 

I=J(Z, R) AA ly te [exp (—y/A,.)—exp (ura) | C Hee (2) 
VA IN EAS A 

At low energies J (Z, R) is a complicated function of x, and it is convenient 
to integrate eqn. (1) numerically. In order to find the appropriate values 
of J(Z, R) we have constructed from the differential energy spectrum of 
the alpha-particles shown in fig. 5, the differential range spectrum at the 
top of the atmosphere for the various charge components considered. We 
have assumed that the spectra of all these particles have the same variation 
with energy per nucleon. 
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These differential range spectra are shown in fig. 8 for groups of 
particles whose charges are centred about those of carbon, magnesium ang 
calcium, and whose fluxes (at Minnesota) were taken as 20, 3 and 2 
particles/m2/ster/sec respectively. The curve shown for the alpha-particles 
corresponds to 320 alpha-particles/m*/ster/sec. 
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The results of these corrections are shown in table 3. It can be seen 
that the corrections are all less than 20°%. As a result it has not been con- 
sidered necessary to modify the curves shown in fig. 8 for the effects of 
secondary production. 

(b) The production of He? particles in stripping collisions of He*-nuclei. 

We have made a similar calculation to the one above to estimate the 
correction that will have to be applied for this process. In making this 
calculation we have assumed that 10°, of all interactions in air will be of 
the He‘ stripping to He* type. This figure is based on an observed figure 
of 6-3°/,(44) in emulsion, obtained from the combined data of Appa Rao 
et al. (1956) and Waddington (1956 a). 

In this calculation the value of 4R will have changed, since instead of 
accepting a He*nucleus with residual range between 6 and 32mm, a 
He?-nucleus will be accepted if its residual range lies beween 6 and 36 mm. 
Such He?-nuclei will be produced by a He*-nuclei with a 4R of 40 mm. 

Using the alpha-particle curve given in fig. 8 we find that the correction 
due to this stripping process amounts to about 2% in each of the groups 
defined in this experiment. 

(c) The production of slow helium nuclei by singly and doubly charged 
nuclei in interactions. 

It has been shown by Yasin (1955) that helium nuclei with energies 
greater than 80 Mev (not evaporation particles) emitted from disinte- 
grations initiated by singly or doubly charged particles are predominantly 
He?-nuclei. ‘This conclusion is supported by inspecting the relative 
numbers of He*- and He*-nuclei—identified from their residual ranges— 
observed in emulsions exposed to the K+ beam of the Bevatron. It is 
therefore theoretically possible to distinguish the He-nuclei produced in 
this manner from those of primary cosmic ray origin, which are presumably 
mostly He*-nuclei. However, in specific cases it is rarely possible to 
distinguish between He®- and He*-nuclei ending in emulsions, and it was 
therefore necessary to calculate the rate of production of these He-nuclei 
from other experimental data. 

We have used the results of Sorensen (1952) and of Dahanayake et al. 
(unpublished, but see Dahanayake et al. 1955). Sorensen considered He- 
nuclei with residual ranges greater than 0-7 mm emitted from stars 
found in emulsions exposed at 68 000 feet ; while Dahanayake et al. took 
a sample of those particles with residual ranges greater than 2 mm emitted 
from stars having two or more shower particles (n,>2) found in emulsions 
exposed at 85000 feet. The energy distributions of the observed He- 
nuclei in the two experiments have been shown by Yasin (1955) to be 
closely similar, and to follow the relation : 

const. 
N(£)dE= aT dH. eM” oe eS) 

Dahanayake et al. observed 90 He-nuclei, of which 32 hadresidualranges 
greater than 6mm, emitted from 4742 stars having n,>2. When a 
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correction was made for the geometrical selection used, it was found that 
about 150 He-nuclei with residual ranges greater than 6mm had been 
emitted from these stars. In the emulsions used in the present experi- 
ment the density of stars with n,>2 was 73 stars per cm? of emulsion. 
Thus the production rate, for this type of star only, was 2:3 He-nuclei 
per cm’ of emulsion. 

Unfortunately, it is not possible to find the total production rate of He- 
nuclei from Sorensen’s results, since the total number of stars examined 
was not specified. Hence we have had to use the production rate derived 
from the results of Dahanayake et al. and correct for the production of 
He-nuclei by stars withn,<2. To make this correction we have examined 
the observed star distributions in the emulsions used by Sorensen, as given 
by Camerini et a/. (1952), and that in the emulsions used in the present 
experiment. We find that they are not appreciably different. From 
Sorensen’s data it appears that for each star with n,>2 that emits a He- 
nucleus there are 2-8 with n,<2. Thus the total production rate in 
these emulsions of He-nuclei with residual ranges greater than 6 mm is 
8-7 He-nuclei per cm® of emulsion. 

In order to obtain an estimate of the total production rate in air from 
that in emulsion we have considered the numbers of black and grey 
prongs \,, emitted from the helium producing stars. In emulsion the 
elements carbon, nitrogen and oxygen (CNO) make up 25°% of the total 
geometrical cross section, but only 15°, of the total interaction cross 
section (Lock et al. 1955, for 950 Mev protons). From the data of 
Camerini et al. 75°, of all stars have NV, <7, and so simulate stars in light 
elementst. Thus, using the experimentally observed interaction cross 
section, only one-fifth of those stars with N,<7 are produced in light 
elements. Sorensen’s data shows that only 30° of the He-nuclei come 
from stars with V,<7. Therefore, assuming that in small stars light 
and heavy nuceli are equally capable of producing He-nuclei, the 
production from CNO in emulsion was 0-52 He-nuclei per cm? of 
emulsion, or, 0:75 He-nuclei per g of CNO, since CNO have a density 
of 0-7 g/cm? of emulsion. 

In work involving emulsions, there is generally a scanning bias against 
the detection of stars with small V,. But, since Sorensen’s results show 
that the efficiency of production of He-nuclei falls off sharply as NV, 
decreases, and, since a star with a 6mm He-nucleus is less likely to be 
missed than a typical small N,, star, we have assumed that for these 
purposes the effects of any scanning loss of small stars may be neglected. 

In calculating the production rate in CNO, we have assumed. that in 
small stars light nuclei are equally as likely to produce He-nuclei as are 
silver and bromine nuclei. This assumption appears improbable, and we 
therefore believe that the production rate for light nuclei given above may 
well be an overestimate of the true value. 


+ Due to the scanning bias against interactions with hydrogen, we have 
assumed that all stars in light elements will be in CNO. 
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In calculating the number of He-nuclei produced above the emulsions, 
we have used the production rate given previously. In addition, 
Sorensen has shown that the emission of He-nuclei is not isotropic, the 
median angle of the angular distribution about the primary being about 
50°, whereas the primary particles are nearly isotropic over the upper 
hemisphere. For this reason we have taken the solid angle that He- 
nuclei are emitted into as being only 37 steradians. 

Making the above assumptions, we find that in none of the groups of 
particles defined in the present experiment are more than 1 or 2% of the 
observed alpha-particles attributable to production of He-nuclei by singlv 
or doubly charged particles. It appears, therefore, that this mode of 
production may be neglected. 

That this conclusion is correct may be justified experimentally. 
Waddington (1956 b) working on a stack of emulsions exposed at a cut- 
off energy of about 650 mev per nucleon, and an altitude of 106 000 feet, | 
searched an area containing an estimated 230 fast alpha-particles, and 
found only one with a residual range less than 32mm. Dainton et al. 
(1952) observed in emulsions exposed at a similar cut-off energy but at 
60 000 feet, 9 slow alpha-particles in an area containing 70 heavy primary 
particles. These 70 particles are estimated to correspond to 1400 fast alpha- 
particles. The few slow alpha-particles observed in these experiments 
can all be attributed to the fragmentation of heavy particles rather than 
to production by the interaction of singly or doubly charged particles. 


ACKNOWLEDGMENTS 


We wish to thank R. B. Thorness and Ray Maas for their help in the 
design and construction of the Saskatoon equipment ; W. Huch for his 
help with the balloon flights; Dawn Copeland, Pamela Curry, Biruta 
Sommers, Janet Witchell and Leona Wurm for doing the very difficult 
and tedious scanning involved in this experiment. 

_ We are grateful to Professor C. F. Powell for the hospitality of his 
laboratory. One of us, C. J. W., wishes to thank the Royal Society for 
the award of a Mackinnon Research Studentship. 


REFERENCES 


AppER Rao, M. V. K., Dantzt, R. R., and NreiaKkontan, K. A., 1956, Proc. Ind. 
Acad., 48, 181. 

Branvt, H. L., and Peters, B., 1950, Phys. Rev., 80, 943. 

CaMERINI, U., Daviss, J. H., Fowumr, P. H., Franzinertt, C., MurruEap, H., 
Te. W. O., Perkins, D. H., and Yexurrett, G., 1951, Phil. Mag., 42, 

1. 

DAHANAYAKE, C., Francois, P. E., Fusmmoto, Y., IREpDALE, P.,.WappiNe@Ton, 
C. J., and Yasin, M., 1955, Nuovo Cim., Ser. X, 1, 888. 

Datnton, A. D., Fowier, P. H., and Kent, D. W., 1952, Phil. Mag., 48, 729. 

Euis, R. A., Gortiims, M. B., and Van Auten, J. A., 1954, Phys. Rev., 95, 
304(a) 

Jory, F.S., 1956, Phys. Rev., 102, 1167. 

Linsey, J., 1955, Phys. Rev., 97, 1292. 


Cosmic Ray Spectrum of Alpha-Particles 175 


Lock, W. O., Marcu, P. V., MurrHEap, H., and Rossmr, W. G. V., 1955, Proc. 
Roy. Soc. A, 230, 215. 

Neuer, H. V., and Stern, E. A., 1955, Phys. Rev., 98, 845. 

Ney, E. P., 1955, reported at Mexico Conference. 

Smeson, J. A., Fenton, K. B., Katzman, J., and Ross, D. C., 1956, Phys. Rev., 
102, 1648. 

SorEnssEn, 8. O., 1951, Thesis, University of Oslo (see also Phil. Mag., 42, 188). 

WanppineTon, C. J., 1954, Phil. Mag., 45, 1312 ; 1956 a, Ibid., 1, 105 ; 1956 b, 
Nuovo Cim., 3, 930. 

WEBBER, W. R., and McDonaxp, F. B., 1955, Phys. Rev., 100, 1460. 

Yasin, M., 1955, PA.D. Thesis, University of Bristol. 


pe LO | 


A Simple Treatment of Meson—Nucleon Scattering { 


By 8S. F. Epwarps and P. T. MatrHEws 
Department of Mathematical Physics, University of Birmingham 


[Received April 26, 1956, Revised August 10, 1956] 


ABSTRACT 


It is shown that the p-wave meson-nucleon scattering cross sections 
can be very simply explained by calculating the first order effective 
potential between meson and nucleon and the scattering in this potential. 
Some comments are made on the S-wave problem and virtual heavy 
meson effects. 


§ 1. InTRODUCTION 


Tue oldest and simplest application of field theory to the low energy 
nucleon-nucleon problem is the adiabatic approximation. An effective 
potential (the Yukawa potential) is calculated by Born approximation. 
This is then inserted in the non-relativistic Schrédinger equation for the 
relative motion of the two nucleons, which is then solved exactly for 
the binding energy and scattering. We show here that the qualitative 
features of the successful description of the p-wave meson—nucleon 
scattering, given by other authors on the basis of the no-recoil pseudo- 
vector theory (Sartori and Wataghin (1954), Chew and Low (1956)), 
may be derived in this, conceptually extremely simple, way. The 
effective potential for the meson is determined by Born approximation. 
This is then inserted in the relativistic Schrodinger (i.e. Klein—Gordon) 
equation for the meson, and the scattering in this potential calculated 
exactly. We shall refer to this as the potential approximation. 

This procedure does not preserve the crossing symmetry. However, 
if the equation is symmetrized with respect to crossing and solved as 
an expansion in the energy one reproduces the Pocalte of the above 
authors. 

If the interaction is supplemented by the ¢? term in the non-relativistic 
approximation to y; theory, there is an S-wave scattering, which is given 
exactly by the potential approximation. Taking the same coupling 
constant and cut-off as in the p-wave case, one obtains an S-wave of the 
required order of magnitude, but lacking the observed dependence on 
isotopic spin. 

In the final section the possibility is investigated that this isotopic 
spin dependence is a consequence of virtual K-meson effects. It appears 
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that these may well give a non-negligible contribution, but do not appear 
likely to give a simple explanation of the effect. 

The equation from which we start is derived in closed form in the 
Appendix, which thus provides a more formal basis for the outlook 
adopted in § 2. 

It is perhaps worth remarking that the Tamm—Dancoff approximation 
is not equivalent to the scattering in an effective potential, although 
it has been described in these terms (see for example, Gell-Mann and 
Watson 1954). The Tamm—Dancoff method excludes the effect of virtual 
meson pair creation, which may take place in any potential. (This is 
discussed in more detail by Feldman and Matthews 1956.) This 
approximation replaces an algebraic equation by an integral equation, 
which has been solved approximately by Chew (1953). 


§ 2. THE p-WavE SCATTERING 
The scattering at energy w, (w?=p*+p?), is determined by the equation 


(p?-+p?—w*)(p | ¢,)—V,(P» Q)(q]¢)=0. . . . (2.1) 


Integration or summation over repeated variables is implied. The 
potential V is determined by Born approximation. If the interaction 
energy density is 


H (x)= —)j,(x)¢,(u)=—g0(x)T,0. Vd,(v), . . . (2.2) 
then ; 
VAP» 9)=(27)°19(P), J(—Q)/m, - - - - ~ (28) 
where 
PADIS DIT ee) ae a (A) 


This is obtained by perturbation theory from the familiar ‘ crossed ’ 
and ‘ uncrossed ’ graphs. The cross section is proportional to the square 
of the modulus of the ¢-matrix, where the t-matrix satisfiest 


tra(Ps 4)=Vire(P» 4) + Ve (Ps k)(h? +p? —w? te) *tye(Ky q). (2-5) 
Introduce t,,“(w) and V,,*9(w), defined by 
tre(Ps 4) =Palirs* (@) Ips 
VelP> V=PaV rs (w)dgp  - - - - + + (2.6) 


ns er 
+ This equation is proved thus : if di” is the incoming wave function, satis- 

‘fying (p?-+2—w?) din—0, then by (2.1), ¢ satisfies 

F d=90+ (p?+w—w? ie) Vo. ai atampee ia (C) 
The differential cross section, o, is given by o~|¢ |?, where 

p=din+ (p?+> w?—te)“1 t dit. 
Thus téit=V¢d. Hence, multiplying (a) by V, 
t=V-EV (p?+p?—w® te). 
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"These matrices then satisfy the matrix equation (not an integral equation) 
# OSV CPN Vd Pe Ae ee me oe 

where d@”(w) is the diagonal matrix 

7 1 k*d?k 
Os =) (k?+-p2—w*—te) ” 
The matrices ¢ and V may be diagonalized by projection operators. 
For isotopic spin the projection operators for the 4 and % states are 
respectively 


(2.8) 


Pi ys=(1/3)7,755 

P= Ors— (1/3) 7,75: ee C22) 
and similarly for total angular momentum. If P,, («=J7'), are the 
projection operators for states of angular momentum and isotopic spin 
4J and 47 (i.e. a=11, 13, 31, 33), 


(w= > Pi lale ss ke eee 

and, expressing (2.3) in terms of (2.9), 
V=(9?/(27)8n) SnPy + - - + + (211) 
where 21,;=—8, %43=%313=—2, N3,—4. Substituting mto (2.7), and 


taking matrix elements for states of definite 7’ and J, gives the algebraic 
equation 
Me Gy 
ae w—gn, fl" (w)| (20) 
The factor d”(w) is cubically divergent. Mass renormalization consists 
of replacing it by 


(ts) 


cath: (an) 0. Oa ne (2.14) 
which gives t(w) a pole at w=0, corresponding to the neutron. This 
approach agrees with the familiar renormalization of perturbation 
theory. If d is cut off at A, it can be simply evaluated : 


A(w)= (27/3) {2Aw?—7(p3—ip3)}. . . . . (2.15) 
Finally in terms of the phase shifts 
L/t,(w)~p* exp [—76,]/sind,.. . . . . (2.16) 


The constant of proportionality can be determined by equating the 
imaginary parts of (2.13) and (2.16). Then equating real parts one has 


prootd, , 3 (4x 2Aw Bp 
= his cae ee —"|. Mae (80) 


Making the approximation of dropping the final term, one can write 
this in the effective range form 


Ayp* cot 0, /o==1—-wlc),, “epee ace ee loans) 

where Teagan 
Ag=+ 3(£) 5/0 bs) gegeaeeaee Rep ehepa mean elo) 
and oOo— 6777/0 g2A. oe 8 eS ° (2.20) 


A resonance occurs at w=w, if n, is positive, that is for «—33. 
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§ 3. UnITaRITY AND CROSSING 


The approximation of § 2 does not preserve the crossing symmetry of 
t,,(@). Tt we interpret V,, of (2.7) as the exact potential (derived in the 
Appendix), the equation is exact and is unchanged if we write 


tralo)=V (eo) +31 .y(w)d"t,.(@) + V(—cw)d"t,.(—a)]. . (3.1) 
The imaginary part of ¢~1 is fully determined by unitary 
Imi-t==— 4/347, 2 we we > (3.2) 


To obtain the effective range formula, which automatically satisfies 
crossing symmetry, we expand (3.1) for the real part of t1(w)/w as a 
power series in w, approximating for V,, as in §2. (Note that we are 
still handling an algebraic, not an integral equation.) This yields 


2R iy —-1 5 : 
a ae, oth rs a {" = “es FA EO) 
where A” is the numerical matrix determined by 
boi@) j= (a) ee fee oe USL) 
This reproduces (2.18) with the definition of w, modified to 
ee OTC GAP serene, A a) (9.0) 


where C,=(—6, 0, 6) in agreement with the above authors. The 
experimental data for the 33 wave are fitted by 
4 =0-08 2, 
A= ps. AEP Area 4a Pay te AO) 


The other phases are then small and negative. 


§ 4. Toe S-Wave 
If the non-relativistic interaction is regarded as an approximation to 
the relativistic ps—ps theory, there is an additional term in the Lagrangian 
(g defined as before) of 
Le nite oc Ta ONE d ie SCe) 
This gives rise to an S-wave scattering, which can be calculated exactly 
by the ‘ potential ’ method : 


—2q2m . 
t,(w)= J , yh one eeeattl aan) 


1+ 29?md(w)/(27)8 
where 


Heyes [ pot ie ak, ede (0 


d(w) is linearly divergent, and it is consistent with our approximation to 
simulate the renormalization of the relativistic theory with a cut-off A. 
This leads to the expression for the S-phase shift 


(4.4) 
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Taking g? and A from the p-wave data, this gives an S-wave scattering 
of the observed magnitude, but the strong dependence on isotopic spin 
observed experimentally is completely lacking in this model. 

This discrepancy might be corrected in a relativistic scheme, but 
another possibility arises in the contribution from virtual K-mesons. 
Since 7, @ and = mesons are all strongly coupled to nucleons, they must 
be strongly coupled to each other. For simplicity we assume there is 
a direct local coupling : . 

| Ly=mg yr (FP + Bay, ee (4.8) 
where y and 9 are the 7 and @ meson fields, respectively. We also 
assume 

Lo=g(WAs+Aps® . 2... (4.6) 
and 
La=Gar?, (7Z,9% + Dahl9%). 2 2... (4,7) 
consistent with the Gell-Mann—Nishijima scheme. The coupling constants 
g; appearing in these interactions are pure numbers and we take g,?/47—1, 
which appears to be consistent with the observed production rates. 
A rough calculation of the resulting effective meson—nucleon potentials 
of order 91°," and ¢,7g,” indicates that these effects might give corrections 
to meson—nucleon scattering of the order of 10°. However, the factor 
m included in L, for dimensional reasons, might well be replaced by the 
K-meson or 7-meson mass, so 10° should be regarded as an upper limit. 
Of this correction, crossing symmetry requires that the 7t-dependent 
part is further reduced by a factor of (w/K-meson mass). Thus although 
heavy mesons may play a part in a detailed theory of 7-meson scattering, 
they do not offer a simple answer to the 7-dependence of the S-wave. . 


A PAP ESN: DEX 
We derive here in closed form the equation for the meson propagator, 
4’, in the presence of the source j,(2). 
Consider the Lagrangian 
L=—3¢,(—-D+ Pv be +(jt3 be 2 - - . (Al) 
where J(x) is a classical field. The propagator in Schwinger’s (1951) 
notation is 
A’,(x, y= ((f,(@), b(Y))4)— (d(x) )<bs(Y) D- « st 2) 


Using the equation of motion for ¢ 


(pie = (jr Ps) ye Gi ee 


Sa »/8I,—18 55. siete: , sarens) 
But OIE _ 4dr) 344.) 
Wim © ERs ar” 
EW yy | 
= chy Ae EE PoE CAL ES 
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Thus 
[((—L1++#?)8 1-89, 0/8 ($14 p= — 78,» =. (A) 
Now 
8 /94bi)=t (Jr Je 2—bCr CG), «=». (A 8) 
where, for example, 
P $j) jy) exp [tf Lda*]3d 
CG AX), IY Ee (aed) 


the P denoting that the operators j,(~) must operate in their chronological 
order. Expanding (A 6) in powers of j gives, to order j?. 
d¢j,(z) /9<b(Y) = —t1j(X), JY) O(ty—te)s ° (A 8) 


where 


Substituting this into (A 5), and making a Fourier transform, leads 
directly to eqns. (2.1)—(2.3) for the meson wave function. 
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SUMMARY 

It is shown that, for a shock wave advancing through any physically 
real barotropic medium, the nature of the reflected wave is uniquely 
determined by the relative shock impedances of the medium through 
which the incident wave passes and the medium upon which it falls. A 
simple criterion is found for determining the nature of the reflection of a 
detonation wave at the end of a block of explosive. Particular attention 
is paid to the examples of an explosive in contact with a gas, water or 
solid surface. Certain approximations, which may be useful when there 
is a reflected shock wave, are discussed. 


§ 1. INTRODUCTION 

As far as the author is aware the only published information on the 
reflection and transmission of shock waves at boundaries is contained in 
papers by Paterson (1948), Taub (1947) and in studies made at the Institute 
of Aerophysics in the University of Toronto (Bitondo et a/. 1950, Bitondo 
1950, and Ford and Glass 1954). These are all concerned with shock 
waves advancing through ideal gas and meeting another ideal gas. It is 
shown below that for a shock wave-advancing through any physically 
real barotropic medium, the nature of the reflected wave is uniquely 
determined by the relative shock impedances of the medium through which 
the incident wave passes and the medium upon which it falls. A simple 
criterion is found for determining the nature of the reflection of a deton- 
ation wave at the end of a block of explosive. Particular attention is 
paid to the cases where the explosive is in contact with a gas, water or a 
solid surface. ‘The paper ends with a discussion of approximations which 
may be useful when there is a reflected shock wave. 


§ 2. Equations oF Motion For THE LNcIDENT Wave 
If the incident shock wave and the boundary upon which it falls were 
plane and parallel and were contained in an infinite rigid unbreakable 
tube, and if viscosity, heat-conduction and all boundary effects were 
absent, the motion in each phase would be one-dimensional. The assump- 
tion of one-dimensional flow will be made in the following paragraphs, 
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for, although one-dimensional flow cannot be exactly realized in practice 
(and in fact, there is a steep decay of pressure behind a detonation wave) 
the conclusions drawn for the interfaces will be valid at least for the 
initial instant in which the reflection occurs, over an area of the surface of 
contact which will depend upon the area over which the boundary may be 
assumed to be plane. Thus the equations will indicate the nature of the 
phenomena which occur and also the initial magnitudes of the physical 
quantities. Suppose that the velocity of the incident waves is Vj. 
Suppose that this wave is moving through a medium in which the equili- 
brium pressure and density are respectively py and py. Suppose that 
these values rise instantaneously to p, and p, on the passage of the 
incident wave. Let uw, and c, be the particle velocity and local speed of 
sound behind the incident wave. The conservation of mass requires that 


Per spk Vi—u}) . . . . . . . 5 (1) 
and the conservation of momentum is expressed by the equation 
P1—Po=PoV 11 = popi'1"| (P1— Po) ites eS) 


If the incident wave is a detonation wave, the Chapman—Jouget condition 
requires that 
VW yricGase cea) Bile oa) a SENS) 


and if it may be supposed that there is an effective ratio of specific heats y 
for the products of detonation in the neighbourhood of the detonation 
wave, then 

Be a Dy ce et Te gee ee ie CE) 
In the detonation phenomena p, may be neglected in comparison with py. 
Tt then follows easily from eqns. (1)—(4) that 


+1 
ES Po 
ah re Se Ban Pm C5) 
2 Py eer 


i : 
(y+ ])po ¥P1 


§ 3. Equations oF Morion AFTER REFLECTION OF THE WAVE 

That a shock wave is transmitted to the medium ahead of the incident 
wave is physically obvious and will not be discussed theoretically. The 
nature of the reflected wave has to be determined. 

Suppose that the shock wave transmitted to the medium upon which 
the incident wave falls has velocity V,. Let 7, be the density behind it 
and 7, the (equilibrium) density ahead of it. Let wy be the beat 
imparted to the medium by this shock wave and p, the pressure of the 
wave. Conservation of mass requires that 


%oV a=Tal V g—Us) OS) aan A. ome) 
and conservation of momentum that 


Pe Wary lg (7) 
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Daa Kol sls | ae fols Mes a) ee ee es 


In the infinitesimal layer of the first medium immediately in contact with 
the second the pressure and velocity must be p, and w, respectively, 
equal to the pressure and velocity set up in the second medium by the — 
transmitted wave. Let p, be the density in this layer. The boundary is 
an example of a contact discontinuity. 

First suppose that the reflected wave is a shock wave moving with 
velocity W, in the direction opposite to that of the transmitted shock wave. 
The motion (supposed steady) in the first medium immediately behind 
this shock wave is then governed by the equations 


po(Uot W )=pi(ui+ We) » Sees oe RED 
and 
Pam sa Pipe (Uy —Ug)?=po(Uet We)(uy—Ue) . - (10) 
[Pa (Or 


representing the laws of conservation of mass and momentum, and the 
equation of conservation of energy, which in the particular case of a gas 
with effective ratio of specific heats y takes the form 


ce = a as 
PAP Geely Gaalag? fs 


where z,=p,/p,. Writew,+ W,—X,; then X, is the speed of the reflected 
shock wave relative to the medium downstream of it and is necessarily 
positive. It is easy to show from (2), (7), (9) and (10) that 


(Po—P1)/(P1— Po) =p 2X al V s—poV 4) /poV (70 V 2+ po Xo): 


Hence a reflected shock wave, requiring p,>,, requires also that 
TV 2>poVi. This result is not dependent upon the form of the energy 
equation nor has p, been neglected in deriving it : it is therefore true for 
any media and for any type of incident shock wave. 

Next consider the case when a rarefaction wave is reflected at the 
contact discontinuity. The only difference in the equations occurs from 
the change of conditions to be satisfied behind this wave. There is no 
longer a discontinuity in the pressure in the first medium after the instant 
of impact of the incident shock wave upon the second medium, but only 
a discontinuity in the derivative of the pressure at the head of the wave 
as it moves upstream. The eqns. (9) and (10) therefore cease to apply. 
The condition for a receding rarefaction wave in any barotropic medium, 
i.e. one in which the pressure depends only upon density, is (see, for 
example, Courant and Friedrichs 1948), 


wy=u+| Ca, =. i vs ee Le) 
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where c?=dp/dp. On substitution for U, from preceding equations it 
follows that 


(oV 2 poV s1=Pa—P1—7oV a| “oe D 
Ps 
showing that if p is a monotonic increasing function of p (as it is for 
any physically real barotropic medium) a reflected rarefaction wave 
(P2</) is only possible if 7)V,<p)V 4. 

The only condition on the media required by the above proof is that 
the incident shock wave should be passing through a barotropic medium. 
The nature of the reflection is thus determined completely by the relative 
values of the shock impedance of the two media, the shock impedance of a 
medium being defined as the product of the equilibrium density of the 
medium and the speed of the shock wave moving downstream through it. 

Expressed in this way, the criterion for a reflected shock or rarefaction 
wave is very simple, but it requires knowledge of the speed of the trans- 
mitted shock wave, which will not be known until the whole reflection 
and transmission phenomenon has already been determined. In the 
following paragraphs a closer examination of the equations is made for 
an incident detonation wave and there is a discussion of the wave reflected 
after impact upon an ideal gas, water or a solid surface. 


§ 4. APPLICATION TO THE REFLECTION OF A DETONATION WAVE 


When the wave reflected into the detonation products, treated as ideal 
gas with an effective ratio of specific heats y, is a shock wave, 


273 (Z2—1)? =p (uy — U9)" {(y—1) +(y+1)z} 
from (10) and (11). By means of (8) and (5) this may be written in the 
form 
yy ety Vy L204) P= 0. (18) 
where Z?=z», ZQ(Z)=u,/u, and Q(Z)=+ {(y+1)(p9/"7o)(1—10/1"2) }*”?- [ 
depends only upon Z for given media and given values of p,, po.] 
Equation (13) is therefore the equation determining the pressure behind 
the reflected shock wave when the incident wave is a detonation wave. 
When w,, uw, and c in (12) are expressed in terms of po, p, and densities. 
’ the equation for the pressure behind a reflected rarefaction wave is 
obtained. By means of (5) and the assumption as before of the adiabatic 
relationship in the product gases of the detonation, Z is in this case the 
solution of the equation 
g(Z)=ZQ(Z)—14 (2y/(y—I (AP ™M™—-)=0 . . . (14) 
where Q(Z) is defined as above. 
The type of reflection may be determined in general by seeking a root 
of f(Z)=0 with Z>1 for a reflected shock wave or a root of g(Z)=0 with 
Z<1 for a reflected rarefaction wave. 


+ If the boundary is in motion before the shock wave strikes it, the shock 
impedance is defined as the product of the density of the medium into which 
the shock waveismovingand the speed of the shock wave relative to thismedium. 
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Now g'(Z)=d(ZQ)/dZ+-2Z-1/” which is positive for all positive values 
of Z since, as is easily seen, dQ/dZ>0 for given media with fixed pj, Py ; 
but g(0)<0 and g(1)—Q(1)—1. Hence there is a solution with a rare- 
faction wave if, and only if, Q(1)>1. In this case, whenever Z>1, 
ZQ(Z)>1 and hence u,>u,. It follows from (4) that no shock wave 
solution ( Do> /p1) is possible and the solution with a reflected rarefaction 
wave therefore exists and is the unique solution of the problem. 

If Q(1)<1, then ZQ(Z) increases with Z until it reaches the value unity 
at Z=Z,, say. There is no shock wave solution when Z>Z,. However, 
f(1)<0 and f(Z,)>0 and when Q(1)<1 it is clear that f(Z) is a strictly 
monotonic aioreiwe function of Z in the range 1<Z<Z, ; thus there is 
a root, Z=Z, of f(Z)=0 in this range and it is the only physically valid 
root of (13). 

The nature of the reflection is therefore determined uniquely by the 
value of Q(1). The result may be expressed as follows : 

There is a refiected shock wave or rarefaction wave according as 


(+1) (1-2) Soret). i. .8 ees 


where 7, is the density in the second medium corresponding to the 
passage through it of a shock wave with the pressure of the incident 
wave. 

EHxamoples 
(a) Explosive in contact with an ideal gas. 

If the medium ahead of the explosive is an ideal gas with ratio of specific 
heats J’, the increase of density across the transmitted shock wave is given 
by 

1[to= (LE —-1)+-£+1) p/p / (L+1) +L —1) p/p} 
For a shock wave with pressure py>o, 11/fo==([-+1)/([—1). Then 
Q?(1)=2(po/7o) (v-+1)/(£+1). For combinations of common explosives 
and gases Q(1) has a value very much greater than unity, showing that a 
reflected rarefaction wave is the rule. 


(6) Hxplosive in contact with water. 

An equation of state for water given in Cole’s book on underwater 
explosions (Cole 1948) and which has been found by Schall (1950) to be in 
good agreement with experimental work, is the following : 


where p, S, 7’ and r are respectively pressure, specific entropy, temperature 
and density, n is an index which may be given the value 7-15 with good 
accuracy, B(S) is a slowly changing function of the entropy, and the 
density 7 is expressed as a function r(Z’, p) of temperature and pressure 
(p=0 corresponds to atmospheric conditions). For the purpose of com- 
putations of pressures and velocities it is often useful to replace the above 
formula by 

p= BS 1 ae eee eee LG) 
where B is taken as a constant and put equal to 3- 047 creed This 
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form of the equation of state has been used in well known researches by 
Polachek and Seeger (1949) on shock waves in water-like substances and 
in a paper by Finkelstein (1947) on the reflection of a shock wave in water. 

For a shock wave with a pressure P1>150 000 atmospheres, 
7, >1:729 g/em®. Then Q?(1)=0-41p9(y+1)>1 (in general) and it follows 
that the reflected wave is a wave of rarefaction. 

As a matter of interest eqn. (14) may be solved for the case examined 
experimentally by Schall, in which an explosive (PETN) was detonated 
in contact with water. The pressure p, was found experimentally to be 
25 x 10% dynes/cm? and py was 1-6 g/cm’. If y=3 is taken for the explo- 
Sive gases, for want of a better value, then a numerical solution of ( 14) 
gives a value of Z corresponding to p,=16-9x 10!°dynes/em®. The 
value obtained by Schall from his experiments was p,=19-2x 101° 
dynes/em?. The error is approximately 12%, which is not surprisingly 
high in view of the assumptions about the value of y and the equation of 
state of water at these high pressures. In fact, the values of 7 derived 
from Cole’s equation of state are really too high, hence Q? is also too high 
and consequently Z and p, are underestimated. 


(c) Explosive in contact with a solid. 


For the pressures set up by the detonation of an explosive it is not 
possible to treat a solid target as elastic. An equation of state for steel 
and lead suggested some years ago by Pack et al. (1948) has proved success- 
ful for comparison of theory with experiment in several instances, and is 
expressed in such a form that it may be used for any solid for which com- 
pressibility data are available for moderately high pressures. or metals 
in general 7,¢)>p,V, so that if the shock wave transmitted to a metal 
were to move into it with a supersonic velocity V.(>c ), it would follow 
at once from the general result on shock impedance that a shock wave 
would be reflected. However, it was shown in the paper by Pack et al., 
quoted above, that the shock wave velocity in a metal may well be less 
than the velocity of sound. (This was found to be true theoretically and 
in experiments for a steel target.) Thus a general statement cannot be 
made without further investigation, and until more is known about the 
equations of state of solids (including metals). It seems likely, however, 
that a reflected shock wave will be the rule for most if not all combinations 
of explosive and metal, since (y-+1)p)<7 if the metal is not a light one, 
and even for a light metal such as aluminium the factor (1—r/71) is 
strongly damping on the left-hand side of the inequality in (15). 


§ 5. ABLARD’S APPROXIMATION 


The practical importance of calculations on the pressure set up by 
reflection of a shock wave lies in attempts which have been made to relate 
the velocity of the rear surface of a target block, which can be SS se 
without too great difficulty, to the pressure applied to its front. A 
theoretical analysis illustrating the connection for a simple wave system 
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is given in Part II. In the interpretation of experiments in which a 
detonation wave is reflected from a metal block, Ablard’s approximation 
is sometimes used, namely p,X,=p)V,. ‘Since, for a detonation wave, 
po V1=p1¢s from (1) and (3), this implies that p,X,=p,c,. An examination 
of the equations of motion given above for a reflected shock wave (not 
necessarily a detonation wave) shows that 


. peXa=pt| {(— 1)+(y+1)z hey | 


1 1) 
For z,=1+« («<1), p Kepied +e oe et. 


When e=0-315 and y=3, p,X, differs from p,c, by 10%. To judge the 
error in another instance it may be noted that for the experiments on 
steel and lead reported by Pack et al. the calculated values of z, were 1:67 
and 1-59 respectively ; in these instances Ablard’s approximation is 20% 
in error. 

Ablard’s approximation may be compared with the physically more 
understandable approximation that the reflected shock wave moves back 
with the local speed of sound relative to the gas in contact with the solid. 
It may be shown that 


1/2 
X=eal {+14 0a 272 | ) 


which, when z,=1-+« (e<1), gives 


y+] 
Xye,{1— “hy 


et (vat Cee \ 


32y* 

Thus X,—c, has the same percentage error as Ablard’s approximation to 
the first order in «, and is slightly better to the second order in «. Some 
use of this approximation will be made in Part II of this work where, in 
discussion of the wave motion arising from the impact of a shock wave on 
an elastic target, the value of « is small. 
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SUMMARY 


An examination is made of the variation with time of the pressure and 
velocity in an elastic target of finite thickness subject to the impact of 
ashock wave. The cycle of reflections and transmissions at the boundaries 
of the target is described in detail. 


§ 1. INTRODUCTION 


A PROBLEM Of some practical importance concerns the motion of an 
elastic target of finite thickness subject to the impact of a shock wave. 
The system of reflected and transmitted waves which is created by such an 
impact is examined in detail below for a one-dimensional system, and 
formulae are found for the variation with time of the pressure and velocity 
in the target. It is assumed that the shock impedance of the target is 
greater than that of the ‘ outer medium’ through which the shock wave 
approaches the target, so that a shock wave is reflected as well as trans- 
mitted (see Part I, Pack 1956). The outer medium is supposed to be 
barotropic, i.e. in continuous flow it has an equation of state of the form 
p=p(p), where p and p are respectively the pressure and density. The 
wave system in a target which is finite in all its dimensions is too com- 
plicated for analysis at present ; calculations based upon a one-dimensional 
model are, however, very useful, since they represent the actual conditions 
for a limited interval of time depending upon the ratio of width to depth 
of the target, the area over which the incident shock wave may be con- 
sidered to be plane and the distance behind this shock wave over which the 
conditions are uniform. 

The shock wave will be supposed to be of such a magnitude that the 
pressure transmitted to the target block remains below the elastic limit 
ofthe block. When the medium ahead of the block is water, this restricts 
the magnitude of the impinging shock wave to about one-half of the elastic 
limit of the block, since the pressure is approximately doubled by impact 
(see, for example, Cole 1948). The magnitude is therefore restricted to 
about 1500 atmospheres for a shock wave passing through water and 
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falling upon a target made of steel ; the restriction is greater when air is 
the medium through which the shock wave advances, since the multiplica- 
tion on impact is more important. A theoretical treatment of the motion 
when plastic flow is present could be made, but the difficulties would be 
greatly enhanced, and a first attempt on the problem may well be con- 
cerned only with the propagation of elastic stresses in the target. 


§2. WavE Motion IN A PLATE 


Suppose that at time ¢ the plane section S originally at distance x from 
the front surface has been displaced through a distance &. Its displace- 
ment from the original position of the front surface is then a+€. The 
displacement of the plane section originally at distance #-+6a is 
a+é+6x+(d€/ex)dx to the first order in small quantities 62. Thus, 
the new width of the section originally of thickness 6x is 6x(1+0€/dz). 
The engineering strain is therefore 0é/dx. If the pressure in the target 
is p(x, t), then the resultant pressure on the section of the target, opposing 
the motion of the target, is (Op/dx)dz. The velocity of S is u=d€/dt. 
Let ry be the equilibrium density of the target. Then the equation of 
motion for the section of thickness 6z is 


ae ap 
ery aay hee Ate he ten wee mers (1) 


The pressure 7 is related to the strain ¢,—0£/0x by Hooke’s Law. When 
the lateral dimensions do not enter the problem (plane strain), there is no 
strain except in the longitudinal direction, i.e. e,=e,—0. If vis Poisson’s 
ratio and p,=vp, are the principal lateral pressures, 
HKe,=—p-+2vp, ; 
0= —Dy upp), 
giving He,——(1-+v) (l—2v)p/(1—v). Hence 
df _ —(1+v) (I—2v) p_ —(1+»)p 
0x l—y EE 3(l—v)K 
where K is the modulus of compression (bulk modulus), and eqn. (1) 
becomes 
G7 3K(1—v) 07& 
ro ae => el Se aia ° ee OS ec Ay (2) 
Write 3K (1—v)/ry (1+-v)=c?, where c is the general speed of sound for 
dilatational waves ; then (2) becomes 


ae Oa sth x. kG) 


This is the familiar form of the wave equation in one dimension. Since 
the independent variables here are x and t, displacements (and hence 
pressures, etc.) are propagated with the speed of sound relative to the 
equilibrium positions of the particles of the target. It follows that a pressure 


Reflection and Transmission. of Shock Waves : II 191 


pulse crosses a target of equilibrium thickness 1 in time I/c, whatever the 
motion of the target. The eqn. (3) has characteristics 


dx : 
ae +c, ie. x+ct=constant, 


and on these characteristics certain conditions of compatibility are 
required to be satisfied (see Courant and Friedrichs 1948). Thus 


(4) 


The density r in the section S is given by the equation of continuity 


TOL Tor (1+ 9) A 


on x—ct=constant, rycu+p=constant, 
and on x-++-ct=constant, rycuw—p=constant. 


Ar dé 
H Lea SS 
ence paar a ciao) 
where Ar=r—ry. 


§ 3. THe PassaGE OF THE SHOCK WAVE THROUGH THE TARGET 


The front surface of the block in what follows is the surface which, in 
the equilibrium position, is facing the shock wave. The time ¢ is measured 
from the instant when the shock wave reaches the front face. Let 
there be an excess pressure p, in the shock wave. When the shock wave 
strikes the block it is reflected with pressure p, (say), and a shock wave 
of pressure p, is transmitted to the block. The value of the pressure 
Pz in terms of p, is determined by the requirement that the particle 
velocity uw, set up in the block must be equal to the particle velocity in 
the ‘ cuter medium ’ adjacent to the front face. The formulae governing 
the reflection phenomena have been given in Part I (Pack 1957). The 
velocity wu, set up in the target by the pressure p, is equal to p,/rgc, as 
is immediately seen from a characteristics diagram in the (x, ¢) plane by 
following a characteristic of the family x-+-ct=constant which comes from 
any part of the target which remains in equilibrium ; on such character- 
istics the value of the constant in the condition of compatibility is zero, 
since u=0 and the pressure p in the target is neglected in the equilibrium 
state of the target. 

In order to show the motion as it appears to a stationary observer, the 
coordinates used in the figure are x+é, for the displacement of the target 
relative to axes fixed in space, y for distance in the outer medium measured 
(positive in the same direction as x positive) from the original position of 
the boundary of the target. 

The x-coordinates of the target are traversed with velocity ¢ by the 
shock wave front. The 2+ é-coordinates are therefore crossed with 
velocity c-u,; the transmitted shock wave is represented by the line 
OA in the figure. The reflected shock wave is represented by OL. 

The characteristics in the outer medium have the equations dy/dt=u +a, 
a being the local speed of sound (see Courant and Friedrichs 1948). 


O02 
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On each side of OL the values of w and a are constant, although the values 
are different on the two sides. Thus the characteristics are straight 
lines in the outer medium with gradients equal to the local values of wa. 


t Pd 
«©°e) m 


OBDE is THE PATH OF THE FRONT FACE OF 
THE Target, PACE THE ParH OF THE REAR Face 


+] 


For pressures such as are considered in this paper, w<a. Also, when the 
multiplication of pressure on reflection is not too great, the value of a, 
is not very different from the speed W, of the receding shock wave ; 
this follows from the approximation Z,=u,+W,=a, discussed in 
Part I. Thus one family of characteristics between OL and OB (the 
front face of the target) is nearly parallel to OL. However, a, is actually 
greater than W,+w,, representing the physical fact that any disturbance 
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behind a shock wave is propagated faster than the shock wave and 
therefore eventually overtakes it. 

The wave in the target, represented by OA, reaches the rear surface 
and is reflected. The pressure in the medium x>/ will be neglected in 
comparison with the pressure in the target. Thus the boundary con- 
dition on x=/ is p=0. 

On a characteristic «—ct=constant in the region OAB, 


TyCU+ P= CUy+ Po. 
In the region ABC p=0 ; hence, in this region 


Pe 
ToC 


The arrival of the reflected wave AB at the front face of the block is the 
signal to the outer medium indicating the finiteness of the target. The 
rarefaction (py=0) which has passed through the target at the time t, 
— leads to a reflected shock wave in the target and a transmitted rarefaction 
wave in the outer medium. In view of the steady flow behind the receding 
shock wave OL this rarefaction wave is a simple wave (see Courant and 
Friedrichs 1948), represented in the characteristics diagram by a fan 
MBN. A characteristic such as « Syd penetrating this fan becomes 
straight again after crossing BN (the last wavelet of the expansion at B). 
The pressure in the reflected wave BC and the transmitted wave BN is 
P4 (<p,), and the value of p,, together with the velocity of the particles 
following these waves, may be obtained from the characteristic relation- 
ships (see Courant and Friedrichs 1948), thus : 


P ad 
on afyd, ut | a ta ee, ee See OTE) 
p? 


ps being the density in the outer medium corresponding to the velocity 
Ww, and pressure py ; 
on A'B’, r,cu—p=constant=2rceu, . . . . . (66) 


The pressure p, and the velocity uw, are obtained by solving (6a) and 
(6 6) simultaneously. 

From (66), u4=2u.+,/"9¢. Only a small change in the velocity is 
possible in the solid target ; but small changes in velocity in a slowly 
moving gas are associated only with relatively small changes of pressure. 
We can therefore calculate p, approximately by writing pa=p2 (1—«) 
and neglecting higher powers of « than the first. . 

If the outer medium is a gas with p=(kp)’, k being a constant (since the 
entropy is constant behind a uniform shock wave), then 
2a 


al 


9 
) ? L404 
FS 4 a+ 5 Slat 
Toc y—1 


which leads to the result 
2yPo/ToCAe 
1+-ypo/ToCde 
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For a steel target, c—6 x 10° cm/sec and r9=7-8 g/em*. Ifp.=1 kilobar, 
a,=3 X 104 cm/sec and y=1-4, then 


TCtg ’ 


Hence we may well expand « in powers of p,/7)ca, and we have 
2yP > | 
pe=p,( 1-222). 2 eee aay 
The change in p is therefore of the order of 2% at this reflection. If the 
outer medium is water with the approximate equation of state p= B(p”—1), 
in which B=3-047 kilobars and y=7:15, then the formula equivalent to 


(7) is 
ee fi te al, 


TCAs 


For a steel target, if p,=1 kilobar as before, a,=1-7 x 10° cm/sec, and 


2y(Pet+ 5) 
TCAs 


== 0-074. 


The change in pressure is therefore of the order of 7°%. 


$4. Tur CYCLE oF EVENTS 


The phenomena described in §3 repeat themselves over and Over 
again; the proportionate fall of pressure at the front surface varies 
slightly at each step until such time as the first rarefaction (at B) has 
reached the shock wave OL and its reflection returned to the front 
surface to end the succession of steady states which have prevailed at this 
surface between each simple wave. It was shown in Part I that 
Ut W=a, if z,=p,/p, is not too much greater than unity; thus the 
rarefaction wave from B (which is the first to return) is almost parallel 
to OLand since the gradient of the characteristics in the target is in general 
so much steeper than the gradients in the outer medium, there will be 
many reflections of the type discussed above before the analysis is no 
longer valid. A calculation indicates that for shock waves in water, of 
strength up to 1500 atmospheres, there are at least 13 rarefaction waves 
(i.e. the target is crossed by waves in both directions at least 13 times) 
before the shock wave is affected by the rarefaction. The cycle of motion 
of the target will continue unchanged until this wave returns to the target. 
The velocity of the rear face of the target increases at intervals of 2J/c 
by an amount equal to p/ryc, where p is the latest pressure to arrive at the 
rear surface and to be removed by this free boundary. The velocity 
after n reflections is therefore 


U=(PotPstPat -.- +Pnys)/toe- 
When the outer medium is water the ratio of p,,,,/p,, tends to a constant 


limit, since «> 2yB/rgcay (~0-6 for a steel target) where a, is the speed of 
sound corresponding to atmospheric pressure (p=0). The target has a 
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limiting speed and the ‘ acceleration ’ of the target (as measured by the 
non-uniform jump in speed per double-traverse of the target by the 
wave system) is diminishing and tending to zero in the limit. When the 
outer medium is air the same phenomena will occur, but in this case 
e—0 as p.—0. 
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SUMMARY 


_ A combination of lens spectrometer and diffusion cloud chamber was 
used to study the difference in the scattering of 0-3 Mev positrons and 
electrons by argon. 611 electron and 506 positron scatters were observed 
with projected angles $>30°. The ratio of the number of electron 
scatters per unit track length to that of positrons was 1-21 with a standard 
error of 6% ; the theoretical value is 1:18. The measured absolute value 
of the cross section agreed with the theoretical value within the experi- 
mental error of about 10%. 


§ 1. INTRODUCTION 


SEVERAL physical causes make the cross section for the scattering of 
fast electrons and positrons by atoms differ from the simple Rutherford 
formula. The various causes show themselves most clearly at different 
values of the relevant variables. Current interest attaches mainly to the 
finite nucleus effect which is best investigated at very high energies (Hahn 
et al. 1956), and to the radiative effect which should be detectable at 
energies of a few Mev (Elton and Robertson 1952). The effect of spin 
can be detected at energies of the order 1 Mev and should be fully accounted 
for by the well known theory of Mott (1931), provided the scattering angles 
are sufficiently large to avoid the screening effect of the orbital electrons. 
Recent experiments leave little doubt of the correctness of the Mott 
cross section for electron scattering, both in its variation with angle and 
energy, and in its absolute value (see, for example, Pettus et al. 1956, 
Chapman et al. 1955). 

The spin-orbit interaction is responsible for the difference in scattering 
between electrons and positrons (Massey 1942). Consequently, measure- 
ment of the ratio of electron to positron scattering is an alternative way 
of checking the Mott theory. A few such measurements have been 
attempted, of which the most successful is that of Ellis and Henderson 
(1955). ‘Their work was done at energies of 0-7 and 1-4 Mev, using gold 
and aluminium foils. For the three angles 22-8°, 34-5° and 47-5°, it was 
shown that, within an error of 5°/,, the ratio of electron to positron scat- 
termg was correctly given by theory. References to earlier work are 
given by Ellis and Henderson. Since their paper, Bosley and Hughes 
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(1955) have published results consistent with theory for the scattering of 
10 Mev particles in nuclear emulsion. 

The present work gives additional evidence for the electron positron 
scattering difference by another technique and at different values of the 
relevant variables. 


§ 2. THe EXPERIMENT 


2.1. Preliminary Considerations 


For values of Z such that Z/137<1, the series expressions given by 
Mott’s theory can be reduced to a much simpler equation for the cross 
section. In the form given by McKinley and Feshbach (1948), this 
formula is 

1(0)= ae aay 436[1—6? sin?10 in 30 ind 

= Can By sec*}6[1—f? sin?}6+7a8 sin 36(1—sin }6)]. (1) 
The third term in the bracket is taken positive for electrons and negative 
for positrons and is responsible for the difference in scattering. This 
equation should be correct to at least 1° for argon, Z=18. 

It is simpler and more accurate in cloud chamber work to measure, not 
6, but the orthogonal projection of 6 on to a plane containing the incident 
particle direction. Let ¢ be the projection of 6, and P(¢)dd be the number 


bs 
of scatters through ¢ per unit path length; then i P(¢)dé can be 
ie 


computed very easily from eqn. (1) and the tables of O’Ceallaigh (1950), 
if the values of NZ? and f are known from the observations. 

Consider equal lengths of positron and electron track, and scattering 
angles greater than some lower limit ¢,._ Let the total number of electron 
scatters exceed that of positrons by 6N ; and let NV be the total number of 
positron and electron scatters added together. If it is required to detect 
dN, it is desirable that the ratio of dN to the statistical uncertainty in V 
should be as large as possible, i.e. dN/./N should be large. For a given Z 
and track length, this ratio increases as 8 decreases. It is therefore better 
to work with low velocity particles but not, of course, so low that other 
factors such as screening, energy loss, and multiple scattering become 
intolerable. Fortunately, a suitable compromise velocity can be chosen 
within the spectrum of *4Cu which was desirable as a source for other 
reasons as described below. 

Recent papers on screening (Gunnersen 1952, Mohr and Tassie ives) 
suggest that screening can safely be ignored for the values of b> 30 
used in the present experiment. The smallness of the screening effect 
is also indicated by the fact that the accurate measurement of electron 
scattering by Pettus e¢ al. (1956) agreed with unscreened theory down to 
60° at 0-2 Mev in gold. In this connection it should be remembered that 
the ¢’s in the range 30° to 40° will be due to O's anywhere in the range 
30° up to 90°. Incidentally, since the positron-electron difference 
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increases with 6 up to about 80° at the energy of this experiment, it will 
be somewhat accentuated by the use of ¢ instead of 6. 


2.2. Apparatus and Procedure 


The scattering occurred in the argon filling of a diffusion cloud chamber. 
Stereoscopic photography was carried out by conventional methods and 
projected angles were measured on full size reprojectedimages. The errors 
to which such experiments are subject have often been described and need 
not be detailed here ; a recent account has been given by Allen e¢ al. 
(1956). ‘To reduce some of the errors the following methods were adopted. 
No attempt was made to determine the incident energy by track curvature 
measurement ; instead, the particles were injected into the cloud chamber 
from a solenoidal lens spectrometer through a connecting tube containing. 
collimating slits. The spectrometer was calibrated by reference to the 
F-line of the thorium internal conversion spectrum and provided a beam 
of accurately known energy with a resolution of 4% in Hp. Instead of a 
Wilson cloud chamber, a diffusion chamber capable of giving one photo- 
graph every 10 seconds was used. This enabled much more economical 
use to be made of the short lived radioactive source. 

For the source, ®4Cu, of half-life 12:8 hours, was chosen because con- 
siderable importance was attached to scattering the electrons and 
positrons under identical conditions and this was materially assisted by 
a source which emitted both particles with comparable energies and 
intensities. The method of obtaining an output of either positrons or 
electrons from the lens has been described previously (Cusack and Stott 
1955). In the present work, batches of 90 pictures of positron tracks were 
taken alternately with 90 of electron tracks ; in between batches, 10 con- 
trol pictures were taken with the lens current switched off. The source 
was covered with a ray-tight shutter except for the 150 milliseconds before 
each photograph. By this procedure it was ensured that the electron 
and positron scattering occurred under identical circumstances. It is 
not impossible that conditions in a diffusion chamber should change 
slightly during 24 hours continuous work, but such changes would occur 
over much longer periods than the 15 minutes required for one batch of 
pictures and would therefore affect both types of particle equally. Hp 
was chosen at the value for which the beam intensity was the same for 
positrons and electrons so that equal numbers of pictures would give 
equal track lengths. In °*Cu the peaks of the positron and electron spectra 
are at Hp=2000 and 1400 gauss cm respectively, but there are two elec- 
trons to one positron approximately and equal intensities were obtained 
from the spectrometer at Hp=2250 gauss cm. 

At the cloud chamber end of the collimating tube the lens vacuum was 
separated from the atmosphere by a 0-75 mg/cm? terylene foil window 
supported on a brass grid of 80°% transmission. The window in the cloud 
chamber wall was made of a similar foil and the two windows were 
separated by about 2 mm of air. The multiple scattering and energy loss 
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of the particles in the windows were small but would in any case be 
independent of the sign of charge to a high degree of approximation. 

The considerations of the preceding two paragraphs make it practically 
certain that no systematic difference between electron and positron 
results could be introduced before scattering occurred. As a further 
precaution, the scanning and measurement of the photographs was done 
by observers who did not know which were positron and which electron 
pictures. The ratio of positron to electron scatters was therefore subject 
to statistical errors only. Among the errors which tend to cancel out 
when the ratio of positron to electron scattering is considered are : errors 
in NZ? due to uncertainty in the exact value of the supersaturation in the 
sensitive layer ; errors due to scanning inefficiency even if the latter is a 
function of ¢ as, indeed, it probably is ; errors arising from the fact that 
all the incident particles do not move exactly right-angles to the camera 
axis as they should if O’Ceallaigh’s orthogonal projection tables are used. 

A more detailed discussion of the apparatus and of the suitability of a 
glass walled atmospheric pressure diffusion chamber to this type of work 
can be found elsewhere (Bowness 1956). 


2.3. The Observations 


1533 positron and 1563 electron pictures were taken in a single run. 
The current through the lens coil and the pressure and temperature in the 
diffusion chamber were kept constant and measured. From these 
observations the following quantities were obtained. 

(i) The total track length. To find this, certain selection criteria were 
applied to the tracks. The diffusion chamber is continuously sensitive 
to background from cosmic radiation and annihilation rays emanating 
from the source. Most of the latter had been removed by a suitable 
distribution of lead in and around the spectrometer. Of the remaining 
background, a considerable fraction consisted of dense highly contorted 
electron tracks which could easily be distinguished and ignored. All 
other tracks were counted if they were more than 1 cm long and if their 
initial direction made an angle of less than 45° with the incident beam 
direction. These two criteria were applied with the help of a rigid 
Perspex arm which could both slide and rotate in the plane of reprojection 
about a point corresponding to the focus of the spectrometer. To 
indicate the critical length and direction, the end of the arm carried 
fiducial marks which could be superimposed on any track in any part of 
the image. The criteria were arrived at after trials based on the examina- 
tion of the 10% of the pictures which were taken with the lens focusing 
current switched off. These pictures usually contained one or two pieces 
of track and after applying the above criteria to them it was found that 
background accounted for 4% of the accepted track length. Finally, 
tracks showing electron-electron collisions were accepted only as far 

e fork. ; 
- ‘as was not necessary to measure the length of every single track. All 
tracks were measured on a sample of the pictures which was chosen 
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sufficiently large to make the statistical error in the total track length 
small compared with that in the observed number of scatters. This 
condition required that all accepted tracks on 30% of the pictures should 
be measured. The results are given in table 1. 

(ii) The particle velocity. The lens current was measured with a 
potentiometer and it was found from the calibration that Hp= 
2246 gauss cm for the focused particles. Small corrections for the average 
energy loss in the windows and the cloud chamber gas were calculated 
from Landau’s formula. The value of 8 for insertion into eqn. (1) was 
0:783 which corresponds to Hp=2154 gauss cm or energy 0-311 Mev. 

(iii) NZ2. The chamber was at atmospheric pressure throughout and 
the condensant was methyl alcohol. The temperature of the sensitive 
layer was known from auxiliary measurements with thermocouples. 
It is not easy to find directly the vapour pressure but it can be estimated 
from the degree of supersaturation known to be required for track forma- 
tion. However, any slight uncertainty here is not important because the 


Table 1. Analysis of Accepted Track Lengths 


Mianbecet |Namber anal ere g Total 


Particle fra i per frame | track length 
_ frames sample fen (meties) 

Electrons 1563 482 27-66 432-3 

Positrons 1533 454 28-15 431-5 


condensant only contributes 2-4°% of the value of NZ?, and, as pointed 
out above, does not enter the comparative measurements. The final 
value of NZ? was 8-44 x 10?4+1° per cm3. 

(iv) The scattering angles. Scattering angles ¢>30° were measured 
directly on the reprojected images with a protractor. The mean of four 
readings was taken in each case and the error in measuring any one of 
them should not exceed 2° and will often be less. The quantity dN/./N, 
referred to above increases as 4 decreases, nevertheless it was decided not 
to measure angles down to 20° or less because smaller angles are less 
easy to detect with certainty and might also involve consideration of 
screening. The result of the angle measurement is given in table 2 
together with the expected angular distribution of the same total number 
of particles. A little calculation shows that, of the observed p'8>30°, 
approximately 70° were derived from @’s>40°, 46% from @’s>50°, and 
30%, from 6’s>60°. Examination of the control pictures shows that 
45°, of the scatters might have been due to background. 


§ 3. DISCUSSION AND COMPARISON wITH THEORY 


Evidence for the excess of electron over positron scattering can be 
obtained from the ratio of the total numbers of scatters with p> 30° 
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for equal track lengths. Denote the ratio of electron to positron scatters 
by r. Then 


r(obs.)= eae Xx seh == 1-205 
506 432-3 
The standard statistical error in each of the track lengths is 1:6% and 
in the numbers of scatters 4-3°4. The standard error in r(obs.) is there- 
fore 6%. 
From eqn. (1) and O’Ceallaigh’s tables, r(theor.)=1-181, giving 
r(obs.)/r(theor.)=1-02+6%. The observed difference between positron 


and electron scattering is therefore accounted for satisfactorily by the 
Mott theory. 


Table 2. Angular Distribution of Scatters 


Electrons Positrons Both 


Angular 
eae bas |e theors | Obes |\Thear |\COball+Theor. 


30-40 286 | 279-1 222 | 234-4 508 | 513-5 
40-50 150 | 127-2 135 | 104-5 285) 231-7 
50-60 SI 67-9 YS) 55-0 136 | 122-9 
60-70 23 40-3 30 32-4 53 PPE 
70-80 20 25:8 17 20:7 37 46-5 
80-90 23 17-5 a 14-0 30 31:5 
90-100 7 12-4 7 10-2 14 22-5 
100-110 “6 Was aa a ii 7:6 12 16-8 
110-120 7 al 5) 5:9 12 13-0 
120-150 5 14-4 15 12-4 20 26:8 
150-180 4 10-1 6 9-0 10 19-1 
30-180 Gli |" GLL-O DUG OUGO aE Tie iO 


Using the values of NZ? and 8, and the track lengths, the total expected 
number of scatters can be calculated and compared with the number 
found. This provides a check on the absolute value of the cross section 
I(6). For ¢’s>30°, the ratio of observed to expected numbers was 
0:92-+4% for electrons, 0:90-+5% for positrons and 0-91-+-3% for both 
together. The errors quoted are statistical standard errors, but of course 
these ratios are also subject to some systematic errors which are absent 
from r(obs.). The most important error is probably the number of scatters 
missed in the scanning. Although the pictures were examined twice, 
independently, with the use of both cameras, and every effort made to 
count all ¢’s, the possibility cannot be excluded that some scatters were 
missed, particularly in the range 30° to 40°. It should be recalled that 
O’Ceallaigh’s tables are given as functions of p=A,/2a where ’,=the 
minimum track length accepted after scatters and 2a is the depth of the 
illuminated region. Since the attempt was made to count all scatters, 
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the tables were entered at p=0. However, scatters with very short 4, are 
the more likely to be missed and if, for example, A, is put equal to 2 mm, 
the difference between observed and expected numbers practically 
vanishes. The reasonable conclusion is, therefore, that the observations 
show agreement with the theoretical absolute value of the cross section to 
within an experimental error of about 10%. The figure shows the 
observed and theoretical angular distributions, although this type of 


THEORETICAL 
4 RELATION 


° 
30 40 50 60 70 80 90 100 II0 120 130 140 [50 160 170 
PROJECTED ANGLES = DEGREES 


Angular Distribution of Scatters. 


experiment cannot be expected to provide a close check on the distribution 
shape. The points show the measured number of scatters, electrons and 
positrons together, in 10° ranges of projected angles centred on the points 
and the standard errors are indicated. The full line shows the theoretical 
distribution of the same number of scatters. There is rough agreement 
between the points and the curve although the observed ratio of the 
number of 30° to 40° scatters to that of 40° to 50° ones is smaller than 
theory predicts. This is probably due to a combination of statistical 
error with scanning inefficiency in the lower range. 


Electron and Positron Single Scattering in Argon 203 


ACKNOWLEDGMENTS 


The authors wish to thank Professor J. D. Bernal for facilities in his 
department ; the Central Research Fund of the University of London for 
financial support ; the technical staff of Birkbeck College for making 
apparatus ; and Miss J. Locke and Mr. R. Miles for their help with some 
of the routine measurements. 


REFERENCES 


ALLEN, K. R., Lrpsicas, N. M., Masor, D., and Patuirs, K., 1956, Proc. Phys. 
Soc. A, 69, 141. 

Bostey, W., and Hueuss, I., Phil. Mag., 46, 1281. 

Bownsss, C., 1956, Ph.D. Thesis, University of London. 

CHapMaNn, K. R., Matsukawa, E., Ross, P. H., and Stewarpson, E. A., 1955, 
Proc. Phys. Soc. A, 68, 928. 

Cusack, N., and Srort, P., 1955, Phil. Mag., 46, 632. 

Eis, J., and HENDERSON, C., 1955, Proc. Roy. Soc. A, 229, 260. 

Exton, L., and Ropertson, H., 1952, Proc. Phys. Soc. A, 65, 145. 

GUNNERSEN, E. M., 1952, Austral. J. Sci. Res. A, 5, 258. 

Haun, B., RAVENHALL, D. G., and HorstapTsER, R., 1956, Phys. Rev., 101, 1131. 

Masset, H.S. W., 1942, Proc. Roy. Soc. A, 181, 14. 

McKintey, W. A., and Fresupacn, H., 1948, Phys. Rev., 74, 1759. 

Monr, C. B. O., and Tassig, L. J., 1954, Proc. Phys. Soc. A, 67, 711. 

Mort, N. F., 1931, Proc. Roy. Soc. A, 235, 429. 

O’CraLLaiIcH, C., Proc. Roy. Irish Acad., 53A, 58, 133. 

Pertus, W., Buosszr, H. G., and Hererorp, F. L., 1956, Phys. Rev., 101, 17. 


The Elastic Scattering of Protons by Protons at 925 MeV { 


By P. J. Duke, W. O. Lock and P. V. Marcu 


University of Birmingham 


W. M. Greson and J. G. McEwen 
Queen’s University, Belfast 


and 


I. 8S. Hugues and H. MurrHHApD 
University of Glasgow 


[Received July 31, 1956] 


SUMMARY 

[Ilford G5 emulsions have been exposed to a 925 Mev beam from the 
Birmingham proton synchroton. Sixty-five examples of the elastic 
scattering of protons by free protons in the emulsion have been found in 
1539 metres of track. After removal of events having azimuthal angles 
w>45° where % is the angle between the plane of the two secondary 
particles and the plane of the emulsion, the remaining 44 events yielded 
a total elastic cross section of 17-+3 mb. The angular distribution of all 
the events shows a strong forward peak in the centre of mass system. The 
shape of this distribution has been compared with the predictions of a 
simple optical model and a value of (0-9-L0-1) x 10-15 em has been found 
for the interaction radius of the proton. The results are compared with 
other experimental work. 


§ 1. INTRODUCTION 
In this paper we present a detailed account of a study of the elastic 
scattering of protons by free protons in Ilford G5 nuclear research emulsions. 
A preliminary report of this experiment has already been published (Duke 
etal. 1955b). Subsequently the mean energy of the proton beam has been 
more precisely determined, and the track length followed in emulsions of 
known. shrinkage factor has been doubled. 

It is intended that the present discussion should also serve as an intro- 
duction to the two following communications. These concern investi- 
gations of the ‘ quasi-elastic ’ scattering of protons with protons bound 
inside a nucleus (McEwen, Gibson and Duke 1957) and of the inelastic 
scattering of protons with free protons in the emulsion (Hughes, March; 
Muirhead and Lock 1957). 


+ Communicated by P. B. Moon, F.R.S. 
{ Now at Department of Natural Philosophy, University of Glasgow. 
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The results of experiments using protons with energies up to 450 Mev 
have been tabulated by Berretta et al. (1954). The distinctive features 
of these results are the constant value of the total elastic cross section 
above 100 Mev and the isotropy of the differential cross section above 20° 
in the centre of mass system (C-system). Indications of a rise in the cross 
section at angles of less than 30° in the C-system resulting from a purely 
nuclear interaction (as opposed to that resulting from Coulomb repulsion) 
were first observed by Hartzler and Siegel (1954) at 428 mev. This devi- 
ation from isotropy was confirmed at 437 Mev (Sutton et al. 1955). 
Following the construction of the large proton accelerators in the U.S.A., 
the U.S.S.R., and Great Britain, the maximum proton energy at which 
elastic scattering experiments have been carried out has been increased to 
about 6 Bev. The published results indicate that the anisotropy of the 
scattermg becomes more pronounced as the energy of the protons is 
increased. There is also an indication that the value of the total elastic 
cross section begins to decrease slowly with energy above 400 Mev. 
These results will be reviewed in more detail later and compared with 
those of the present experiment. 


§ 2, EXPERIMENTAL PROCEDURE 


The experiment was performed using protons scattered out from the 
circulating beam of the Birmingham synchrotron (Moon et al. 1955). The 
target material was the free hydrogen normally present in Ilford G5 


Fig. 1 


Plan of part of the synchrotron vacuum box and pumping manifolds showing 
the inflector tip and the trajectory of a scattered proton beam. 
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emulsions. The protons emerged through a circular glass window which 
was clamped to the wall of the vacuum box by a steel annulus, as shown 
in fig. 1. The trajectories of these scattered protons in the stray magnetic 
field were computed graphically and are shown in the same figure. The 
exposure position at P along the trajectory was chosen to give a maximum 
value for the ratio of the main beam intensity to that of the background 
of lower energy particles. 

In order to verify experimentally that the beam which had passed 
directly through the window was reaching the plates, standard x-ray 
films were exposed perpendicular to the beam at several points along the 
trajectory. An exposure of about 100 pulses of protons of full energy 
(about 105 protons/cm?/pulse) was sufficient to give a sharp image of the 
proton beam. ‘Two of these x-ray films are shown in figs. 2 (a) and (6). 
That in 2 (a) was exposed at about 15 cm from the window and that in 2 (5) 
at the exposure position. The different areas covered by the beam in (a) 
and (b) were due to the horizontal defocusing and the vertical focusing 
action of the magnetic field. The latter effect, which was also the cause 
of the elliptical shadow of the steel annulus in (a), was not considered when 
the trajectory was computed. Because of this focusing effect and the 
finite extent of the emulsions at P, it was clear from the x-ray films that 
the plates were accepting not only protons of almost the full energy which 
had emerged through the window but also protons of slightly lower energy 
which had passed through the porcelain wall of the step and through the 
steel annulus. The actual extent of these three regions was determined 
from the position of the shadows cast on (a) and (6) by a block of absorbing 
material placed near the window, at A in fig. 1. The fairly sharp line of 
the shadow on film (b) corresponded to the inner limit of the exposure 
position and hence the shadow cast by the same block on (a) indicated the 
region of the step from which the emerging protons were being accepted 
by the emulsions. From the area and thickness of the steel annulus, the 
porcelain wall and the glass window, the mean energy of the beam was 
found to be 925+30Mmev. This value should be compared with the 
previous estimate of 950+3) mev (Duke e¢ al. 1955b). The value of 
925+.30 mev found from the geometry of the experiment is also consistent 
with the value of 870+35 Mev found from measurements on 8 elastic 
scatters in which one of the secondary protons came to rest in the emulsion. 

A few of the recent exposures were made with a beam which was 
scattered out by a carbon target through a similar window further round 
the synchrotron. The relative positions of the target, the window, and 
the emulsions were the same as those described above and so the beam 
energy was assumed to be unchanged. 

The plates were scanned by following the tracks of protons using Cooke 
M4000 microscopes fitted with x45 oil immersion objectives. The 
intensity of the exposures (about 2x 10* particles/em?) corresponded to 
about 8 tracks per field of view (100 square) so that the probability of 
an observer changing from one track to another while following it was 
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negligible. Only those tracks which were within 24° of the mean beam 
direction and of minimum ionization were selected and followed. All 
interactions of protons in the emulsion were recorded ; 2p and 4p stars 
and deviations of the primary proton were subjected to a more detailed 
examination. In this notation a star of type Np has been produced by a 
proton and has N secondary prongs of length greater than 5p. All 
secondary tracks with length less than 5 were classified as recoil tracks. 
A deviation of the primary proton unaccompanied by any observable 
change in grain density was attributed to elastic scattering in the field of 
a large nucleus (Lock e¢ al. 1955). The 4p stars were examined to detect 
any which could have been due to the reaction p+p— p+p +at-+a-. 


§ 3. ANALYSIS OF 2p STaRsS 


The application of the principles of conservation of momentum and 
energy to the elastic collision of a proton with a free proton initially at 
rest yields the following results :— 


(i) The track of the incident proton and the two secondary tracks are 
coplanar, that is, the angle, ¢, between the incident proton and the plane 
of the two secondary tracks is zero. 


(ii) The true angles of scatter, «,, and «,, are related by the expression 
tan a, tan «,=2/H Laie ee REC) 


where £# is the total energy in the laboratory system expressed in units of 
the rest mass of a proton. For protons with an energy of 925-+-30 Mev, 
2/E is equal to 0-670-+-0-008. Thus the true opening angle of the inter- 
action, «,+,, must lie between 783° and 90°. 

These properties were used to separate the true elastic events from 
inelastic events and interactions with complex nuclei. The angles and 
approximate grain densities of the secondary tracks of all 2p stars were 
carefully measured, apart from those stars having a secondary caused by 
an «-particle or a heavier nucleus, or a definite recoil stub. The angle of 
scatter, projected on the horizontal plane, was measured to + 7° using 
the hair line in the eyepiece graticule. The dip angle of each track was 
measured using a X95 oil immersion objective to + 7°. Each time an 
event was measured, the shrinkage factor of the emulsion was determined 
to +1% using a method described elsewhere (Duke et al. 1955 a). The 
error of +3° specified for the dip angle includes that introduced in the 
shrinkage factor determination. The angles «, and «, were obtained from 
the measured angles by : 


cos a—=cos @ cos § cos §)+sin dsinds . . - - + (2) 
and ¢ was given by : 


G08 dy COS 31 COS dy [tan 8, sin 9,—tan 6, sin 6,+tan do sin (0,—6,)] 


Sk aes (a, + ag) (3) 


P2 
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where 8, 5, and 6, are the dip angles of the primary and the two secondary 
tracks and 6, and @, are the projected angles of scatter. A measure of the 
departure from the perfect angular correlation required for an elastic 
collision was obtained by inserting the measured value of the larger « into 
eqn. (1) and hence calculating the value of the smaller «. The difference 
between this value and that obtained by direct measurement was denoted 
by da. For a true elastic event 4a should be zero. The estimated 
average statistical errors for both the measured and the derived quantities 
are given in table 1. The values stated include an allowance for a small 
emulsion distortion. In particular cases systematic errors of about 1° 


Table 1 
Projected angle 6 +15’ 
Dip angle ) +45’ 
True angle of scatter oy +50’ 
Departure from 
angular correlation Au SEN 5’ 
Coplanarity angle d +1° 
Fig. 3 
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Departure from angular correlation Aa 


Graph of the coplanarity angle ¢, plotted against the departure from 
angular correlation, da. 


may be introduced by large distortion effects, a small angle scatter of a 
secondary near the event or random irregularities in the fine focus move- 
ment of the microscope. The statistical error itself would be increased 
if the track being measured left the emulsion after less than about 200 u 
had been traversed. 

All those events which had «++, between 70° and 95° and gave values 
of 6<5° and |4a|< 5° were plotted on the graph of ¢ against da shown 
in fig. 3. The large limits were taken to ensure that no true elastic inter- 
action could be missed. There is a definite clustering of points near 6=0 
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and d4a=0. Semicircular arcs of radii 0-5°, 1-:0°, 1-5° ete. centred on 
da=0 and ¢=0 were drawn on fig. 3. The number of events per unit 
area included between each pair of arcs plotted against the radius of the 
are is shown in fig. 4. The distribution has a roughly Gaussian form 
wnich approaches zero at a radius of 1-5° and continues thereafter as a 
tong tail ’, indicating a random background distribution. The radius of 
1-5 was taken as the cut-off value, and all the events outside the cor- 
: responding semicircle were assumed to be inelastic events or interactions 
with complex nuclei. This value of 1-5° for the cut-off is consistent with 
the estimated errors in ¢ and 4« shown in table 1. 


Fig. 4 


60 
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No. of events per unit area 
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Arc radius 


Graph of the number of events per unit area included between successive arcs, 
drawn on fig. 3, plotted against the arc radius. 


Of the 122 events shown in fig. 3, 57 fall outside and 65 inside the semi- 
circle. These 65 events were taken to be examples of the elastic 
scattering of the incident protons by free protons in the emulsion. 


$4. SCANNING EFFICIENCY 


A histogram of the total number of elastic events in four equal ranges 
of azimuthal angle, ;, defined as the angle between the plane of the two 
secondaries and the emulsion plane, is shown in fig. 5(#). There is a 
definite decrease in the number of events with large values of ¢. Figures 
5(b) and (c) are similar histograms for those events having the smaller 
angle of scatter «, >12° and <12° respectively. The loss is more marked 
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in fig. 5(c) although it is still appreciable in (6). It appears that a con- 
siderable number of events having values of s>45° have been lost in the 
scanning or in the analysis. The lost events are predominantly those 
having «,<12°. 

There are several reasons for the loss of those events having large values 
of sand small values of «,. It is inherently more difficult to detect events 
which take place in a plane almost perpendicular to that of the emulsion. 
In these cases the main component of the small deviation of the proton 
primary is in the dip angle. Such deviations are difficult to observe when 
scanning along the track. Furthermore it is these events which are most 
likely to be affected by errors in shrinkage factor measurement, by the 
shortness of the recoil track, by distortion, or by small angle scattering of 
one of the tracks near the event. 


Fig. 5 
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Histograms of the number of events in equal intervals of ys, the azimuthal 
angle ; (a) all events, (b) events having the smaller angle of scatter a, 
> 12°, (c) events having «,<12°. ‘ 


In order to calculate the total cross section we have assumed that there 
is no loss of events with 4<45°(44 cases) and that the number of events 
with any y is twice the number with p<45°. 
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§ 5. Resuuts 


The 44 elastic events having <45° were found in 1539 metres of track 
length. Two events were subtracted from this number to allow for those 
background events which were included within the semicircle leaving 
42-+-7 elastic collisions with y4<45° and therefore 84+-14 events with any 
% value. The mean free path for the elastic scattering of protons in 
Ilford G5 emulsions is therefore 18-3 +3% metres. The free hydrogen 
content of the emulsions used in the experiment was obtained from data 
supplied by Lford Limited (Waller, private communication), together 
with a knowledge of the humidity at which the plates were stored im- 
mediately prior to their exposure. The mean value of the hydrogen 
content was (3:17-0:06) x 10”? atoms/cm®. From this figure a value of 
17+3 mb was obtained for the total elastic cross section. If the loss of 
events with large values of % was neglected, the 65 events gave a total 
elastic cross section of 13-43 mb. 


Fig. 6 
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Differential cross section plotted as a function of the angle of scatter in the 
centre of mass system. 


The angles «* of the secondary tracks in the C-system with respect to 
the direction of the primary may be calculated using the relationship : 


. eait 
tan —— = [£ UR ae eee ee ee ee ee) 
2 N 2 
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For a 925 Mev proton eqn. (4) reduces to tan $4*=1-22 tana. Values of 
«* have been obtained for all 65 elastic eens The differential cross 
section has been calculated for six equal intervals of solid angle after 
making a correction for the loss of those events having large values of q. 
This correction was made by including in the first angular interval only 
those events which had <45°, giving 23 events or 46 for all values of ¢. 
The total number of events for all 6 intervals was increased to 84 (the 
number used in obtaining the total elastic cross section) by distributing 
the 3 extra events required amongst the other 5 angular intervals in the 
ratio of the number of events actually observed. The result of this 
calculation is shown in table 2 and the values obtained for the differential 
cross section are plotted graphically in fig. 6. 


Table 2 
Angular Range Observed Corrected da 
C-system degrees} number of events | number of events dQ 
in mb/ster. 

0-334 30 46 8:-941-3 
333-48} 22 23-9 4-6 +0-9 
481-60 5 54 1-0+0-4 
60-705 4 4-3 0-8 +0-4 
704-804 2 2-2 0-4+0:3 
803-90 2 2-2 0-4+0:3 


§ 6. Discussion 

The experimental values of the total cross section for elastic scattering 
obtained at energies between 437 and 5700 Mev are givenin table 3. There 
appears to be a tendency for the total elastic cross section to decrease slowly 
with increasing energy. It should be noted that in all the experiments 
carried out using counter techniques, the value obtained for the elastic 
cross section depends on an extrapolation of the observed angular distri- 
bution. Measurements were made only at angles greater than about 
30° in the C-system in all experiments except that of Meshcheryakov 
et al. (1956) who made measurements down to 5°. 

The values obtained in the present work for the differential cross 
section are in general agreement with those obtained by Smith e¢ al. 
(1955) at 1000 Mev. The increase in the value of the scattering cross 
section at forward angles is part of a general tendency which begins at 
about 430 Mev (Sutton e/ al. 1955) and continues up to the present 
experimental limit of about 6 Bev (Cork and Wenzel 1955). Reference 
should be made to the report of Hess (1956) for a complete summary 
of the differential cross section measurements at different energies. 

It has been pointed out by Serber and Rarita (1955) that the results of 
Smith et al. (1955) for the differential cross section can be fitted by the 


expression do J (AR sin «*) 
d§d sin «* | 


(5) 
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where k is the wave number of the incident proton in the C-system, R is 
an iteraction radius of the proton, and ./, is a first-order Bessel function. 
We have fitted eqn. (5) to our results and have obtained a value of 
(0-9+0-1) x 10-18 em for R. Thisisin good agreement with that obtained 
by Batson and Riddiford (1956), Fowler et al. (1956 a, b) and Chen e¢ al. 


Table 3 
Se eee ll ol EI nn 


Energy _ | Total elastic cross 


(Mev) section in mb Technique Reference 
437 23°8+1-2 C Sutton et al. (1955) 
440 23-5+1-2+ C Smith et al. (1955) 
460 24-0 +0:6 C Meshcheryakov et al. (1956) 
560. 26-:0+0:8 GC Meshcheryakov et al. (1956) 
590 25-2 +2-07 C Smith et al. (1955) 
650 +50 26:3+1-8 Cl. Ch. Batson and Riddiford (1956) 
660 23 +2 C Bogachev and Vzoroy (1954) 
660 24-7+1-0 C Meshcheryakov e¢ al. (1956) 
800 21-5+2-07 C Smith et al. (1955) 
810+ 100 24 13 Cl. Ch. Morris et al. (1956) 
: Em Present work. 
Cc Smith et al. (1956) 


Cl. Ch. Fowler et al. (1956 a, b) 
Cl. Ch. Block et al. (1956) 

Cl. Ch. Wright et al. (1956) 
Km Giles (1955) 


+ Normalized to the value of do/dQ at 90° obtained by Sutton et al. (1955) at 
437 Mev ; the values given are those quoted by Hess (1956). 
Key: C, counters ; Cl. Ch., Cloud chambers ; Em, Emulsions. 


(1956) from a study of proton scattering at energies between 600 and 
2700 Mev. Interaction radii of the same order of magnitude have been 
obtained from pion—proton scattering experiments at 1-4 Bev (Hisberg 
et al. 1955, Walker and Crussard 1955). From electron scattering experi- 
ments McAllister and Hofstadter (1956) have obtained a value of 
(0-74-40-24) x 10-18 cm for the r.m.s. radius of both the charge and the 
magnetic moment distributions of the proton. 
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SUMMARY 


A study has been made of the production of 7-mesons in the collisions 
of protons having an energy of 925 mev, with the hydrogen present in 
nuclear emulsions. The results yield the following cross sections : 


p+p—+>a7t+tn+p 27-3 mb, 
at+D 0-6-L0-4 mb, 
7+p+p 6+2 mb, 


mtta—t+p+tp <0-2 mb. 


Angular and momentum distributions have been obtained for particles 
emitted in the first process. The results may be interpreted in terms of 
the production of mesons by the decay of an excited nucleon. 


§1. InTRODUCTION 


RECENT measurements of the magnitude of the total cross sections for 
proton—proton collisions have shown a rise from the nearly constant value 
of about 23 mb over the region from 150 to 400 Mev (Berretta et al. 
1954) to about 48 mb at 900 Mev (Shapiro et al. 1954). Above 900 Mev 
the cross section decreases slowly to at least 2-7 Bev and reaches a value 
of 41 mb at that energy (Shapiro 1956). 

In contrast to these measurements, those on the cross section for the 
elastic scattering of protons have shown a slow decrease from a value of 
25 mb at 400 mev to about 15 mb at 2-75 Bev (Block et al. 1956; Smith 
et al. 1955). 

The difference between the cross sections for the total and the elastic 
scattering represents the cross section for meson production. ‘The present 
paper describes measurements on the characteristics of meson production 
when the incident proton possesses an energy of 925+30mev. The 
measurements were made in nuclear emulsions, which were exposed to 
the external beam from the Birmingham proton synchrotron. 


+ Communicated by the Authors. 
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In table | the possible processes are given for the production of 7-mesons 
which may occur at this energy. In the present experiment the energy 
in the centre of momentum C-system is 429mev; that required for 
the production of one, two or three 7-mesons is about 140, 280 and 420 Mev 
respectively. No attempt was made to find examples of triple production. 

The method of identification of events representing 7-meson production 
utilizes the principles of the conservation of energy and momentum. 
Hence processes involving the emission of two uncharged particles could 
not be distinguished. 


Table 1 
p+p > w*-+n+p (a) 
a4) (5) 
n+p+p (c) 
is ws ol (d) 
CoE widen tse (e) 
a a oR aA) 
na Pp (9) 
3-meson production (A) 


§ 2, EXPERIMENTAL PROCEDURE 
2.1. Preliminary Examination 


The methods used for the exposure and examination of the nuclear 
emulsions are described in the previous paper (Duke ef al.). All events 
having two or four secondary charged particles were examined. Those 
accompanied by the emission of an electron or in which one of the tracks 
was produced by an «-particle or by a nuclear recoil were rejected. 
All events identified as elastic scatterings between free protons were also 
excluded. 

In proton—proton collisions at 925 Mev the maximum angles of emission 
of the nucleons with respect to the primary direction are 57° and 35° 
respectively, in the laboratory system of coordinates, for the emission of 
one and. two 7-mesons. 

If an event satisfied these conditions the angles of emission and the grain 
density of each secondary track were measured. Further events could 
then be rejected since the energy (and hence the grain density) of the 
. protons emitted in processes (a), (c) and (d) had to fall within prescribed 
limits for any permitted angle of emission. The simultaneous production 
of a 7+-meson and a deuteron could be recognized at this stage since 
unique relations exist between the angles of emission and the grain 
densities of the meson and the deuteron; the tracks of the primary 
proton and the two secondary particles must also be coplanar. 
Considerations of the balance of energy and momentum permit the 
identification of process (d) since all four secondary particles are observed. 


a-Mesons in Proton—Proton Collisions at 925 MeV 217 


2.2. The Method of Identification for Hvents 
Involving One Uncharged Particle 


7 The details of the procedure used to balance energy and momentum 
are given in an Appendix. The probability that an event can conserve 
energy and momentum is expressed in terms of a ratio R, defined by 


Guo | 
PACs, 


where Gy is the measured grain density for one secondary track, 4G, 
is the maximum error of this quantity, and @ is the grain density which 
this track must have if energy and momentum are conserved. A gaussian 
distribution is expected for the values of this ratio. The calculated 
distribution of R shows that the majority of genuine examples of 7-meson 
production should have values of R<1-5. Whenever possible, measure- 
ments of the multiple scattering of the secondary tracks were made in 


Fig. 1 


G.D. 


.e) 15 30 45 60 75 90° 
e 


Distribution in grain density and angle with the primary, of secondary tracks 
assumed to be mesons for which a positive indentification, by means of 
multiple scattering measurements, was not possible. Each cross 
represents one track. Proton tracks arising from elastic scattering of 
the primary particle by a free nucleon must lie on, or to the left of, the 
curve drawn. 


order to decrease the limits of error of the energy determinations and 
to check the assignments of identity made on the basis of energy and 
momentum balance. As a result of these measurements new values of R 
were obtained in a few cases. It was found that of 136 of the events with 
R<15, the track assigned to a 7-particle was positively identified as 
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am7-meson in 36 events and a proton in 44 while for the remaining 56 events 
the measurements were inconclusive. 

If the unidentified events are accepted or rejected in the ratio of the 
numbers of identified mesons to protons the cross section for the process 
ptp—-a7t+n+p is 18mb. If all the uncertain events are accepted 
this cross section becomes 27 mb. 

We believe that the latter value is probably the better one for the 
following reasons : 


(i) An examination of the unidentified events shows that in the majority 
of cases the track assigned to a z-meson has a large angle of dip and a low 
erain density. If these tracks were due to protons they would have an 


Fig. 2 


Distribution of the ratio R for events not rejected by multiple scattering 
measurements. In this and the following figures the full line histograms 
represent the distribution for all the events, the dotted line histograms 
those for events in which the 7-meson was positively identified. The 
full line curve is the calculated distribution for all events. 


energy in excess of that allowed for an elastically scattered proton at 
the same angle of emission (see fig. 1). Since the background events 
can only arise from complex nuclei in the emulsion an improbable sequence 
of collisions inside a nucleus would be required to produce such protons. 


(ii) Most of the measurements of multiple scattering were made on 
tracks of low grain density. Ifa proton produces a track of low grain 
density the conservation of momentum requires that it should be emitted 
at a small angle to the direction of the primary proton. Therefore, 
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such protons tend to produce the long tracks in the emulsion which are 
necessary for multiple scattering measurements. Since mesons which 
produce tracks of low grain density may be emitted at any angle, multiple 
scattering measurements are biased in favour of protons. 

It was found that the majority of the 44 events in which the track 
assigned to the meson was in fact a proton, were events in which the range 
of error in the determination of the energy of one of the secondary 
particles was large. 

In view of these considerations we have accepted the 56 events in 
which the multiple scattering measurements were inconclusive, as 
genuine examples of the process p+p—>a++n-+p. 

The distribution of R for all events in which R is less than 5, except 
those rejected by multiple scattering measurements, is shown in fig. 2. 
The calculated distribution has been normalized to the total number of 
events with # less than 1-5. 

The method of identification we have used for events in which one 
meson is produced is more selective for the production of 7°- than for 
a@*-mesons, ‘This is so because energy conservation restricts the possible 
energy values of the neutral particle, and if this particle is a meson the 
range of possible momenta is smaller than that for a neutron. Thus 
the region over which a momentum balance is possible is restricted for the 
production of 7°-mesons. For the production of 7°-mesons, however, 
there exist no criteria equivalent to the identification of 7+-mesons by 
measurements of multiple scattering and by the angle-grain density 
plot (fig. 1). 

§ 3. Resunts 


A total of 453 and 298 events having two and four secondary particles 
respectively were found. In table 2 we display the numbers of accepted 
events remaining at each stage of elimination. These events 


Table 2 
Number of events 2 Prongs 4 Prongs 
Found 453 298 
After preliminary examination 238 0 
With R<1-5 160 — 
After multiple scattering LIEGE — 


+ This includes all events involving single meson production. 


were found in a total length of 1076 metres of track of the 925 Mev 
protons; the hydrogen content at the time of exposure was 
(3-17-+0-06) x 10” atoms/cm’; the cross sections for processes observed 
could thus be found and are listed in table 3. 
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The angular and momentum distributions in the C-system for the 
7-mesons, neutrons and protons are given in figs. 3, 4 and 5. The 
experimental angular distributions were found to be symmetrical about 
90° as is to be expected for a collision between two identical particles. 
To improve the statistical accuracy of the results the distributions have 
been folded about this angle. For all the distributions, we have plotted 


Table 3 
Process No. of events Cross-section 
pt+p- a*+n+p 92 Yi BSA Oy 
a+D Yy 0-6 +0-4 mb. 
7°+p+p 22, 6--2 mb 
aT 4-are 4 Dae 0 <0-2 mb 
Fig. 3 
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Angular distributions in the C-system of the secondary particles from the 
reaction p+p > at+p-+n. 


values of the quantities such that energy and momentum are conserved 
exactly (‘ consistency ’ values—see Appendix), rather than the measured 
values. The small differences between the measured and ‘ consistency ’ 
values are due to the experimental errors but the distributions of both 
sets of values are very similar for all the quantities. 

A satisfactory feature of the results is that the distribution of the 
iia angles shows no evidence of loss of events when this angle is 
arge. 
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Fig. 4 
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Momentum distributions in the C-system of 7+-mesons from the reaction 
p+p+zt+p+n. The full line curve represents the form of the 
distribution calculated from purely statistical considerations. 


Fig. 5 
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distributions in the C-system of the proton and 
Be ine” fecictibn p+p >at+p+n. The full line curve represents: the 
form of the distribution calculated from purely statistical considerations, 
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§ 4. DIscussION 
4.1. Cross Sections 

The measured cross section for inelastic scattering is 34+-4mb. That 
for elastic scattering, obtained from the same exposure and examination 
of the emulsions, is 17-43 mb (Duke et al. 1956). These results compare 
favourably with those obtained using counter techniques (Shapiro 1956 
Smith et al. 1955), but less satisfactorily with those found by using cloud 
chambers (Block et al. 1956). They are displayed in table 4; in the 
region of energies surveyed the total cross section for the scattering of 
free protons remains nearly constant at a value of 48 mb (Shapiro 1956). 


Table 4 
Cross section in mb 


Technique 
Elastic Inelastic 


Counter 26=o3 


Cloud chamber 2A-ee 

Emulsions 34+4 
Counter 2943-5 

Cloud chamber 2 Deo 


The values quoted by Block et al. were obtained by a normalization 
of the observed ratio of inelastic to elastic scatterings to the measurements 
of total cross sections made by Shapiro et al. The cross sections for 
inelastic scattering from counter experiments represent the difference 
between the measured values for total and elastic scattering. 

The results displayed in table 4 suggest that the events accepted by 
us as examples of the production of mesons from free hydrogen are not 
seriously contaminated by spurious examples arising from the more 
complex atomic constituents of the emulsion. 

Unfortunately, firm theoretical predications cannot at present be made 
concerning either the absolute or relative magnitudes of these cross 
sections. The observed angular distributions for the elastic scattering 
between free protons suggest that the scattering is of a diffraction 
character for incident protons with energies between 650 Mev and 3 Bev 
(Smith et al. 1955; Cester et al. 1956). Over this range of energies the 
observed distributions may be successfully fitted by an expression of the 
form [J ,(kR sin @)/sin 0]? where J, is a first-order Bessel function, & is 
the wave number of the incident proton in the C-system and R is an 
interaction radius for the proton (Serber and Rarita 1955). 

Despite the successful representation of the experimental results by 
this simple expression, an examination of table 4 strongly suggests that 
the cross sections cannot be simply described in terms of optical ‘ black 
body’ scattering, since, with the exception of the cloud chamber result 
at 800 Mey, the cross section for elastic scattering is invariably less than 
that for inelastic scattering. | 
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It is of interest to attempt to relate the concept of optical scattering 
to the statistical theory of Fermi (1950, 1953). This theory assumes 
the formation of an excited compound system which may decay into 
its parent particles or into a mixture of mesons and nucleons. The former 
mode of decay will represent compound elastic scattering which would 
interfere coherently with the diffraction scattering. If a compound 
system is formed and decays in the manner proposed by Fermi then the 
observed shape of the differential cross section for elastic Scattering 
suggests either that the interference takes place in such a manner that 
the characteristic pattern for diffraction scattering is preserved, or that 
the compound elastic scattering is negligible. 

Since the probability for the occurrence of any process may be 
represented by an expression of the form | H |?W where H is the interaction 
term for the process and W is a statistical weighting factor covering the 
density of final states, either or both of the above proposals could be 
satisfied by the form and magnitude of | A |?. 

We will compare the relative cross sections for the decay of the 
compound system by assuming the formation of the following three 
states, (a) two nucleons, (b) two nucleons and a 7-meson, (c) a nucleon 
and an excited nucleon. 

State (c) represents the assumption that the production of the meson 
proceeds through the formation of a nucleon in an excited state with 
spatial and isotopic spin quantum numbers of 3/2 (Peaslee 1954). The 
energy available in the C-system for all the experiments listed in table 4 
is well above the threshold for the production of this state. 

We will consider the extreme assumptions that either the states (a) and 
(6) completely dominate (c), or that states (a) and (c) dominate (b). The 
first assumption represents the original statistical theory of Fermi. 
For this assumption 


6 9-43 (Block 1956). 
Oo 


a 


For the second assumption we may write 


o_ [H.W Ue reece 
ie, Pe, Wi Ratlaeme 

where p=density of momentum states available, f= factor to allow for 
conservation of isotopic spin and, S=weighting factor for the spins of 
the systems and where the factor 2 occurs since indistinguishable particles 
appear in the elastic scattering. Inserting numerical values and assuming 
that the mass of the excited nucleon is 1-27 proton masses (see end of this 
section), the weights may be written as 

Le =0:66X 1x 2x 2=2°6. 

Wa 

If we take | H, ? =|H,[? then o,/o,=2°6. Thus it may be seen 

that inelastic scattering is favoured if it is assumed that the production 
of the meson proceeds through the 3/2, 3/2 state. This is in accord with 


Q2 


and 
Oo 
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the suggestion made above that compound elastic scattering should 
be small, but it must be borne in mind that we have assumed | H, ?=|H, |? 
which may not necessarily be true. 


4.2. Relative Yield of 7+- and 7°-Mesons 


Our measured ratio for the yield of 7+- to 7°-mesons is 4:2+1:0. In 
table 5 we list published data for this ratio obtained in experiments 
over the energy range from 650 to 925 mev. The data at higher energies 
suffer from uncertainties through the production of two a-mesons. 


Table 5 

Mnerey. Technique Ratio Reference 
in Mev 

650 Cloud chamber 3°8+1-9 (a) 

—0-9 

660 Counters 3-1+1-0 (b) 

800 Cloud chamber 1748 (c) 

925 Emulsions 4-2+1-0 Present experiment 


(a) Batson and Riddiford (1956). 
(b) Dzhieliopov e¢ al. (1955), Tyapkin et al. (1955). 
(c) Block et al. (1956). 

It can be seen from table 5 that the results show anomalously large 
fluctuations. It should be borne in mind that the measurement from 
the present experiment is almost certainly contaminated by events from 
interactions in complex nuclei, which could conceivably affect this 
ratio. The ratio of 17-8 quoted by Block et al. must be treated with 
caution since positive assignments of events as the production of 7+- or 
7°-mesons could not be made in about 25°, of the events found. 

In view of these reservations no serious comparison can be made 
between the statistical model (scheme (b) above), which predicts ratios of 
1, 3 or co, and the excited nucleon model which predicts a ratio of 5. 


4.3. Angular and Momentum Distributions for p+p—>a7++n-+p 

The experimental angular distributions of the secondary particles 
in the centre of momentum system are approximately isotropic. The 
distribution of the angles between the 7-meson and the proton, shown in 
fig. 6, shows a stronger correlation than would be expected from the 
statistical theory.t There isno marked correlation in the angle between 
the 7-meson and the neutron. The correlations are readily explained 
in terms of an excited nucleon since nine-tenths of the 7-mesons would 
be produced by the decay of the excited nucleon into a 7+-meson and 
a proton and only one-tenth by the decay into a 7+-meson and a neutron. 


{ Unless it is specificially stated in the following discussion we shall assume 
that the statistical theory and momentum space considerations predict the 
same results, . 
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The curves superimposed on the experimental momentum distributions 
(figs. 4 and 5) are those expected from considerations of the momentum 
space available to each particle on the statistical theory (Block 1956). 
It can be seen that there is a tendency for the mesons to be emitted with 
low values of momentum. 

The statistical theory predicts that the momentum distributions of 
the neutrons and of the protons should be identical. The experimental 
distributions appear to differ, that for the neutrons shows a broad peak 
about 400 mev while that for the protons is wider and has a tendency 
to peak at the maximum possible momentum. 


Fig. 6 
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TT- NEUTRON 
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Distributions of the angles in the C-system between the 7*t-meson and the 
proton and the 7+-meson and the neutron, from the reaction 


Papa 4 Pre D S Pd et 
The full line curve represents the form of the distribution calculated 
from purely statistical considerations. 


When the assumption is made that the 7-meson and one of the nucleons 
result from the disintegration of an excited nucleon it is of interest to 
examine the energy released in the disintegration (Q-value), since 
experiments on the scattering of 7 *-mesons by protons suggest that this 
quantity should be ~160 mev. The experimental eer buen obtained 
for the Q-values together with that predicted by statistical theory (Block 
1956) are shown in fig. 7. It can be seen that if an excited nucleon is 
formed, factors other than those occurring in the scattering of 7 *-mesons 
on protons influence the disintegration of this state since the experimental 
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distribution peaks at about 110 mev. The corresponding distribution 
for 7+-meson and neutron does not show a peak ; this is to be expected 
since no correlation was observed in the spatial angle between these 
particles (fig. 6). If the existence of an excited nucleon is assumed, then 
the angular distribution of the 7-meson and proton in its coordinate 
system may be obtained ; this is displayed in fig. 8. The mesons have 
a tendency to be emitted in the backward direction in this coordinate 
system. In this context we mean the backward direction with respect 


Fig. 7 
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Distribution of Q-values for the decay into a w+-meson and a proton and 
am*-meson and a neutron, of an excited nucleon in the reaction 
p+p > n+N* + 7t+p-+n. 
The curves represent the distribution anticipated from purely statistical 
considerations when no excited nucleon is formed. 


to the direction of motion of the excited nucleon in the C-system. Thus 
they tend to be emitted with low energy in the system of the two colliding 
protons while the protons preferentially appear with the maximum energy 
consistent with the conservation of energy and momentum. The curves 
superimposed on the histogram are the distributions obtained by Ashkin 
et al. (1956) at 150 mev and by Ferrari et al. (1956) at 124 mev in experi- 
ments on the scattering of 7+-mesons by protons. Since the distribution 
of ()-values shows a peak at 100-120 mev, the meson distribution should 
be compared with the results obtained at the lower energy. The 
comparison is surprising. The data given in fig. 8 indicates that meson 
production does not result from the decay of a pure isobar state since 
this would lead to a symmetrical angular distribution about 90° in the 
coordinate system of the excited nucleon. However, the distribution 
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would be modified by the 
presence of any meson pr Ki 
place through state (b). a "ae fae 

For the production of 7°-mesons the numbers available were too few 
to provide useful information on angular or momentum distributions. 


Fig. 8 


Angular distribution of the 7*+-meson from the reaction 
p+p > n+N* + 7++p-+n 
in the coordindate system of the excited nucleon. The curves super- 
imposed on the distribution represent the measured angular distributions 
in the C-system, of 7+-mesons from the scattering of mesons by protons 


at 124 and 150 Mey. 


§ 5. CONCLUSIONS 


The results show that the production of a single 7+-meson with a neutron 
and a proton predominates over all other processes and that the production 
of two 7-mesons rarely occurs in proton—proton collisions at 925 Mev. 

Whilst the topics discussed in § 4 indicate that factors other than the 
availability of momentum space influence the production of 7*-mesons, 
the predictions for the excited nucleon model are not sufficiently unique 
to conclude that firm experimental evidence exists for the formation of 
a 3/2, 3/2 state. A disturbing feature of the experimental data on 
m-meson production in the energy range from 600 mey to 1-5 Bev is the 
lack of agreement between this experiment and that of Block et al. at 
800 mev. Where comparisons can be made our results at 925 Mev appear 
to be in general agreement with those at 650 mev and 1:5 Bev but not 
that at 800mev. In all the experiments except the counter work at 
at 660 mev the statistical accuracy of the results is low, and the marked 
difference between the resujts at 800 mev and ours at 925 Mev could 
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conceivably represent extreme statistical fluctuations. The recent 
evolution of the “hydrogen bubble chamber with its considerable 
advantages in the rate of collection of experimental data, should now make 
possible a more satisfactory study of meson production when the energy 
of the primary proton is in the region of 1 Bev. 
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APPENDIX 
Procedure for Momentum and Energy Balance in Single Meson Production 


If assumptions are made regarding the identity of a secondary particle, 
three values for its energy and momentum may be found corresponding 
to the observed grain density, and its maximum error (1-5 times the 
standard deviation in this instance). 

Since the energy of the incident proton is known it is possible, using 
the principle of the conservation of energy, to obtain nine values for the 
energy of the neutral particle in the laboratory system by combining 
the possible values for the energy of one charged secondary with those of 

the other. Nine values of the momentum and angle of emission of the 
- neutral particle were determined in a similar manner by an application 
of the principle of the conservation of momentum. Thus a further nine 
values can be obtained for the energy of the neutral particle. 

In a three dimensional system with axes corresponding to the laboratory 
energies of the three secondary particles, the nine values for the energy 
of the neutral particle, obtained from energy conservation, lie on three 
parallel, coplanar and straight lines. The nine values obtained from 
momentum balance can be joined by three curves corresponding to the 
three straight lines obtained from energy conservation. If the curves 
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and the straight lines intersect then points on the line of intersection 
represent energies of the secondary particles for which momentum and 
energy are conserved for the measured angles of emission. 

In fig. 9 the three straight lines, and the line of intersection LMN are 
shown ; for clarity the nine points obtained from momentum balance 
are omitted. 

It was found in practice that the point of maximum probability on the 
line of intersection always lay on one of the straight lines, DOE or GOB 
in fig. 9; itis point Min thisinstance. Thus the ratio OM/OD is a measure 
of the probability that energy and momentum are conserved for the event. 


Fig. 9 


Illustration of the energy and momentum balancing technique. 


In practice this procedure was modified by projecting on to a plane 
containing the axes representing the energy of the neutral particle and 
one of the charged particles. The ratio OM/OD is unaffected by the 
change to two dimensions. In order to obtain an acceptance criterion 
related to the measured quantities we define a quantity R : 

@u-G| 
SAG 


Gy is the measured grain density for one secondary track, AGy is the 
maximum error of this quantity and G is the grain density which the 
track must have if energy and momentum are conserved. Thus Cu 
is the grain density corresponding to an energy represented by the point 
O. Similarly G corresponds to the point M and 4G y corresponds to the 
length OD. The energies of the particles corresponding to the point M 


R 
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are those used to obtain the ‘ consistency ’ values of quantities in the 
C-system for which the distributions are plotted in the figures. 

If the variation of the three values of energy of one of the charged 
secondary particles was small compared with that of the other, only the 
central value was taken and thus the nine values obtained for the energy 
of the neutral particle were reduced to three, so simplifying the procedure. 
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SUMMARY 

Quasi-elastic collisions of protons with protons bound in nuclei have 
been examined in G5 emulsions exposed to 925 Mev protons from the 
Birmingham synchrotron. 

The directions of the secondary particles have been used to obtain 
values for two components of the momentum of the bound protons. 
The observed momentum distribution may be described approximately 
by the function exp (—K/7)+0-05 exp (—K/40), where K is the energy 
of the bound proton in Mey. The results show that (2-2+0-3)% of all 
nuclear reactions in the emulsion are quasi-elastic collisions with angle 
greater than 30° in the C-system of the two protons. This figure agrees 
with a calculation based on a semi-transparent nuclear model. The 
calculation indicates that such events contribute about 5°%, of the inter- 
actions with the light nuclei and about 0-7°/ of those with the heavy 
~ nuclei in the emulsion. 


§ 1. INTRODUCTION 


THE search for elastic proton—proton collisions in G5 emulsion exposed 
to 925 mev protons from the Birmingham synchrotron has been described 
by Duke e¢ al. (1957). In the course of this search, measurements were 
made on events with only two emergent particles. Those in which the 
primary and secondary tracks were coplanar and at suitably related 
angles were considered to be elastic collisions with free protons. The 
events which lay close to, but definitely outside, the experimental limits 
for elastic scatterings were interpreted as due mainly to quasi-elastic 
collisions of protons with protons bound in nuclei. “- 

In analysing these events we have assumed that since the binding 
energy is very small compared with the kinetic energy of the incident 
proton, the only difference from a free elastic collision is that the struck 
particle may have appreciable momentum. We have shown that the 
deflection of the protons while leaving the field of the nucleus is 
unimportant in most cases. 


+ Communicated by the Authors. 
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Quasi-elastic collisions of 340 mev protons in light nuclei have previously 
been studied by Wilcox and Moyer (1955), and those of 1-4 Bev 7-mesons 
in nuclear emulsions have been observed by Walker and Crussard (1955). 


§ 2. EXPERIMENTAL PROCEDURE 


During the search for elastic proton—proton collisions, the angles and 
approximate grain densities of the tracks were measured for all disinte- 
grations in which two charged particles were emitted (2p stars) provided : 


1. Neither track was identified as due to an «-particle or heavier 
nucleus. 


2. No sign of an additional track, however short, was seen. 
For the purpose of the present analysis, only those 2p stars which 
satisfied three further criteria were considered : 


3. The projections of the two tracks on the plane perpendicular to the 
direction of the incident particle made an angle greater than 90°. 


4, The event was not identified as due to the reaction p+p—7*-+d. 


5. Neither secondary track was identified from its scattering as due 
to a meson. 


Fig. 1 


Incident track 


The geometry of a 2p star, showing the angles «,, a, and ¢, and the x-, y- and 
2-AXeS. 


In those events which satisfied the above conditions, the measured 
values of the horizontal and vertical projections of the directions of the 
incident and emergent particles were used to give the angles shown in 
fig. 1, «1, w and d. «, and «, are the angles between the secondary 
tracks and the direction of the primary. The formula for ¢, the angle 
between the incident direction and the plane containing the emergent 
tracks, is quoted by Duke e¢ al. (1957), Exact formulae were used 


except when a graphical method employing stereographic projection 
was accurate enough. 
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§ 3. CALCULATIONS 


Separate calculations have been made of the effect of the three 
components of the momentum of the struck proton on the direction of 
the emergent protons. 

If the three tracks are not coplanar, this may be attributed to a 
component of momentum of the struck proton, perpendicular to the 
incident proton direction and in the plane perpendicular to that of the 


two emergent particles, equal to p,=pine tan ?; Pine is the momentum 
of the incident particle. 


Fig. 2 
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The values of «, and «, for all 2p events which satisfy the criteria in § 
a 4 less Phan ie i the dotted curves show the angles calculated for elastic 
p-p events and for quasi-elastic events with p,—=+0-12pinc ; the me 
corresponding to angles ae of 30° and 150° for the two protons in their 
centre-of-mass system is drawn also. Open circles represent events 
considered to be elastic. 
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event was obtained by interpolation between curves of this type, or 
where p, is small by the approximate formula 
TI] BT) 
Pa~Pine tandy, + tama, + (HE —2)0 + HT? * 
T=tan«,tana,, and F is the total energy of the system measured in 
units of the rest-mass of the proton. 


Fig. 3 
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Grain density of individual tracks in 2p stars with ¢<7° and |p,|<0-12pine ; 
ae is the centre-of-mass angle obtained by assuming that the event is a 
quasi-elastic collision. The full curve gives the grain density for elastic 
collisions, while the broken curves give values for quasi-elastic collisions 
with p,=-+0-12pince. Open circles refer to events rejected because the 
grain density of at least one track falls outside the limits p,=+0-24pine. 


A component p,, perpendicular to the incident proton direction, 
and in the plane of the two emergent protons, tends to increase one 
angle and decrease the other. Within the region of C-system angle 
considered, the effect can be distinguished from those of other components 
only by measurements of the energies of the secondary particles. To 
a good approximation the grain densities of the secondaries of a quasi- 
elastic collision depend only on p, and on their angles in the C-system ; 
fig. 3 shows the curves for quasi-elastic collisions with zero p,, and for p, 
corresponding to an energy of 20 Mev. In principle the distribution of 
p, could be obtained from the energies, but experimental errors in the 
grain densities are too large for this measurement to be significant. 
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The C-system for a particular quasi-elastic collision is the system in 
which the centre of mass of the incident proton and the moving struck 
proton is stationary. The angle a. of emission of two protons in this 
system may be obtained from their directions &, and «, in the laboratory 
system, though an unknown component P, is liable to give serious uncer- 
tainty in the assigned value of a when «.<30°. In practice, we have 
plotted lines of constant «, on diagrams resembling fig. 2, and obtained 
values of x for individual events by interpolation. 


§ 4. TREATMENT OF DaTa 


The problem is to separate quasi-elastic events from the background 
of other nuclear events which by chance have also produced two tracks : 
it is similar to the problem of separating the elastic p-p events from 
those now being considered. 

In fig. 2 are plotted the values of «, and «, for all 2p events satisfying 
the criteria given in § 2, and with ¢ less than 7°, ie. |p, |}<0-12pine. 
Elastic p—p events lie within experimental error of the line 


tan «, tan «,.—2/H, 


and quasi-elastic events lie close to it. The rest, due to other processes, 
have angles which are not closely correlated, though small values are 
rather more probable. The majority of quasi-elastic events are expected 
to lie between the lines corresponding to p,=+0-12pine. 

Figure 2 shows also a line of constant angle «, in the C-system: a is 
30° for one proton, 150° for the other. Provided «¢ is over 30°, the values 
of «, and «, give a satisfactory measure of p, which is almost independent 
of the existence of a component p,. When «¢ is less than 30° (or greater 
than 150°), the apparent value of p, is liable to be much affected by 
experimental error in «, Or %, by the existence of a component p,, and 
by the greater probability that the slower proton, with energy of the order 
of 100 mev, has suffered further deflection in crossing the potential 
barrier of the nucleus. Events which could correspond to quasi-elastic 
collisions with c.m.s. angle <30° (or >150°) are therefore not considered. 

The grain density of each track of a 2p star is plotted as a function of 
ae: fig. 3 shows this plot for stars with 6<7° and | p, |<0-12pine. If the 
points corresponding to both tracks of an event do not both lie within 
experimental error of the limits p,—=-+0-24pine the vo is rejected. 

For the remaining events, with #,>30° and ¢< 13-5 , the values of p, 
are obtained from the angles «, and «, by interpolation on exact curves, 
or from the approximate formula given above. They are plotted in 
fig. 4 against the corresponding values of p, obtained frora the igure 
from coplanarity. The semicircle encloses all quasi-elastic events 


with energy of the struck proton less than 20 Mev. The ellipse corresponds 
mits which are set as bounds for elastic collisions 


approximately to the li 
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with free protons by Duke ef al. (1957). The axes of this ellipse give a 
rough measure of the limits of experimental error for each point on the 
diagram. 

§ 5. RESULTS 

The distribution of momentum of the struck protons may be obtained 
from fig. 4. Circles are drawn on the p,-p, diagram, and the number 
of points per unit area is plotted in fig. 5 as a function of the proton 
energy obtained from the resultant momentum in the xz plane. 


Fig. 4 


Pine 
The values of p, and p, for events with ae>30°. The semicircle encloses all 
events with energy of the struck proton less than 20Mev. The ellipse 
encloses the elastic events, which are not plotted individually, and quasi- 
elastic events estimated to number four. 


Apart from the quasi-elastic events, there is a background in the 
p,-p, plot which is thought to consist mainly of events in which a meson 
has been produced but not recognized, and a few cases where one of the 
protons, quasi-elastically scattered in a heavy nucleus, has suffered 
further deflection before emerging; almost all cases where a scattered 
nucleon has formed a deuteron by pick-up are excluded by the grain 
density criterion. In the background events there is little likelihood 
of coplanarity or angular correlation; to estimate their effect on fig. 5 
we have assumed that they include half the events whose equivalent 
resultant momentum is between 0-15pine and 0-24pine (corresponding to 
energy between 30 Mev and 80mev) and that they are distributed 
uniformly over p, and p,. A similar result is obtained if we assume that 
all the events beyond 0-24pine are background. The corrected numbers 
are shown as open circles in fig. 5, 
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a the momentum distribution is Gaussian, it may be described by 
a function N(p,) which is the probability per unit p, of a given value of 
the component p,, and ‘ 


N(p.)dp,=exp (—p,2/p,2)dp» «. . . . « (1) 


Fig. 5 


Number of events per unit area 


0 10 20 30 40 «50 60 70 
Energy K, in Mev. 


The number of points per unit area of fig. 4 plotted in arbitrary units against 
the energy of the struck proton. The open circles represent numbers 
corrected for background, and the full circles the uncorrected numbers ; 
the probable errors of the latter are the same as those shown for the former. 
The full line is obtained from a probability per unit volume in momentum 
space proportional to exp (—K/7)+-0-05 exp (—K/40), while the broken 
line represents a Gaussian distribution falling to 1/e at Ky=11 Mev. The 
probable errors of K for individual events are about } of the horizontal 


spacing between the points plotted. 


Also the probability per unit area of a p,—p, plot is the same function of 
the resultant momentum 4/(p,’+?.") : 
N[V(p.2-+p2) dp dp.=exp [—(p.? +P; )/Po tp dp. + + (2) 
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and the probability per unit volume in momentum space is the same 
function of the resultant momentum | p |: 


N(| p |)dp, dp, dp,=exp [—(p2? +P? +D2")/Po lap « Ip, Ip, 
=exp (—K/K,)dp, dp, dp, 0 eee) 


where K is the kinetic energy of the struck proton. 

Fitting an equation of form (2) to the corrected momentum distribution 
shown in fig. 5, we find that p, is a momentum corresponding to an energy 
K,=11 Mev, with an uncertainty of +3 Mev resulting mainly from the 
correction for background events. The r.m.s. value of a single component 
corresponds to 5-5 Mev, which gives a mean nucleon energy of 16-5 Mev. 
But there is evidence that the distribution cannot be fitted by any 
single Gaussian curve (see § 6). 

If we apply a less rigorous grain density criterion, there is no significant 
change in the momentum distribution obtained, but a somewhat greater 
uncertainty in the background. 

If we take distributions of p, and p, separately from fig. 4, and fit 
equations of type (1) to them, we obtain values of py which agree well 
with each other and with that given by the number per unit area. This 
gives extra confidence in the reliability of the results and the analysis. 


§ 6. DISCUSSION 


Henley (1952) found that of several possibilities a Gaussian momentum 
distribution gave the best agreement with experimental work on 7-meson 
production from carbon by 340 Mev protons (Richman and Wilcox 1950). 
The width of this distribution was chosen to give an average kinetic 
energy of 19-3 Mev, to fit the results of Hadley and York (1950) on 
deuterons formed by the pick-up of protons in carbon nuclei by 90 Mev 
neutrons (Chew and Goldberger 1950). Selove (1956) has examined 
in greater detail the pick-up of neutrons in carbon and beryllium by 
95 Mev protons : the momentum distributions derived from these results 
in the same way appears to be the sum of two Gaussian distributions 
described by exp (—K/7)+0-15 exp (—K/50). 

Cladis ef al. (1952) examined the energy spectrum of single protons, 
incident at 340 Mev, scattered quasi-elastically from carbon and oxygen 
at 40°; they found the best fit to be given by a Gaussian distribution of 
internal momentum, falling to 1/e at a momentum corresponding to about 
16 Mey : the average kinetic energy of the nucleons was therefore about 
24mev. Wilcox and Moyer (1955) detected in coincidence both protons 
from quasi-elastic collisions of 340 Mev protons with protons in various 
target nuclei, and found for beryllium a Gaussian momentum density 
distribution falling to l/e at about 20 Mev; this gives a mean nucleon 
energy of 30 Mev. 

Wattenberg et al. (1956) have fitted the nucleon momentum distributions 
in carbon and oxygen, which they obtained from experiments on the 
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photodisintegration of quasi-deuterons, to Gaussian distributions falling 
to 1/e at 19 Mev. 

If our results are to be fitted to a single Gaussian distribution, this 
must fall to 1/e at 1143 Mev. The distribution found by Selove fits 
_ our results rather less well than does the single Gaussian distribution. 
A better fit, however, is given by a distribution of the form 

exp (— K/7)+-0-05 exp (— K/40). 
It should be noted that this function gives the probability per unit volume 
im momentum space: since it is not a simple Gaussian function, the 
probability per unit area of a p,—p, plot is not given by the same 
expression ; the curve shown in fig. 5 has been obtained by the appropriate 
transformation. 

It is difficult to estimate the background correction exactly: if, as 
an extreme course, we make no correction whatever for background 
events, the best single Gaussian distribution falls to 1/e at about 18 mev, 
and a good fit with two Gaussian distributions is given by 

exp (—K/7)+0-07 exp (—K/50). 

According to calculations based on a semi-transparent nucleus model 
(see Appendix), about 30° of the observed quasi-elastic events should 
come from the heavy nuclei in the emulsion. Thus, although the 
momentum distribution will refer mainly to protons in light nuclei, 
heavy nuclei will make a significant contribution and may be responsible 
for the large proportion of protons with low momentum. 

The mean free path in emulsion for quasi-elastic collisions with a >30° 
is 16+2m, compared with 0-35 m for all nuclear interactions. Thus 
(2:2+0-3)% of all events are observed as quasi-elastic events with angle 
greater than 30° in their centre-of-mass system. 

Scanning efficiency has been discussed by Duke et al. (1957): their 
results show that the chief loss is among the events which we have 
excluded by the criterion %.>30°. The loss of events with a. >30° is 
slight, and does not cause any significant uncertainty in the above figures ; 
the form of the momentum distribution is quite unaffected, since there 
is little correlation between the probability of loss of an event and the 
values of p, or /,. 

The Riecuions outlined in the Appendix give 1-7°% as the probability 
that an event is a quasi-elastic collision with a single prcien, pel both 
protons escaping from the nucleus at angles between 30° and 150° in the 
C-system. In view of the approximations involved, the agreement 
between this value and that observed is satisfactory. The calculations 
indicate that quasi-elastic collisions contribute about 5% of the 
interactions with light nuclei and about 0-7° of the interactions with 


heavy nuclei in the emulsion. 
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APPEND LX 
Calculation of the Cross Section for Quasi-Elastic Events 
By J. G. McEwEn and J. B. A. ENGLAND 


The nucleus is assumed to behave like a uniform sphere in which the 
probability of a moving proton colliding with a nucleon is given by 
a constant coefficient as is the absorption of light in a semi-transparent 
medium. The probability of a particular quasi-elastic event being 
observed is the probability of the initial collision multiplied by the 
probability that each particle will escape without further collision. 
Calculations have been made for different types of collision at different 
points in the nucleus and mean values for the whole nucleus have been 
obtained by numerical integration. The procedure is shorter than the 
Monte Carlo method in this application, and could be modified to deal 
with other nuclear reactions at high energies in which few particles are 
involved. 

Figure 6 represents a plane in a nucleus in which a quasi-elastic collision 
is assumed to take place. The incident proton travels a distance 1, 
in the nucleus before making a collision, and the secondary protons travel 
distances /, and 1, before emerging. The probability that such a collision 
will take place when the primary has travelled a distance between 1, 
and 1,+dl,, and that both secondaries will escape, is 

kydl, exp (—kyl,) exp (—kgly) exp (—keols) 
where k,, k, and ks are the absorption coefficients in nuclear matter for 
the primary and secondaries respectively. It is required to integrate 
this expression over all parts of the nucleus where collisions could occur, 
the directions of the secondary particles having the appropriate angular 
distributions, and the plane of the primary and secondaries at any point 
being allowed to take all orientations. 
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For simplicity quasi-elastic events are treated as being coplanar and 
having the same angular distribution inside the nucleus as free proton— 
proton collisions. Since only events with C-system angle %.>30° were 
observed, the angular distribution was divided into four intervals of equal 
probability between «.—30° and %—=90°, and representative angles 
a, and «, for each interval were used in the calculation. 


Fig. 6 


Direction of 
incident protons 


> 


Diagram of a plane in a nucleus, containing the paths of an incident proton 
(length /,) and two secondary protons (lengths J, and /,). 


The cross section of the nucleus projected on a plane normal to the 
direction of the incident protons is divided into five zones of equal area 
as shown in fig. 7(a). The dotted circles divide each zone in half. 
Representative points are taken, one on each of these circles, and protons 
are considered incident at each point. The paths of these protons through 
the nucleus are shown as lines in fig. 7 (0). Pras 

Along each path points are taken representing intervals of equal 
collision probability (calculated from the absorption coefficient). The 
secondaries from collisions along a given path can lie with equal probability 
in any plane containing the path, but the probability that the secondary 
particles will emerge from the nucleus without further interaction 
varies with the position of the collision and with the angle % pe 
the plane of the collision and the plane through the centre of the nucleus, 
i own in fig. 7 (0). . 
ae aa fo eect : particular path, the path ae in ae 
nucleus J, and /, of the secondary particles from collisions along : a 
path are measured. The collisions are considered to take place at the 
representative points mentioned above, unless one of them lies very near 


the nuclear surface. 
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In practice lines representing the directions of the primary and 
secondaries of each type of event are drawn on tracing paper and placed 
on the appropriate point on the proton path in the given plane as in 
fig. 6. Since paths in only half the plane in fig. 7 (6) are being considered, 
an equal number of events in which the angles are reversed must be 
used. Symmetry then limits 5 to values between 0 and 37. 


Fig. 7 


(a) Cross section of a nucleus, perpendicular to the incident protons, divided into 
five zones of equal area. (6) Paths through the nucleus of five incident 
protons, one representative of each zone. 


The whole set of collisions is not used at each point. To keep the 
calculation reasonably short and at the same time to give it reasonable 
accuracy, two events were considered at each point, one with normal 
angles, and one with reversed. ‘These were selected at random, within 
certain restrictions imposed to avoid extreme fluctuations. 

The probability of escape is calculated for each collision with the 
absorption coefficient appropriate to each track. In the actual calcula- 
tions the absorption coefficient always occurs multiplied by the nuclear 
radius. The value of the product used for the incident proton was that 
found at 860 Mev by Chen e¢ al. (1955), viz. k;R=0-70A1%. The values 
for the secondary protons in each representative event were calculated 
from the known proton—free nucleon total cross sections at lower energies, 
with an allowance for the suppression by the Pauli principle of collisions 
with low momentum transfer (Goldberger 1948). 

The average probability taken over all events along a given proton 
path in one plane varies with the angle % of the plane. The variation 
is small for paths near the centre of the nucleus but large and rapid for 
those near the edge. In order to obtain a reliable average probability 
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over all events, planes and angles, more planes were considered for the 
outer paths than for the inner. 

For each path the average probability per collision in each plane is 
plotted against the angle 4; twice the area under this curve measures 
the contribution from collisions along this path. The total probability 
of escape per collision of an incident proton in the nucleus is the average 
of contributions from all the paths. This total was obtained for light and 
heavy elements separately. 

The accuracy obtained with this method is limited by the failure of 
the uniform model for paths too near the edge of the nucleus. The 
calculations for this region give a rather high probability of quasi-elastic 
collisions. This probability will be reduced by several effects : there is 
believed to be a shortage of protons near the nuclear surface (Hess and 
_ Moyer 1956); also secondary particles produced near the surface are 
especially likely to approach the surface at large angles of incidence, 
and to suffer internal reflection or refraction by the potential barrier. 

To compare the results with experiment we allow for the fact that only 
events with «.>30° were identified, by multiplying by a factor which 
we estimate to be 0-63. This figure includes an allowance for the 
suppression by the Pauli principle of collisions with low momentum 
transfer. Further, since only collisions in which the struck particle is 
a proton are observed as quasi-elastic, we multiply by Z/A. The 
measured values of total (elastic and inelastic) and of elastic cross sections 
for the collisions of protons with free protons at this energy (Chen et al. 
1956, Smith et al. 1955, and Duke et al. 1957) show that mesons are 
produced in two-thirds of these interactions. With the assumption that 
this is also true when the struck proton is bound, the calculations give the 
probability that an event is a quasi-elastic collision with a single proton, 
with both protons escaping from the nucleus at angles between 30° and 
150° in the C-system, as 5% in light nuclei and as 0-7°% in heavy nuclei. 

With the known composition of the emulsion and the absorption cross 
sections obtained by interpolating from the work of Chen et al. (1955) 
the results for light and heavy nuclei combine to give the total Insel aS 
for quasi-elastic events in emulsion : this is 1-7°/, in the range of C-system 
angle considered. 
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ABSTRACT 


A comparison is made between the theories of Wigner (1934), Wigner 
and Seitz (1934), and of Pines (1953, 1955), Bohm and Pines (1951, 1953), 
Pines and Bohm (1952), which both represent different approaches to an 
inclusion of correlation in the wave mechanical treatment of a free electron 
gas. Starting with a determinantal wave function for the electron gas, 
correlation is taken into account by means of a ‘ modulating ’ function 
by which the complete wave function must be multiplied. It can be 
shown that in both theories these modulating functions are nearly the 
same. A more detailed consideration of Wigner’s treatment gives a 
‘cut-off’ wavelength for the collective behaviour which is numerically 
nearly equal to the ‘cut-off > wavelength in the theory of Bohm and 
Pines. 


§ 1. INTRODUCTION 
TuHE Hamiltonian for a gas of N electrons is given by the expression 


H= ¥ (—#)2m. A)+1/2. y Coley acre ae sae) 


t=1 1,j=1 

where the first term denotes the kinetic operator and the second term 
denotes the Coulomb interaction of the electrons. The summation in the 
second term is over all pairs (7,7) of indices, and the prime indicates that all 
pairs with i=) are excluded. If we build up an approximate eigenfunc- 
tion corresponding to the Hamiltonian (1) by taking a product of plane 
waves y;, 

Pe Sie cee (9) 10,1 6s) ato te eee) 


orja determinant of plane waves 

P= etl ig, ass YN} > + 9 0. o 0.18, (9) 
W,, and Y,;"are known as the solutions of the Hartree (H) or Hartree— 
Fock (HF) approximation. With the approximations (2) or (3), the one- 
electron functions y, have to satisfy the Hartree or Hartree-Fock 
equations respectively and for the Hamiltonian (1) a set of plane waves is 


nicated by the Author. : 
Be euent es - Max Planck Institut fiir Hisenforschung, Diisseldorf, 


West Germany. 
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known to give a self-consistent solution of these equations. Assuming 
doubly occupied levels, one for each spin, then no more than two electrons 
can have identical wave functions, i.e. the same wave vector for the plane 
waves, and those which have identical wave vectors belong to electrons 
of opposite spin. 

Obviously the Hartree function Y%, does not include correlation, as the 
wave function of the ith electron is 4, and therefore independent of the 
positions of the other electrons. Some correlation is brought about by 
using the Hartree-Fock function Yy,, but these correlations are purely 
accidental, originating from the use of the determinantal eigenfunction 
and not from a proper inclusion of the Coulomb interaction. As a 
determinant is the simplest antisymmetric combination of one-electron 
functions, it can also be said that these accidental correlations arise from 
the inclusion of the Pauli principle. 

That electron correlation plays an important réle is immediately seen 
by calculating the lowest eigenvalues for a free electron gas using the 
Hartree and Hartree-Fock approximations. In the Hartree approxima- 
tion, this energy is given in atomic units by 


Hy=2-21/r?24+-1-2/r, Ee es ES 
and in the Hartree-Fock approximation 
Egp=221/r2+1-2/r,—0-916/r,,  . . . . « (5) 


Here 7, is the interelectronic spacing defined by J . (47/3) . a)? . r3=L?, 
where a, is the Bohr radius and L the side of the cubical box in which the 
electrons are included. The first and second term in (4) and (5) are the 
kinetic energy and the potential energy in the field of the smeared-out 
charge of the electrons, while the additional term in (5) is the so-called 
‘exchange-energy ’. Both approximations give far too high energies. 
According to the Hartree energy (4) no cohesion would exist in the alkali 
metals, and using the Hartree-Fock energy (5) the cohesion in these 
metals would be much weaker than actually observed. Therefore the 
inclusion of correlation should improve the wave function and con- 
sequently the energy. It has become a common usage to define the 
correlation energy as an energy gain compared to the Hartree-Fock 
energy, so that (5) has to be replaced by 


H=2-21/r2+ 1-2/r,—0-916/1,—E corr. 2s eer Oo 


where £,,,,, denotes the correlation energy. 


§2. IMPROVEMENT OF THE Wave FUNCTION 
Starting with an approximation of the N-electron wave function made 
up by a combination of one-electron functions, for example the Hartree 
function WY, or Hartree-Fock function Yyp, it is easily seen how an im- 
provement can be achieved. We discuss the question using the Hartree 
function Yy. As in this approximation the probability density of finding 
two electrons (7,j) a distance | x,—x;, | apart is completely independent of 
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the distance | x;—x; |, an improvement of the Hartree function could 
obviously be achieved by modulating Y, with a function f(|x;—%; |) 
which has the following properties : 


(a) f(|x;—x; |) is positive for all values of |x;—x, |. 
(6) f(|x;—x;|) is smallest for |x,—x,|=0, does not decrease with 


_ increasing | x;—x,| and approaches a constant value for large | x,—x, ]. 

Such a function reduces the probability density when two electrons 
are close together and does not affect the probability density (excluding 
a normalization factor) when the electrons are far apart. Furthermore 
it can be seen that the potential energy of the electron gas is decreased 
because f(| x;—x; |) decreases when 1/|x,—x,| increases. The suggested 
improvement of the Hartree function then would be 


N 
Je allt (|=; |) 
4j= 
where the prime indicates that i=) is excluded. The same is true if we 
use the Hartree—Fock function Wy, instead of Py. 
To determine the modulating function 


N 
F(x,, os Millers) ty ice EO) 


4, 
the methods of Bohm and Pines and of Wigner have been used up till 
now. They will be dealt with more thoroughly in the next sections. 
Wigner does not give the modulating function explicitly but it can be 
derived from his results. Another attempt has been made by Mayer 
(1955) using the density matrix method but his results cannot be compared 
numerically with the other treatments. 


§3. Tue TREATMENT OF BOHM AND PINES 


The method of these authors (Bohm and Pines 1953, Pines 1953) 
can be shortly summarized as follows : By a series of canonical transforma- 
tions, the Hamiltonian (1) is transformed into the sum of three Hamil- 


tonians : 
Via Vs Be Sie Na od ees . . . ° e . ° (7) 


where Han, is defined by 
N 
H —f2/2m. > A(1—N’/3N)—27Ne? > 1/k? 
i=1 


part. = k<hig 


and describes the behaviour of independent particles, and where 
I6 *k f * 
Hon = > hio/2 . (A, A ptA rank ) 
k<ke 
describes collective oscillations in the collective coordinates A;, A,, 
"res o a collection of individual electrons inter- 
and where H,,. corresponds t 


acting with a weak short-range force (s.r.=short range), 80 that H,,, can be 
neglected in a first approximation, For a complete definition for the 
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letters occurring above see Bohm and Pines (1953). , denotes the absolute 
value of a wave vector smaller than ky, where ky is the absolute value of 
the wave vectors at the top of the Fermi-distribution. That means that 
there are just as many collective degrees of freedom as there are wave 
vectors smaller than &,. The frequency of the oscillations to a first 
approximation is given by 
wrwp=(4rNe?/m) 42. 

Normally, these collective oscillations are not excited as their energy is 
far higher than the energy of an electron at the top of the Fermi- 
distribution, therefore the collective part of the Hamiltonian enters 
merely as a zero-point energy. Introducing the dimensionless parameter 


B=k,/ ko 
which is related to the number of collective degrees of freedom N’ and 
the number of electrons NV by 
f/6=N'/3N, 

Bohm and Pines determine this parameter by minimizing the total energy 
obtained by the splitting-off of the total Hamiltonian into an individual 
and a collective part. 

Under these assumptions, neglecting the weak short-range part H,, of 
the Hamiltonian and taking subsidiary conditions imposed on the com- 


plete wave function into account, the complete wave function is found to 
be 


N 
Y=Wap. []’ exp{—g(x,—x,)/2hwp}, . . . . . (8) 


4,j9=1 


where g(x;— x,) is given by 
g(X;—X,)=27e?- > exp (tk(x,—x,))/h? 2°. . . (9) 
!Ki<ke 


and therefore has just the same form as we have suggested in § 2. 
We wish now to evaluate the modulating factor 


f(x;—x;)=exp (—9(x;—x;)/2fiwp) 
which enters the wave function (8). The argument of the exponential 
function can be transformed into a sine integral and we get 


f(*;—x;)= exp {— (€/2hwp|x;—x,|) . (2/r7) . St (k, |x,—x; l)} 


where Si(x) is defined by 


> 


Si Gels (sin y/y)dy. 
By setting k 
| x, —x; |=. 75; 
where ay is the Bohr radius and introducing the values for the constants 
we finally obtain 


F(ris)= exp {—(2) «(GM « (r.27/rg) « (2/ar) . Si((97/4)2? . B . (r,s/r,))}. (10) 
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This modulating factor has the properties (a) and (b) stated in the preceding 
sections and plots of f(r,;) for the interelectronic spacings r,=2 and r,=4 
are given in figs. 1 and 2. These figures also contain the plots for the 
modulating functions derived from Wigner’s investigations and we will 
deal with them in §4. The values 7,=2 and r,—4 are very close to 
the actual electronic densities found in Al and Na. 


Fig. 1 


j j , Si) 
The modulating function f(7;;) for 7;=2. 
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The modulating function f(7 4) for r= 4. 
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Tt is interesting to note what values the modulating function f(r;;) 
takes for 7,,=0. Calculating these values from (10) we find 

f(O)\=6sp 038287772 Bes ee cen) 
The quantity 6 giving the ‘ cut-off’ wavelength k, for the collective 
coordinates depends on the interelectronic spacing r,. Taking the 
numerical values of 8 given by Bohm and Pines, values of (11) have been 
calculated in the range 7,=1 to 7,=5. The results are given in table 1 


and fig. 3. 
Table 1. Values of f(7,;) for 7,;=0 derived from Bohm and Pines’ results 


~ 
ne) 
= 
= 


Sooo 
OD S180 
lon on en) 


OF Whe 


O9 


0:8 
Bohm and Pines 


Or7 


0:6 


Om 
| 2 %) 4 bY 


»—_ I, —- 


Values of (0). 


From formula (11) and table 1 it can be seen that the value of the 
modulating factor for ,,—0 drops to zero when the electron gas becomes 
very dilute and approaches unity when the density increases. In the 
range of actual metallic densities the factor takes an intermediate value 
between both extremes. For extremely high densities, the modulating 
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factor f(7,;) approaches unit for all values of 7,; in that case electrons behave 
more and more like non-interacting particles. 


$4. Toe TREATMENT oF WIGNER 


The method of Wigner for including correlation is actually based on the 
application of second order Rayleigh-Schrédinger perturbation theory. 
For the wave function without Coulomb interaction Wigner and Seitz 
(1934) starts with a normalized product of two determinants made up 
from plane waves 


Yy(Xq) ~~. hy (X,) Bilyi) -.- daily.) 
1 . 
nt . . . 3 


Pur= (12) 


(4) Ss oS b,(X,) brAY1) ae b(n) 

The first determinant contains all electrons with upward spin and the 
second determinant all with downward spin, so that N—2n is the total 
number of electrons. Such a wave function has been shown to give the 
same energy as a single determinant including spin functions for doubly 
occupied levels (Dellbriick 1930). To include correlation, Wigner (1934) 
assumes that correlation is brought about by the positions of all electrons 
with downward spin influencing every electron with upward spin. That 
means the wave function for the vth electron with upward spin 


b(x)=L-3 . exp (27ivx/L) 
normalized in a box of sidelength L is replaced by a perturbed wave 
function 
UAV ---s¥a3 =H ly; X)=L-8" . (exp (2Qrivx/L) 


+ >o,,- exp (Qrix/L)}=. -. -° . . . (18) 

, . 
where the «,, depend on the positions y,, ... ,y, of all the electrons with 
downward spin and y is an abbreviation for the complete set y,, ... , Yn. 


The major difference in the methods of Bohm and Pines and of Wigner 

lies in the fact that Wigner contents himself with the accidental correlation 

brought about by the determinantal form of the wave function ae 

electrons with parallel spin and introduces correlation only for electrons 
71 i in. cart 

; eee the perturbation method for the ee eee of the 

x,,, in (13) Wigner (1934) finally gets for the perturbed wave function 


‘ bly; x)=H,(x) - 1—h,(yi—¥)— - - - —he(Ya—¥)} (14) 
where h,(y—x) is given by 


exp [27r2(v—yw) (y—x)/L] 
nY—W= > Goo testi, table) 


2v—“ 
oce. 
SAC eeet OA eed ek 
Ee Aree Gti Faced CLE 
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The summation in the first sum is over all wave vectors p. of the Fermi- 
distribution for which 2v—p is an occupied state and in the second sum 
over all w for which 2v—p is unoccupied. The quantities t,,, and t,, are 
defined by 
y= (407h2/mL?) . (u—v) pw, tyy,= (2/87)? . (€?/r,) - e(|u—v| /ko), 
hes 
1 0-415 1-465 
“(= 3a 4032-50) (+o%72b02 (1 -+0%/2-56)8 
and where k, is the absolute value of the wave vectors at the top of the 
Fermi-distribution. As usual, e denotes the electronic charge and r, 
the interelectronic spacing. 

If we replace now the unperturbed wave functions in the first deter- 
minant of the complete wave function (12) by the perturbed ones (14), we 
would not get a wave function which consists of the wave function (12) 
multiplied by a modulating factor, so that a comparison with the Bohm 
and Pines result and the form of our wave function in §2 would not be 
possible. To make this comparison possible we have therefore adopted 
the following procedure: Firstly we have replaced the function 
h,(y—x) in (15) by an averaged mean h,, (y—x). Secondly after inserting 
these averages in (14) and (12), every term in the same column of the first 
determinant in (12) consists of the unperturbed plane wave multiplied 
with the same factor, so that these factors can be taken out of the deter- 
minant. Then the complete wave function gets the form of the un- 
perturbed wave function (12) multiplied by a modulating factor. As 
already mentioned above, the Wigner treatment does not include correla- 
tion of electrons with parallel spi apart from the accidental correlations 
brought about by the use of the determinantal wave function. Therefore 
the modulating factors depend only on the differences | x,—y,,|, where 
x, is the position of the vth electron with upward spin and y,, the position of 
the pth electron with downward spin. 


§ 5. THe AVERAGE COEFFICIENTS AND THE ‘ CUT-OFF’ 
WAVELENGTH 
To find the average value of the function h, (y—x) in (15), averaged 
over all v of the Fermi-distribution, we look at the coefficients of 
exp (271(v—p) (y—x)/L) 
(u—v)? 
in the two sums on the ie side of Pe These coefficients are either 
= (tye tty)? - (e/a) 
if 2y—p is an es sine or 
= (ty thy’ Bie Vu ee eon) . (e/7L) 
when 2v—p is an Na iben ane state. We now sum all the coefficients 


c,, and c,,,’ for which w—y is equal to a fixed value A. The vector A can 
take all values which are smaller than or equal to 2k), where ky is the 
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absolute value of a wave vector at the top of the Fermi-distribution. 
W ithout giving here full details of the rather lengthy, though elementary 
calculations, the summing of the cy, and c,,’ for constant A=.—v 
turns out to give in a good approximation 
1A(B) - (Allo) + (ie) - (LAM) 
ees C,,{=A2 . c= 4 \ 0 16} ° 1™O 2 
=, se i: : ‘ A .[A|"/kg? + €(c) et) 
where <(o) has been written for o2 . €'(c), €'(c) being the function defined 
in §4, where o=|A |/ky and the quantity A is given by 
A==(3a/4) . (2°37 .7)1? .. (kgr,)-*=9-05/k,,. 
Introducing numerical values we get from (16) 
pe ek (| mA |/ko) ie) (| A |/Ao)8 


*~ — (9-05/r,) . (| A|/ko)4 + €(c) 
Returning now to (15) we obtain for the average h,.. (y—x) of the function 
hify—x), averaged over all v ' 


Nay (Y—X)= DC, - exp (QrtA(y—x)/L) . . . . (18) 


(17) 


where the summation runs over all vectors A for which |A| is smaller 
than or equal to 2k,. 

The calculation of the coefficients c, has been carried out for different 
values of the interelectronic spacing 7,. As an example the coefficients for 
the case r,=4 which is very close to the actual interelectronic spacing in 
Na have been written down in table 2 and have been plotted against A in 
fig. 4. 


Table 2. Coefticients c, for r,=4 


A 0-125k, 0:25k, 0:375k, 0-5ko 0-6k 0-7ko 
6q=— 0-94 0-93 0-93 0-99 0-8 0:6 

A= 0-75k,y 0:8k, 09k 1-0ko 1-5k 2-Oko 
C= 0-5 0-4 0-24 0-14 0-01 0 


From these data and formula (17) it is obvious that the coefficients c¢, 
are nearly equal to unity up to a certain wave vector and then tend 
towards zero with increasing A. This means that up to a certain wave 
vector the perturbation measured by the ¢ as remains constant and 
decreases rapidly for larger wave vectors. This behaviour of the coefli- 
cients c, suggests that the summation in (18) should be cut off at a fixed 
wave vector A=A, and the coefficients c, replaced by unity. This yields 

hey x= Sy exp (271A(y—x)b) se. (19) 

lAl<laol 

wavelength 2, will be expressed as a fraction of the 


where the ‘ cut-off ’ RC gO 
wave vector at the top of the Fermi-distribution. 


absolute value ky of a 
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We denote this fraction p’ : 
JA, | =f" . ko, 

In this approximation we have now to determine the cut-off wavelength 
as a function of 7,. We have chosen for A, that value of A for which the 
coefficients c, have decreased to the value 0-5 which obviously is a satis- 
factory assumption. This calculation has been carried out in the range 
of actual metallic densities and has given the values of f’ plotted against 
r,in fig. 5. This figure also contains the values of 8 in the Bohm and Pines 
treatment and it can be seen that the quantities 8 and 6’ which determine 
the ‘cut-off’ wavelength in both treatments have nearly the same 
numerical values. 


Fig. 4 


lO 


O:5 


0:5 ne) 1:5 
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The coefficients c,. 


§6. THe Mopunating Function 
Using the approximation (19) for h,,(y—x) it is now easy to perform 
the summation on the right side of (19). Measuring | y—x| in multiples 
r of the Bohr radius ao, | y—x|=r . a) and by setting r/r,=z we obtain 
hay (¥ —X)= (4/377) . 2-? {2-1 . sin (2a . 0-306 . B’ . z)—27 . 0-306 . B’ 

008 (2ar7.0-306 + Bie) 5 eee ns eee a 20) 
and this result immediately follows from (19) without further approxima- 
tion. This function is positive everywhere, has its largest value for 
v=0 and decreases with increasing r. Remembering now formula (14) 


and the method of our averaging procedure, the modulating function which 
has to be multiplied with the starting wave function (12) and which 
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represents the inclusion of correlation between electrons of opposite spin 
is found to be 


F(x w=TT I—h,, (Yi —,) — ha MY X= oem Navn (21) 


In this form the modulating factor cannot be compared immediately with 
(6) or the Bohm and Pines result (8) ; moreover the perturbation method 
applied by Wigner breaks down when for one x; several of the h,, (y;—x,) 
become large, i.e. when several of the differences | y,;—x,| become simul- 
taneously very small for one (or several) x,;. We therefore confine 
ourselves to the case when for a fixed x; only one of the coordinates y,, 
let us say y,, comes very close to x; while all the others ly,—x,|, ... , 
ly,—*;] are sufficiently large so that the corresponding h,, can be 
neglected. Taking out of the product (21) the term which contains x,, it 
can be written under the restriction just stated (which can be partly 
removed afterwards) 


F(¥1— Xi) =1—- Igy (Y1—%). 2 se e+ (22) 
This expression has to be compared with f(r;;) given in formula (10) which 
has been derived from the result of Bohm and Pines. 


Fig. 5 
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The quantities B and f’. 


In the limit |y,—x,;|~0 it can be easily deduced from (20) that 
(y,—;,) tends to 
OF Wt a eee eee ey 


and this has to be compared with the corresponding result of Bohm and 


S2 
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Pines given in formula (11) and table 1. Both results have been plotted 
against the interelectronic spacing r, in fig. 3 and it can be seen that they 
nearly coincide in the range of actual metallic densities. 

To complete the comparison in the case where for a fixed x; just one 
of the differences | y,—x,| becomes small while all the others are large 
we have plotted the complete function (22) f(y,—x,)=/(7) in figs. 1 and 2 
for the cases 7,=2 and7,—4. These figures also contain the function (10) 
f(7,;) derived from Bohm and Pines’ work. In the range from r=0 to 
r=2 (r=r,,) both theories give nearly the same values for the modulating 
function, while for larger values of ry the Wigner treatment gives higher 
values. 

As we have already mentioned above, the Wigner treatment breaks 
down when for one (or several) x; several of the y,—x, become simul- 
taneously very small and therefore the h,, (y;—x,) large, ie. when 


> hav —%) Jee ee ee 
j= 


becomes large. But this difficulty can be overcome to a certain extent. 
From the case looked at above in which just one of the h,, became large 
while all the others were negligibly small, it can be seen that the Wigner 
treatment represents a good approximation when 


Me 
= 


col ¥j =X 0-D ee eee ne 


The numerical value depends, of course, on the degree of approxima- 
tion desired. We have therefore to deal with the case when this inequality 
is no longer fulfilled, i.e. when one electron comes close together with m 
electrons, m>1, of opposite spin. We have to find the modulating factor 
brought about by this configuration. Obviously there are two possible 
cases. Hither the single electron has an x-coordinate and the m electrons 
have y-coordinates or vice versa. According to (21) the corresponding 
factor in the modulating function in the first case is given by 


m 
USA Fase) Se ee es 
j= 
and in the second case by 
m 
TT] (Pes (ys —263) ee en OD 


We have assumed that (25) is not fulfilled so that (26) no longer represents 
a good approximation. But in the range of actual metallic densities h,.. 
is always <0-5. Therefore (27) still represents a reasonable approxima- 
tion, while (26) is no longer valid. It was just the lack of symmetry in 
the positional coordinates x; and y, of formula (21) which has enabled us 
to avoid the breakdown of the Wigner treatment when (24) becomes large. 
However, when the h,, involved in (26) and (27) are small, both formulae 
give nearly the same numerical value. 
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§7. A ChassicaL ANALOGUE 


Although it does not mean very much we would like to mention here a 
relation of the modulating function with a classically defined quantity. 
In the case of two electrons which are the distance r apart, the corre- 
sponding term in the modulating function for r—0 is given by 


T1O=8 5 
in the Wigner treatment [see formula (23)| and by 
f(0)=exp (—0-352 . r,1/2 . B) 
in the work of Bohm and Pines [see formula (11)].. The quantity /(0) is 


related to a classically defined quantity in the following way. The mean 
kinetic energy of an electron ina Fermi gas is given in atomic units by 


—9).9] /p 2 

€kin— - 21/r, 5 
If we look at an electron of this kinetic energy as a classical particle and 
let it run against an electron at rest, it comes to rest at a distance 7, 


(measured in multiples of the Bohr radius a)) where the whole of its kinetic 
energy has been transformed into potential energy. Then the relation 


Pelt ==T | ook 


holds. If we now plot 7y/7, against (0) for different values of the inter- 
electronic spacing we see that a linear relation exists between both 
quantities. The result has been illustrated in fig. 6. 


Fig. 6 
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ABSTRACT 


An analysis is made of liquid immiscibility in metal systems. It is 
shown that the simple Hildebrand rule, which holds for non-polar liquids, 
can be modified to apply to metal systems, by the introduction of a new 
term which allows for the increased energy due to a difference in electro- 
negativity of the components. The electronegativity term involves the 
maximum number of bonds which can form between two given metals, 
and is considered to be a function of their relative atomic sizes. On the 
theory suggested, out of a total of 529 alloy systems considered, 426 can 
be accounted for compared with 312 on the simple Hildebrand rule. Of 
the exceptions, several can be explained on the basis of errors in the caleu- 
lation of the energy factors involved, and others can be brought into 
agreement on the assumption that the electronegativities of some specified 
elements vary with the elements with which they are alloyed. The 
application of the modified hypothesis to ternary systems is discussed. 


$1. InTRoDUCTION 


THE factors affecting the alloying behaviour of the commoner metals 
and metalloids have been considered in great detail in the last two decades, 
and for many binary metallic systems, it is possible to make reasonable 
estimates of the extent of the solid solubilities and likelihood of compound 
formation. Much of the pioneer work in this field has been carried out by 
Hume-Rothery, Raynor and their collaborators and their conclusions 
have been recently summarized (Hume-Rothery and Raynor 1954). The 
principal factors which govern the extent of the solid solubility of one 
metal in another are the difference in their atomic diameters, their relative 
valencies and the difference in their electro-chemical nature. If the 
difference in atomic size between two elements is greater than about 14%, 
then the strain in the lattice produced by substitution of one type of atom 
for the other becomes so large that the solubility range is restricted. The 
extent of the solubility also depends on the lattice structures of the two 
elements, since a marked difference in the type of bond formed in the pure 
metals can lead to restricted solubilities ; other conditions being satisfied, 
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complete mutual solubility can only occur when the components have 
similar crystal structures. Hume-Rothery also postulated that other 
conditions being equal, a metal of lower valency is more likely to dissolve 
one of higher valency than vice versa. This generalization is not always 
followed, especially in the case of polyvalent elements, and attention has 
to be paid to the Brillouin-zone characteristics of the solvent element. 
If there is a large difference in electrochemical nature between the two 
elements, then a strong tendency to form stable compounds can be 
expected and this will decrease the extent of solid solubility, even if the 
size factor is favourable. Compounds or stable intermediate phases are 
also formed in systems for which the relative atomic sizes are such that 
a certain type of packing is easily achieved ; the structure of the new 
phase may be either completely different from those of the parent lattices 
or based on them as in a superlattice. 

In general, it can be assumed that for a system showing marked solubility 
in the solid state, there is no a priori reason for restricted solubility in the 
liquid condition, so that if the characteristics of two elements are such that 
solid solution is predicted, then no liquid immiscibility would be expected. 
There are many binary systems which show restricted solid solution but 
are miscible in all proportions in the liquid phase. This is due to two 
main differences between the liquid and solid states. The distribution of 
atoms in a liquid is less restricted than in the solid so that the liquid can 
accommodate atoms of different sizes more easily and hence the range of 
relative atomic size for complete liquid miscibility is much wider than 
that for solid solution. Secondly, although a high electrochemical factor 
can restrict solid solution, it may help to attract two elements together in 
the liquid state to give liquid miscibility. 

The object of the present work was (a) to establish the requirements for 
complete liquid miscibility of two elements and (b) to suggest a working 
basis for the prediction of the effect of a third element on the range of 
immiscibility of two given elements. 


§2. Review oF Previous Work on Miscipiuiry or Liguip MEtTats 


Only two serious attempts to obtain a general theory of liquid metal 
immiscibility have been reported. Axon (1948) based his work on the 
relation between the size factor and the difference between the melting 
points of the two metals, while Hildebrand and Scott (1950) applied the 
relationships for immiscibility in non-polar liquids to liquid metal systems. 

Axon attempted to systematize the knowledge of the factors influencing 
the form of binary equilibrium diagrams, and considered systems in which 
no intermediate phases are found, neglecting those that give complete 
solid solubility. The systems were of two types, viz.: those showing 
liquid immiscibility and a eutectic respectively. A diagram was con- 
structed in which each system was represented by a point having co- 
ordinates representing a size factor and a temperature factor. These 
were defined as follows. 
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} ei 
Size factor— —— x 100, 
where L and 8 are the closest distances of approach of the larger and 
smaller atoms respectively in the pure elements ) 
Temperature factor = aan | 
(9,465) 7 
where 6, and 6, are the melting points of the two metals on the absolute 
scale and Jé is the difference between them. 
The diagram obtained could be divided into three regions as illustrated 
in fig. 1. Region A was limited to a size factor of 30 % and a temperature 
factor of 0-175, and only included systems of a simple eutectic type in 
which the eutectic composition did not vary much from 50 atm. %. For 
region B, representing systems having a size factor greater than 30%, all 
but one gave a miscibility gap. Systems in region C were approximately 
equally divided between the two types. For region C it was noted that 
the eutectic composition moved progressively away from 50 atm. °% as the 
temperature factor increased beyond 0-175. 


Fig. 1 
B 
oc 
2) 
om 
O 
= 
+ 30 
N } 
wn 
C 
0.175 
TEMPERATURE FACTOR 


Axon’s diagram for immiscible and eutectic-type metallic systems, 


This approach is of interest but has been limited to about 80 binary 
systems of a relatively simple type, and does not include those in which 
liquid immiscibility has been reported over a limited region of composition. 
Further, no account is taken of the electrochemical nature of the elements 
forming the alloy and no suggestions were made of extending the analysis 
to ternary systems. It could be surmized from Axon’s hypothesis that 
the range of miscibility of two metals could be increased by the addition of 
a third element for which both the atomic size and the melting pomt were 
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intermediate between the corresponding properties for the other two. 
The validity of this has not been examined. 

The other general approach to the problem was that of Bildebeata and 
Scott (1950) who calculated that the excess free energy of formation of a 
liquid solution was very similar to the theoretical expression for the energy 
of mixing of a regular solution. The energy of mixing 4H ™ is given by 


ARM = v[ (S£)"— (S2)" al babe 


where V is the average atomic volume of the solution, V, and Vz the 
atomic volumes of components A and B, 4H," and JE," their heats of 
vapourization, and 4, and ¢, the volume fractions of the two components. 
Hildebrand postulated that when the heat of mixing is sufficiently great, 
separation into two liquid phases will occur and the condition for complete 
miscibility was deduced as 


AE 1/2 AH.” \¥212 
V+ Vo)| ( aa —( 7) | <2RT. 


The term (4£,"/V 4)! was considered to be a measure of the binding 
energy of component A and was called the * solubility parameter ’ 6. 

The condition for immiscibility therefore becomes 

3(VatV 5g) (84-53)? >2RT. 

This relationship proved to be a most useful guide to the behaviour of 
non-polar liquids. When applied to liquid metal systems, Hildebrand 
and Scott (1950) found that of 47 immiscible binary systems, 40 obeyed 
the relationship, 2 were borderline cases and 5 showed marked disagree- 
ment. On the other hand, nearly 50°, of the alloys which are known to 
be miscible in the liquid state, would be immiscible according to Hilde- 
brand’s rule. These include many technologically important systems 
such as Cu-Sn, Cu-Zn, Mg—Al, Al-Zn, AI-Ni, ete. for which the nature 
of the equilibrium diagrams is in little doubt. Many of these exceptions 
to the immiscibility rule show compound formation, suggesting that there 
may be a relatively high electrochemical attraction between the elements 
which is not taken into account by the Hildebrand relationship. If a 
factor could be introduced into the expression to allow for the electro- 
chemical factor, the exceptions to the simple rule might be removed. 
This possibility is discussed in the following sections. 


§ 3. HLECTRONEGATIVITY CONCEPTS 


There are numerous possible ways of defining the chemical affinity 
between two elements ; the main difficulty is to find which is appropriate 
to the atoms in the liquid metal state. One simple method of classifying 
the elements is from their relative position in a table of electrode potentials 
but, since the potentials refer to the ions in aqueous solution, it is unlikely 
that these apply to liquid metals. In fact a study of various binary 
systems quickly revealed that the difference in electrode potential between 
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two metals could not. be simply correlated with the behaviour of their 
molten alloys. Other methods of expressing electronegativities were 
therefore considered and included the following. 

(a) Mulliken (1934 a, 1934b, 1935, 1949) defined the electronegativity 
(X) for an univalent atom as the mean of the ionization potential J and 
the electron affinity Z. 

ie. X=}([+ £) electron volts. 

There is a lack of knowledge of the values of # for various elements 
and Mulliken’s electronegativity concept has received little application. 

(6) From consideration of the heats of formation of halides of various 
elements, Pauling (1932) developed an electronegativity or X scale such 
that the square of the difference in the numbers assigned to any two 
~ elements (X ,—X,)* is roughly proportional to the heat 4 accompanying 
the formation of a bond between them. If D(AB) represents the dissocia- 
tion energy of a bond between A and B, then according to Pauling : 

D(AB)=3[D(AA)—D(BB)]+ 4. 
4 is found empirically to be additive, and X, and X, are chosen so 
that their difference is proportional to 1/ J. 

(c) Gordy (1946) examined the electronegativity of some thirty elements 
as given by Pauling, and found that X could be expressed in terms of the 
number of valency electrons n’, and the single covalent radius of the 
atom 7 according to the expression 

DAES = (n’+1)+0-50. 

This expression can be derived by defining the electronegativity of a 
neutral atom as the potential Ze/r, where Ze is the effective nuclear charge 
of the atom acting on a valency electron at a distance r from the nucleus 
equal to the single-bond covalent radius of the atom. Assuming that the 
electrons in the closed shells exert their full screening power and those in 
the valency shell have a screening constant of 0-5, Z is found to be equal 
to 0-5 (n’+1). On this basis, the electronegativity can be expressed as 
[K, (n’+1)/r+K,], where both K, and Ky, are constants, the electro- 
negativity scale being arbitrary. The most marked deviations from Gordy 8 
rule were the values for copper, silver and gold when »’ was taken as unity, 
although if the Pauling valencies for these elements are substituted in the 
Gordy equation, the X values agree closely with those estimated by 
Pauling. ; 

When the electronegativities are recalculated from Gordy’s equation 
using the Pauling valencies throughout, there are many marked differences 
with the experimental values as given by Pauling. There is no doubt that 
the electronegativity of a given element depends on the valency it assumes 
in a particular system, but there are no obvious grounds for considering 
the Pauling valency for some elements and not for others. 

Altshuller (1954) has considered the electronegativity values for oe 
copper, zine and gallium series of sub-group elements as calculated by 
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Gordy. He suggests that the values for copper, silver and gold are ‘too 
high, and those for zinc, mercury and cadmium are too low. Recalcula- 
tion of the electronegativities from the dissociation energies of the homo- 
polar halides and hydrides give the following values compared with those 
quoted by Gordy. A} 


Hlement Cu Au | Zn | Cd | He | Gao} tn ) at . 


Altshuller 
Gordy 


(d) Walsh (1948) found that the force constant controlling the vibration 
of a diatomic hydride is roughly proportional to the square of the first 
ionization potential of the element, and suggested that the strength of the . 
bond between hydrogen and a given element should be taken as a measure 
of electronegativity. In general, there is a reasonable relationship 
between the values given by Walsh and Pauling but there are differences 
which are most marked for elements such as Be, Mg, Hg, Si, ete. Walsh 
associated these exceptions with a ns? pair of electrons in the outer sheli 
of the normal state of the atom, and postulated that the bond strength 
with hydrogen for these elements should be related, not to the ionization 
energy of the ground state, but to the lower energy required to ionize an 
electron from an excited valency state. Some support for this is obtained 
from a plot of the force constants of diatomic A—O molecules against the 
ionization potential of A. Instead of a single curve as for the A—H bonds, 
the points for each group of the Periodic Table lie on a characteristic 
curve. Therefore, although the hydride bond strength is a useful 
guide to electronegativity in some cases, it does not provide a general 
basis for comparison of various elements. 

(e) Haissinsky (1946) recognized that the electronegativity of a given 
element will vary to some extent with the element with which it is com- 
bined ; this can be appreciated from the difference in the electronegativity 
values derived from the oxides and bromides as compared with values 
obtained from the fluorides and chlorides. There will also be a tendency 
for the electronegativity to increase with the valency and Haissinsky 
deduced that this would be less marked for the metalloid elements and 
noble metals, From a consideration of the heats of formation of the 
chlorides, bromides and iodides, Haissinsky developed a series of electro- 
negativities for a number of elements, which was extended by a process of 
extrapolation and interpolation from group to group and period to period 
of the Periodic Table. The electronegativity values are generally 
in good agreement with those of Pauling, the most notable exception 
being gold and mercury. 
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§4. CHotce or ELEcTRONEGATIVITY VALUES 


The chosen values of electronegativity for 56 elements are given in 
table 1 together with the source and the valency state where appropriats: 
For many elements, the values given by Pauling and Haissinsky are in 
reasonable agreement and have therefore been adopted for the present 
work. With few exceptions, the Haissinsky values have only been used 
where a Pauling figure is not available. In cases in which Haissinskv 
differentiates the electronegativity for various valency states of a given 
element, the value selected is that giving the best agreement with the 
data available, e.g., Mn™, MoY!, TH! VY and W!Y; in the case of Mn, 
Mo and V, this corresponds to the most probable valency in alloys 
but the TI"? and WY! states are more likely than TI! and W!Y. The 
choice of electronegativity value is discussed further in later sections. 


§ 5. PROCEDURE 


Data on alloy systems have been taken from various sources including 
Smithell’s Metals Reference Book, A.S.M. Metals Handbook, Metal 
Abstracts and the books by Hansen (1936) and Janecke (1940) on alloy 
constitution. The systems have been divided into two main categories, 
viz. elements miscible in aii concentrations in the liquid state and those 
in which some liquid immiscibility has been reported, whether extensive 
or limited to a narrow range of composition. For some immiscible 
systems, there is little doubt as to their authenticity but for others, either 
the data is contradictory or the work was carried out many years ago 
with relatively impure materials and inadequate techniques. Systems 
which may have been reported incorrectly as immiscible are therefore 
listed separately. For ease of classification, the systems have been 
tabulated according to their group in the Periodic Table. 

The Hildebrand factor }(V4+V,)(64—5,)? was estimated for each 
binary system from the solubility parameters of Hildebrand and Scott 
_ where available (see table 1) otherwise the parameter of a given element 
was calculated from the most probable heat of vapourization. The 
solubility parameter of the element of lower melting point was extra- 
polated to the temperature at which the second element melts, assuming 
the variation of the solubility parameter with temperature to be linear and 
approximately the same for elements in a given group. As a guide to the 
temperature variation of the solubility parameters of elements in a given 
group, the data given by Hildebrand and Scott (1950) were plotted and 
extended. ioe ot) 

Although at normal pressure it is difficult to assign liquid immiscibility 
to a system in which one metal boils before the other melts, such systems 
have been included in the analysis, since increasing the pressure 
to sufficiently raise the boiling point of one element should have little 
effect on its solubility parameter as defined by Hildebrand. For con- 
venience, the difference between the Hildebrand factor and the energy 
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Table 2. Data for Systems with Authentic Immiscibility 


Hildebrand | Electronegativity 


System excess factor k 
(cal/g mol) (electron volts) 
Li-Na 6,300 0-01 27:3 
Li-K 24,500 0-04 26:6 
Li-Rb 35,800 0-04. 38:8 
Li-Cs 49 000 0-09 23-6 
Li-Fe 34,200 0-49 3:0 
Na-Mg 1,100 0-09 0-53 
Na-Ca —(1,500) 0-01 a 
Na-Zn 6,700 0-09 3-2 
Na-—Cd 300 0-04 0-33 
Na—Al 42,000 0:36 5-1 
Na-La 19,800 0-09 9-2 
Na-—Fe 96.400 0-64 6-5 
K-Mg 17,900 0-16 4-9 
K-Zn 30,600 0-16 8-3 
K-Cd 11,700 0-09 5:6 
K-Al 107,800 0-49 9-54 
K-—Pb 25,800 0-49 2-3 
K-Fe 215,200 0-81 11-5 
Rb-Fe 272,000 0-81 14-6 
Cs—Fe 320,000 1-00 13-9 
Cu-Tl 20,200 0-09 14-1 
Cu=Pb 28,000 0-49 2-5 
Cu-Cr —(8,800) 0-49 — 
Cu-U —(5,100) 0-81 — 
Cu-—Se 50,800 0-04 55:1 
Cu—Te 65,600 0-01 284-4. 
Cu-Co —(2,300) 0-25 = 
Ag-Ta 37,700 0-25 6-5 
Ag-V 7,100 0-01 30:8 
Ag—Cr —(3,300) 0-16 -— 
Ag—-W 66,400 0-09 32:0 
Ag-U —(1,700) 0-36 = 
Ag-Se 24.700 0-25 4-3 
Ag—Fe 4,700 0-04 5-1 
Ag-Co 15,100 0-04 16-4 
Ag-Ni. 9,600 0-04 10-4 
Ag-Ir 36,700 0-04 39-8 
Au—Mo 7,600 1-00 0:33 
Au-W 21,400 2-25 0-41 
Au-Se 41,300 0-49 Ba 
Be-Mg 72,400 0-09 349-0 
Be-Sn 34,500 0-04. 374-0 
Mg—Fe 41,600 0-25 72 
Mg-W 98,800 0-16 26'8 
Ca-Cd | —(4,500) 0-01 — 
Ca-W 163,400 0-36 19-7 
Ca—Fe 74,400 0-49 6-6 
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Table 2—(continued) 


Hildebrand | Electronegativity 
System excess factor k 
(cal/g mol) (electron volts) 
Ba-Fe 131,000 0-64 8-9 
Sr—Fe 124,000 0-49 11-0 
Zn-Tl — (2,300) 0-49 _ 
' Zn—Pb —(2,600) 0-09 — 
Zn—Bi —(2,000) 0-36 — 
Zn—W 105,400 0-16 28-6 
Cd-Al 15,600 0-16 4-2 
Cd-Ga 8,800 0-09 4-2 
Cd-Cr 28,000 0-16 7-6 
Cd-Fe 48,100 0-36 5:8 
Hg-Ga 23,400 0-16 6-3 
Hg-Si 29,100 0-64 1-97 
Hg—Nb 123,000 0-64 8:33 
Al-In 5,700 0-01 24-7 
ALT] 12,100 0-16 3:3 
Al-Pb 10,900 0-00 oe) 
Al-Bi 15,500 0-09 7:5 
La-Ti 18,000 0-16 4-9 
La—Mn 15,500 0-04 16-8 
Ga-Tl 6,200 0-25 Re 
Ga-Pb 4.900 0-01 21-2 
Ga—Bi 8,200 0-16 2-2 
Ga-Te 28,900 0-49 2-6 
In-Te 14,700 0-49 1:3 
TLSi 10.900 0-01 43-4 
Tl-Ge 7,300 0-04 (os) 
TIL-As 2.300 0-O1 10-0 
Tl-Se 3,900 0-25 0-68 
Tl-Te 4.700 0-04 5-1 
TI-Mn 22,700 0-25 3-9 
Tl-Fe 41,300 0-04 44-8 
TI-Co 62,400 0-04 67-6 
TI-Ni 50,300 0-04 54-5 
Th-U — (4,500) 0-04 = 
Si-Pb 10,500 0-09 5-1 
Si-Bi 15,700 0-00 00 
Sn—Cr 3,300 0-04 3-6 
Sn—W 72,800 0-01 3156-0 
Sn—Fe 18,600 0-00 oO 
Pb-Zr 21,500 0-01 - 93-2 
Pb-Cr 19,500 0-00 oe) 
Pb-Se 5,300 0-81 0-28 
Pb—-Mn 20,200 0-01 87:6 
Pb-Fe 40,100 0-04 43-5 
Pb-Co 63,300 0-04 68-6 
Pb-Ni 30,900 0-04 55°2 
Bi-Cr 26,900 0-09 13-0 
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~ Table 2—(continued) 


Hildebrand Electronegativity 
System excess factor 
(cal/g mol) (electron volts) 


Bi-Se 1,500 0-36 


Bi-Mn 27,300 0-16 
Bi-Fe 51,000 0-01 
Bi-Co 78,200 0-01 
U-Ce 8,500 0-04 


Table 3. Data for Miscible Systems 
ges So EE ee a a ee 


Alloy Hildebrand | Electronegativity 


system e€Xcess factor k 
(cal/g mol) (electron volts) 
Li-Cu 24,100 1:44 0-73 
Li-Ag 30,500 0:81 16 
Li-Hg 4,300 0-00 00 
Li-Al 8,100 0:25 ihe 
Na-K 2,700 0-01 11-7 
Na-Rb 6,300 0-01 27-3 
Na-Cs 12,100 0-04 1321 
Na-Au 60,100 4:84 0-52 
Na-Ga 28,400 0:25 4-9 
Na-Tl 3,700 1-00 0-16 
Na-—Ge 28,100 0-64 1-9 
Na-Sn 19,800 0-64 1:3 
Na-—Pb 6,000 0-36 0:72 
Na-Sb 12,900 1-44 0-39 
Na-Bi 4,300 0-81 0:23 
K-TI 22,400 1-21 0-80 
K-Sn 58,000 0-81 Sol 
K-Sb 42,400 1-69 1-1 
K-Bi 22,100 1-00 0-96 
Rb-Hg 3,800 0-04 4-1 
Cs-Hg 8,400 0-09 4-0 
Cu-Mg 31,500 1-00 1-37 
Cu-Ca 60,600 1-44 1-8 
Cu-Zn 11,200 1-00 0-49 
Cu—Cd 34,600 1-21 1-2 
Cu-Hg 54,800 1-44 [+7 
Cu—La 17,100 1:00 0-74 
Cu-—Ca 4,300 0-64 0-29 
Cu=in 20,000 0-64 1-4 
Cu-Ge 4,100 0-25 0-71 
Cu-Sn 10,600 0-25 1:8 
Cu-As 8,600 0-04 9-3 
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Table 3—(continued) 


Hildebrand | Electronegativity 
@XCESS factor k 
(cal/g mol) (electron volts) 


Alloy 
system 


Cu-Sb 17,000 0-01 
Cu-Bi 34,800 0-16 
Cu—Ce ~ 13,600 ico 
Cu-—Pr 17,300 1-21 
Ag-—Be 13,500 0-16 
Ag-Mg 10,000 0-49 
Ag—Ca 22,800 0-81 
Ag-Sr 47 ,000 0:81 
Ag-Zn 3,900 0-49 
Ag—Cd 0-64 
Ag-Hg ; 0-81 
Ag-—In 0-25 
Ag-Tl 0-00 
Ag-—Pb 0-16 
Ag-Sb 0-04 
Ag-Bi 0-01 
Ag-Pt 0-04 
Au-Be 
Au-Mg 
Au-Ca 
Au-Zn 
Au-Cd 
Au-Hg 
Au-La 
Au-In 

Au-Tl 

Au-Sn 
Au-—Pb 
Au-As 
Au-Sb 
Au-Bi 

Au-Te 
Au—Co 
Au-Ni 
Au-Rh 
Au-—Pt 
Au-—Ce 
Au-—Pr 
Be-Al 
Be-Si 

Mg-Sr 

Mg-—Ba 
Mg-Hg 
Mg—Al 
Mg-Ga 
Mg-Si 

Mg-Ge 
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Table 3—(continued) 


Hildebrand Electronegativity 
excess factor 
(cal/g mol) (electron volts) 


24,200 0-04 
62,200 0-25 
24,300 0:25 
33,900 0-64 
6,700 0-49 
49,700 0-25 
22,700 0-49 
5,200 0-25 
59,600 0-36 
28,800 0-64 
8,300 0-36 
6,600 0-04 
4,400 0-09 
7,400 0-16 
14,100 0-16 
6,900 0-25 
12,800 
8,900 
10,200 
23,600 
34,000 
31,900 
39,500 
300 
1,300 
4,700 
3,900 
62,600 
77,800 
31,100 
6,200 
18,500 
8,000 
6,600 
68,500 
24,600 
1,500 
600 
4,200 
5,100 
25,700 
45,000 
28,000 
38,200 
2,200 
9,100 
7,700 
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Table 3—(continued) 


Hildebrand | Electronegativity 


Alloy excess factor k 
system (cal/g mol) (electron volts) 
Al-Pt 11,300 0-36 ee 
Al-Ce 1,800 0-16 0-49 
Al-Pr 400 0-16 0-11 
La-Ni 37,000 0-25 6-4 
Ga-Ni 12,600 0-09 6-1 
In-Zr 2,500 0-04 2-7 
TI-Pt 64,100 0-04 69-5 
Ti-Nb 3,700 0-04 4-0 
Ti-Ta 20,200 0-04 21-9 
Ti-Mo 8,000 0:25 1-4 
Ti-W 22,200 0-00 oe) 
Ti-Co 900 0-01 3°9 
Zr-Sn 2,700 0-01 11:7 
Zr-Nb 4,100 0-04 4-4 
Zr—Mo 9,000 0-25 1-6 
Zr-Sb 17,400 0:25 3-0 
Zr—Ru 19,600 0-25 3°4 
Th-V 2,600 0-49 0-23 
Th-Nb 27,700 0-49 25 
Th—-Bi 27,300 0-49 2°4 
Th-Se 36,200 1-69 0-93 
Th—Te 43-800 1-00 1-9 
Th-Fe 2,300 0-36 0-28 
Th—Co 10,600 0-36 1:3 
Th-Ni 8,700 0:36 1-0 
Si-Ta 29,500 0-16 8-0 
Si-Sb 4,900 0-09 2-4 
Si-Mo 16,000 0-09 7:7 
Si-Co 7,700 0-01 33-4 
Si-Ni 5,000 0-01 21-7 
Si-Ce 900 0:49 0:08 
Ge-Bi 3,200 0-01 13-9 
Ge-Te 28,600 0-16 7°38 
Ge-Fe 3,600 0-00 oo) 
Ge-Co 9,600 0-00 eo) 
Ge-Ni 7,600 0-00 CO 
Sn—Pb 1,000 0-04 Veal 
Sn-Bi 3,200 0-01 13-9 
Sn—U 7,400 0-16 2-0 
Sn—-Te 19,900 0:16 5:4 
Sn—Mn 4,700 0-09 2-3 
Sn—Co 31,700 0-00 oe) 
Sn-Ni 29,100 0-00 CO 
Sn-Pt 33,100 0-16 9-0 
Pb-Te 5,700 0-36 0:69 
Pb-—Pd 24,100 0-25 4-2 
Ph—Pt 66,800 0:36 8:0 
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- Table 3—(continued) 
eee a ee 


Alloy Hildebrand Electronegativity 


excess factor k 
com (cal/g mol) (electron volts) 
V—W 2,700 0-04 2.9 
Nb-U 3,400 0-25 0-59 
Ta-U 12,500 0-01 54-2 
As—Mn 2,600 0:36 0-31 
As—Fe 14,200 0-09 6-8 
As-—Co 12,200 0-09 5:0 
As-Ni 19,900 0-09 9-6 
As-Pt 29,100 0-01 126-0 
Sb-Cr 18,300 0-36 2-2 
Sb-Se 10,000 0-09 4-8 
Sb-Te 10,700 0-00 fo) 
Sb—Mn 13,500 0-49 15 
Sb—Fe 32,600 0-16 8-8 
Sb-Co 48,600 0-16 13-2 
Sb-Ni 45,600 0-16 12-4 
Sb-—Pd 16,700 0-01 79.4 
Sb—Pt 57,000 0-00 oe) 
Bi-U 28,800 0-25 5-0 
Bi-Ni 65,600 0-01 284-0 
Bi-Rh 81,000 0-09 39-0 
Bi-Pd 21,900 0-04 as et 
Cr—-W 6,700 0-01 29-1 
Cr—Re 7,300 0-49 0-65 
Cr-Ir 8,800 0:36 0:82 
Mo-U 8,000 0-64 0-54 
W-U 23,100 0-09 11-1 
W-Pd 5,100 0-16 1-4 
Se—Fe 57,400 0-49 Hel 
Fe—Ce 27,600 0-36 3:3 
Co—Ce 38,800 0:36 4-7 
Ni-Ce 35,700 0-36 4-3 
Ni-—Pr 38,900 0:36 4-7 


Table 4. Data for Systems of Doubtful Immiscibility 


Alloy Hildebrand | Electronegativity 


system excess factor k 
(cal/g mol) (electron volts) 

Na-—Ag 35,900 1-00 1-6 

Cu-V —.(6,900) 0:16 -— 

Cu—Mo —(800) 0-01 — 

Cu-W 3,200 0-36 0-39 

Cu-Fe — (6,200) 0°25 — 


Ag-Ba 52,000 1-00 22:5 
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Table 4—(continued) 


Hildebrand | Electronegativity 
eXCeSs factor 
(cal/g mol) (electron volts) 


Alloy 
system 


Ag-Ti —(5,900) 0-09 
Ag—Mo 39,800 0:04 
Ag-Te 37,300 0-04 
Ag—Mn —(3,200) 0:25 
Be-U 3,300 0-04 
Mg—La 400 0-00 
Mg-U 29,200 0-01 
Mg-Ni 52,600 0-25 
Zn—Se 5,200 1-44 
Cd-Si 11,900 ; 
CdSe 100 

In-Si 3,900 

Sn—Se 13,800 

Pb—-W 131,000 

Pb-U 25,800 

As-Te 21,400 
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term 2R7' (where T° absolute is the melting point of the more refractory 
component) for each system will be termed the ‘excess Hildebrand 
factor’ and this quantity has been listed for immiscible, miscible and 
doubtful systems in tables 2, 3 and 4. Miscible systems for which the 
Hildebrand rule holds, i.e., where excess Hildebrand factor is considered 
to be negative, have not been included in the tables. There are a few 
immiscible metal systems in which the excess Hildebrand factor is nega- 
tive ; these have been listed for completeness in table 2. 


§6. RESULTS oF PLOTTING THE Excess HILDEBRAND FacTOR AGAINST 
ELECTRONEGATIVITY 


If it is assumed that the failure of the Hildebrand rule when applied to 
metal systems is due to an electronegativity factor, then immiscibility 
would only be expected for a high excess Hildebrand factor with a low 
electronegativity. Conversely, if the electronegativity were high, then 
in general, this would counteract a high excess Hildebrand factor, and 
miscibility would be obtained.. To test this hypothesis, the electro- 
negativity against excess Hildebrand factor has been plotted in fig. 2 for 
low values of these parameters and in fig. 3 for the remaining ranges. 

Undoubtedly, although there are several exceptions, the dots repre- 
senting immiscible systems tend to lie above the crosses representing 
miscible alloys. An arbitrary line has been drawn through the origin in 
figs. 2 and 3, which divides each diagram so that there is a greater number 
of dots in the top left-hand portion and a greater number of crosses in the 
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oles right-hand corner, The slope of the dividing line in each diagram 
ee i el Over a fairly wide range without markedly affecting the relative 
stribution of miscible and immiscible systems. Accepting that the 


assumptions are correct, it is profitable t : vod 
slope of this line. P e to examine the significance of the 
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Plot of excess Hildebrand factor against electronegativity factor (low values). 
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§7. SuGGESTED MopiFIcaTion TO HILDEBRAND’S HYPOTHESIS 
According to Hildebrand, immiscibility in non-polar liquids should 
occur if 
3(VatVz) (84-52)? >2RT. 
On the present theory, for metallic systems in which there is an electro- 
negativity factor, it is necessary to introduce a new term Hz, which repre- 
sents the binding energy, so that the condition for immiscibility becomes 


4(V +V 5) (84-82)? Lp >2 KP. 
The binding energy, Hz which is required for the above expression is the 


combined energy of all the bonds which can form between the two elements. 
A simple picture of a bond is not possible as resonance always occurs to 
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some extent between the various types which include ionic, covalent and 
metallic. If however the mean fraction of each type of bond is repre- 
sented by f;, fg and fy respectively and H;, Hg and Hy are the energies 
associated with one bond if wholly of one type, then the binding energy 
for n bonds is given by 
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Plot of excess Hildebrand factor against electronegativity factor (high values). 


At the moment, there appears to be no accurate method of assessing 
the energy #y, associated with a metallic bond but it is generally con- 
sidered that it is more closely related to the covalent than the ionic bond. 
If, as a first approximation, we can assume H,=E,, then the expression 


0.4 0.6 0:8 1.0 1.2 


for the binding energy becomes 


By=n[Eo(fotfu) tH. 


1.4 


1.6 
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By definition, f;+-fy-+-fo=1, so that the above expression can be 
rewritten 


Lz=n[Eo(1—f,)+fE 7. 


If this represents the covalent-ionic resonance energy of n bonds then 
by definition (see Pauling), it is approximately equal to (X ,—X,)?; 
so that : 

E,~n(X4—X pz)? electron volts 


or 23 060 n (X,—X >)? calories/g. 
Thus to a first approximation, the condition for immiscibility becomes 


2(V4tVp)(84—83)?—23 060 n(X ,—X p)*>2 RT. 
Rearranging this expression, immiscibility may be expected if 


3(VatVp)(84—8p)?—2 RT 
23 060(X,—Xp) 

The above relationship states that if two elements can form a sufficient 
number of A—B type bonds, then the energy associated with bonding will 
counter-balance the binding energy associated with the elements as 
expressed by Hildebrand’s solubility parameter, and liquid miscibility 
will occur. If two elements cannot form at least k A—B type bonds, 
where & is given by the expression 

2(VatV2z)(4—8z)"—2 RT 
23 060(X,—Xp)? 


then the two liquids will be immiscible. 

In addition to assigning a precise energy to a given bond, there is also 
the problem of deciding the maximum number of bonds per atom which 
two metals can form: i.e. the maximum value of n. Presumably the 
maximum number of bonds formed by a given metal is directly related to 
the number of electrons available for bonding, i.e. to its maximum valency. 
The valencies of the elements considered vary from 1 to 6 so that, as a 
first approximation, we can assume that if the ratio of the excess Hilde- 
brand factor to the electronegativity factor is greater than 6, then 
immiscibility could be expected. On the other hand, if this ratio is less 
than 1, the metals should be miscible. (The arbitrary dividing lines in 
figs. 2 and 3 correspond to a k of about 5 which is reasonable on the above 
basis.) Ifthe ratio is between 1 and 6, the type of system obtained may 
be affected not only by the valencies of the two elements, but also by their 
relative atomic sizes which will determine the possible types of packing 
and hence the coordination number. The upper limit for the average 
number of nearest neighbours for both solids and liquids is 12 (although 
this number may be exceeded locally) so that on the average, the maximum 
number of near neighbours of unlike atoms is 6, which appears to be an 
additional reason for choosing this as the maximum number of bonds. 

The possible types of packing and hence the maximum number of 
ponds that can form in a given system will depend on the relative sizes of 


>n. 
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the atoms or ions. To see if any general guide could be obtained on the 
effect of relative atomic size on the critical k value (which will determine 
whether immiscibility occurs or not), the calculated k values in the range 
1 to 6 for both miscible and immiscible systems were plotted against the 
percentage difference in the Goldschmidt atomic diameter (based on the 
smaller atom). As for figs. 2 and 3, miscible systems were represented 
by crosses and immiscible systems by dots, and the plot is shown in fig. 4. 
A curve has been drawn which encloses the maximum ratio of dots to 


Fig. 4 
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Plot of k value against difference in atomic size. 


crosses, 1.e. of immiscible to miscible systems. Only 7 of the 24 immiscible 
systems are outside and 16 out of 80 miscible systems are inside the area 
enclosed by the curve. Little reliance can be placed on the shape of 
the dividing curve since errors in both the excess Hildebrand and electro- 
negativity factors can appreciably affect the calculated value of k. In 
addition, the percentage difference in atomic size has been calculated from 
the Goldschmidt diameters, which are estimated from the distance apart 
of the atoms in the pure crystalline solid. At the best, therefore, this 
size difference will only approximate to that appertaining to alloys in the 
liquid state, although the error in calculating the atomic ratio is not likely 
to be significant compared with other errors. The type of packing may 
also be affected by the relative compressibility of the atoms. 

Two conclusions may be drawn from fig. 4. The critical k value 
appears to reach a minimum of two for a percentage difference in size of 
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about 33%, implying that the coordination number is lowest for such 
systems, which is not unreasonable. It is noteworthy that only two 
immiscible systems have a k value between 1 and 2 whereas there are 
38 miscible systems which fall into this range. Figure 4 suggests that the 
critical & value is 6 for systems where the difference in atomic size of the 
two elements is less than 10° but a size difference of from 10 to 60% 
appears to reduce the number of bonds that can be formed. The situation 
is restored again when the atomic sizes differ by more than about 60%, 
which is reasonable as this is the limit at which interstitial packing may 
occur. As discussed by Barrer (1948), the coordination number for unlike 
atoms in interstitial solutions can be six which could permit six bonds. 


§ 8. CONSIDERATION OF EXCEPTIONS 


From the above discussion, a useful working basis for assessing the 
type of liquid equilibrium between two metals is a calculated k value of 
less than 2 for complete miscibility and greater than 6 for some immisci- 
bility. If the & value is between 2 and 6, then the incidence of 
immiscibility also depends on the relative atomic sizes and although no 
rigid relationship can be given at this stage, the curve of fig. 4 can be 
taken as a guide. On this basis, out of a total of 529 alloy systems 
considered, 426 can be accounted for compared with 312 on the simple 
Hildebrand rule. 

Detailed comparison is as follows :— 


Immiscible Miscible 
systems systems 
(i) Total number considered. 99 430 
(ii) Number conforming to simple Hilde- ; ; 
brand rule. 83 224 
(iii) Number conforming to simple modified 
rule. i.e. k>6 for immiscibility, k<2 es 
for miscibility. au 328 


(iv) Number with & from 2 to 6. ; 

(v) Number of type (iv) conforming to size : ‘ 
rule (fig. 4). . 15 2 

(vi) Total conforming to modified rule, i.e. 
sum of types (iii) and (v). 


In addition to the above 529 systems, there were 22 for which the 
experimental evidence was either conflicting or insufficient to draw any 
definite conclusions as to the existence of liquid immiscibility (see table 4). 
The present theory suggests that immiscibility could be ear died in oe 
following 8 systems, Ag—Ba, Ag—Mo, Be-U, Mg-—U, Pb-W , Se As-Te, 
and Mg-Ni. For the remaining 14 systems, the analysis indicates no 
immiscibility, which is in general agreement with more recent investiga- 
tions. 
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A question to be decided is whether the exceptions to the modified rule 
are real or explicable by errors involved in the calculation of k. In view 
of the uncertainties in the experimental determination of the latent heat 
of vapourization of a metal and in the extrapolation of the solubility 
parameter to high temperatures, for many systems the value of (84—53) 
could reasonably be in error by -5. The effect of this on the excess 
Hildebrand factor increases with increase in the difference between the 
solubility parameters and has been calculated for various values of 
(8,—5,). Table 5 gives the maximum and minimum errors for the 
excess Hildebrand factor assuming a combined error of 5 in the solubility 
parameters. 


Table 5. Errors in the Hildebrand Factor for an 
error of 5 in (64,—6,) 


Maximum errors in 
(54 —dz) Hildebrand factor 
(cal/g atom) 


15 +1,750 —1,250 
30 +3,250 —2,750 
50 45,250 —4,750 
65 +.6,750 —6,250 
80 +8250 F750 
100 +10,250 —9,750 


When possible errors of this magnitude are taken into consideration, 
the k value of many miscible systems which appear to be exceptions can 
be brought into agreement although there are still some which cannot be 
explained on this basis. The electronegativity values are generally 
expressed to the nearest 0-1 electron volt so that the possible error in 
(X,—X,) is 40-1. The effect of this error on k decreases with increase 
in the electronegativity factor but can be very important for small 
differences in electronegativity. Further, as previously discussed, the 
electronegativity of a given element depends on the valency it assumes 
and can therefore vary with the nature of the alloying element. For 
example, the electronic behaviour of the transitional elements is dependent 
on their environment and on this basis, their electronegativities cannot be 
represented by single numbers. Other elements whose electronegativities 
may change with environment include mercury, thallium, tin, lead, 
arsenic, antimony, and bismuth; at the moment no satisfactory basis 
exists for assessing the probable variation. An additional complication 
which may arise, especially in the case of the transitional elements, is 
that when the metallic component of the bonding is large, the error 
introduced by assuming that the energy of a metallic bond is similar to 
that of a covalent bond may become significant. This may be a further 
reason for the large number of exceptions among the transitional elements. 
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§9. Exrension or Hypornusis To TERNARY SYSTEMS 


It is often of technological importance to decrease the range of liquid 
immiscibility of a given pair of metals by alloying with a third element. 
There are few rules which enable one to choose the third element to any 
degree of certainty, but it is suggested that the modified Hildebrand 
hypothesis elaborated above can be used to advantage. To extend the 
range of miscibility of two metals A and B, the choice of the third element 
C should be such that its solubility parameter 5g should preferably be in 
the range of 5, ->d, or not greatly different from 5, or 8, and the electro- 
negativity X, should be as different as possible to X, and X zg Little 
data is available on the effect of a third element on the immiscibility range 
of two metals, but the following three examples from published data may 
be given as not inconsistent with the above hypothesis :—cadmium 
additions to bismuth—zine alloys, tin to lead-zine alloys and zirconium 
to uranium—thorium alloys, all of which give increased miscibility. The 
electronegativity values and solubility parameters at room temperature, 
given below for the three systems: 


(a) Cd Bi Zn 
) 45 48 58 

x 1-1 1:8 1-2 

| (b) Sn Pb Zn 
) 65 51 58 

x 1-7 1-5 1-2 

(c) Zr U Th 
r) 94 128 85 

3% 1:6 1:3 1-1 


§ 10. Discussion 


This theoretical analysis was made primarily to determine whether the 
simple Hildebrand immiscibility rule based on the concept of the solubility 
parameter, could be applied to metal systems, provided that due attention 
is paid to the decreased energy of the system arising from the difference in 
electronegativity of its components. Restrictions on the accuracy of 
the calculations were imposed by the lack of precise information on the 
properties which determine the solubility parameter and. its variation bie 
temperature, and by the difficulty of providing a precise figure for the 
electronegativity factor for a given system. Despite this, it has been 
shown that the idea of solubility parameter can be applied to liquid metal 
systems when combined with Pauling’s erect omega ayy concept based 
on the binding energy of two elements. Briefly, of 88 ons nie pee, 
which obey the simple Hildebrand rule, 72 agree with the modified hypo- 


thesis: of 206 immiscible systems chosen because they do not agree with 
\ 
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Hildebrand’s rule, 130 can be brought into line by allowing for electro- 
chemical differences. Several exceptions can be explained allowing for 
reasonable errors in the physical constants used in the calculations or by 
assuming that the electronegativities of certain elements depend on the 
metal with which they are alloyed. The dependence of the valency on 
the alloying element is not new to alloy theory but at the moment, it is 
difficult to see on what basis one chooses the electronegativity in a given 
case. 

A general picture of the liquid alloy state can be deduced from the 
above analysis. For non-polar liquids, the solubility parameter as 
defined by Hildebrand appears to be a satisfactory method of calculating 
the type of equilibrium between constituents. If two such liquids are 
mixed, then provided that their solubility parameters are not too different, 
the liquids will be mutually soluble in all proportions ; when the solubility 
parameters differ appreciably, a miscibility gap will be observed. Similar 
basic principles apply to liquid metal systems but the situation is com- 
plicated by an attraction between unlike atoms due to differences in their 
ionization potentials and electron affinities. Thus although the solubility 
parameters of two elements may be sufficiently different to give immisci- 
bility in the absence of an electronegativity factor, if they can form enough 
bonds of sufficient strength, miscibility will occur. The major problem 
associated with this hypothesis is to determine the maximum number of 
bonds which can form between two given elements. If the metals have 
the same size, they can be packed in the same way as the pure metals and 
the maximum number of bonds would appear to be six. As the size 
difference increases, simple packing becomes more difficult, the number of 
close atoms of each type decreases and the number of possible bonds will 
decrease to two when one element is about a third as large again as the 
other. With further increase in the size difference to about 60°%, the 
atoms can be packed into the interstices between the other type of atoms 
and the number of possible bonds is about six. 

Unfortunately, in the present state of knowledge, it is not possible to 
predict with absolute certainty whether two metals will be miscible or 
immiscible in the liquid state. For a given system, the possibility of 
miscibility depends on the estimated value for / and the accuracy to which 
this can be done. If the solubility parameters and their variation with 
temperature, and the electronegativity values were accurately known, it 
would be possible to get a reliable value for k. A knowledge of the 
bond formation between two elements with variation in their relative size 
proportions would enable one, not only to predict immiscibility, but to 
calculate the ranges of composition and temperature over which it occurs. 
It is considered, however, that the above type of analysis forms the best 
basis available on which to estimate the possibility of immiscibility. 
Further, it can be used as a guide to the best choice of element to add 
to a given system showing immiscibility to increase the range of 
miscibility. 
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The hypothesis as suggested above obviously requires further develop- 
ment. No distinction has been made between systems in which the 
elements are virtually completely immiscible and those giving immisci- 
bility over a narrow range only. An analysis of elements which give both 
limited liquid immiscibility and compounds would also be of interest. The 
value of the modified Hildebrand rule lies mainly in the success that can 
be obtained in predicting immiscibility in systems which have not been 
investigated experimentally. 
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CORRESPONDENCE 


On the Naor Relations between Madelung Constants for Cubic 
Ionic Lattices 


By F. G. Fumi and M. P. Tost 
Cattedra di Fisica Teorica dell’ Universita, Palermo, Italy 


[Received August 10, 1956] 


BerTavut (1954) has proved by the use of the Patterson function that 
linear relations exist between the Madelung constants for some cubic 
ionic lattices. He mentions that these relations were first found by 
Naor through an extension of the Frank method of lattice summation 
(Frank 1950), but Naor’s proof does not appear to have been published 
(Bertaut, private communication). Here we give two simple proofs of 
Naor’s relations. 
The Madelung constant «, of a cubic ionic lattice referred to the lattice 
parameter a is defined as 
—tha p3 €,O\r5) 2 1a) wig oe eee 


where ¢(r,) is the electrostatic lattice potential at the position of the pth 
ion of the unit cell}, «, is the charge of this ion in units of the smallest 
ionic charge in the lattice, and z is the number of molecules in the 
unit cell. 

To prove the Naor relations we write d(r,) as a sum of the electrostatic 
potentials in r, due to the simple cubic lattices which compose the crystal 
(Hund 1935): following Hund, we denote by ¢&(ayz) the electrostatic 
potential in the point xa, ya, za (a, lattice constant) due to a simple cubic 
lattice of positive unit charges, neutralized by a uniform distribution of 
negative charge. For the CsCl-type of lattice one has, with obvious 
notations, 

%a(C8Cl) = —Sal dose M*)—dose(X-)]= —adbosa(M*). . . (2) 
However, 


Souci M*)=(000) P95) ©. 2. 2. . (8) 
where (000) is the self-potential of the ‘ neutralized ’ simple cubic lattice 
of metal ions M+ and —:j(4$}4) is the potential due to the ‘neutralized ’ 


simple cubic lattice of halogen ions X-. So, using ne Hund identity 
(Birman 1955, Fumiand Tosi,to appear) 15(000)=3y( (400) + 375($40)+4(444) 


one obtains oa 4(CsCl) = 3a [eh (400) yt ab( £10)]. ‘ , c 4 " 


+ We call wnet cell the cell which contains an ion for each Bravais sub-lattice 
of the structure, namely the smallest lattice cell, 
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Similarly, if one uses the cubic symmetry of (xyz), the Hund identity 
for (000) and the identity (444) =44(444), one finds 


222 ), 
a&4(NaCl)=— 24 bxac(M*) —$nao(X~)]= —6arp($$0) (5) 
%,(zincblende)= — fal pzns(M+)—dgng(X~)]= —Balys(400) + 2 
0 q( Cal s)= — da[ 2p cqp,(M*+)— 2b oqp,(X-)]= — 6a 2yp( 300) + 34 
20.4(Cu20)= —3a[4$ou,0(M+)—4b.n,0(X--)] 
) 


a) ed 
ne Nie 
io 
— ~— 
= 

—~ nme 
~I [=>) 
~—— ae 


= —6a[3y(300)-++4y(340) +A] T. ©. > (8) 
These expressions for the «,’s lead at once to the Naor relations 
a 4(Cal,)=4«,(CsCl)+«,(NaCl) ek eee (oh) 
a4(Cal’,)= —a,(NaCl)+4«,(zincblende). . . . (10) 


From eqns. (4)—(8) it is also clear that no other independent, linear relation 
exists between the Madelung constants of the structures that we have 
considered. 

The Naor relations can also be derived very simply by considering the 
CaF’,-type of lattice with lattice constant a alternatively as a superposition 
of a CsCl lattice of lattice constant 4a and of a NaCl lattice with lattice 
constant @ or as a superposition of a zincblende lattice with double charges 
and of a NaCl lattice, both with lattice constant a. In the first case 
one has 


Poar,(M*)=2$06(M")+onao(Mt) 2. 2 6 6. . (I) 
boar, (X~)=2bo.(X~)+ byaciltet)=2¢os0(X7) - . . (12) 

from which eqn. (9) follows directly. Similarly in the second case 
doar, Mt) = 2b2ns(M*) + byaciltat) = 2b2ns(M*) - + (18) 
poar,(X—-)=2bzns(X~)+onanilMt) . . . - « - . (14) 


which give directly eqn. (10). 
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+ This equation is given explicitly in Hund’s paper, and lead Hund to 
discover the incorrectness of Emersleben’s value for a,(Cu,0). However, 
the incorrect value [a,,(Cu,0) =9:5044] is still quoted in the standard literature 
(e.g. M. Born and K. Huang, Dynamical Theory of Oa! Lattices (Oxford 
1954), p. 387) in place of the correct value a, (CuO) =10-2594. 
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Iv has usually been assumed that the isotopic spin impurities in highly 
excited states of a light nucleus will be much larger than those in the low 
lying states, due to the proximity of states of different isotopic spin 
(Gell-Mann and Telegdi 1953, Wilkinson 1953, MacDonald 1956). This 
is valid only if the matrix elements of the Coulomb interaction connecting 
these nearby states to the highly excited states are as large as are the 
matrix elements involved for the low lying states. 

It appears, however, that the large matrix elements of the Coulomb 
interaction are those coupling states belonging to different configurations. 
For an infinitely long range potential, all matrix elements connecting 
states of different isotopic spin would vanish, due to the orthogonality of 
the states. For the Coulomb interaction, the orthogonality shows up 
more for states of the same configuration, where there is cancellation 
amongst the different two-particle matrix elements, than for states 
belonging to different configurations, where the two-particle matrix 
elements are smaller but do not tend to cancel. As an example, we give 


Po Py ( ¥oH® ¥;) 
(Mev) 
¥(1s*1 p§[44]000) lstlp'1aaiee 2p)100) —0-48 


1s[3]420, ip*/44]000 2s)100) —0-36 


Y((1s41p7[43]441, 1d)001) s‘lp'[43]441, 1d)101) —0-02 


¥(( 

(| 

v(1 

((1s!1p7[421]441, 1d)101) 0-01 
Y((1s*1 p7[43]441, 28)101) 0-02 
¥((1s41p7[43]441, 2d)101) 0-11 
¥(((1s'lp*[42 010, 2p)441, 1d)101) =0-23 
(((1s41p*[42]012,,, 2p)441, 14)101) —0-35 
P(((I1st1p%[33|111, 2p)441, Id 


in the table some matrix elements connecting 7=0 and T=1 states of 
2C in Z-S coupling ; the triads of numbers are values of 7'SL, and the 
states) are appropriate to the ground state of 2C and to a 7’~0 state of 
about 20 Mev excitation. 

We have used 


He= — : rs) (Te) 
i<j 7] 
and radial wave functions for an oscillator potential well, with parameter 
a=2:2x10-8 em. The matrix elements connecting states differing in 
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the excitation of two particles are small compared with those connecting 
states differing in only one particle, as was pointed out by MacDonald 
(1956). 

On the basis of these matrix elements and those similarly calculated 
for other nuclei, it would seem that in general the main isotopic spin 
admixtures to a state are states of different configurations, with the amount 
of admixture independent of the excitation of the state. Thus isotopic 
spin may remain a useful quantum number up to quite high excitation 
energies. 

The experimental evidence on isotopic spin impurities is mainly in 
accordance with this view. To the extent that (d, d’) and (d, «) reactions 
proceed through a compound nucleus, the observed absence of transitions 
to final states with 7’ values different from that of the initial nucleus 
indicates absence of isotopic spin mixing in the states of the compound 
nucleus, which are at high excitation. Goward and Wilkins (1955) have 
interpreted their ’C(y,«) results in terms of very small isotopic spin 
impurities in the excited states of #2C. The impurities of various states 
of light nuclei as obtained from electric dipole transitions and competing 
reactions have been given by Wilkinson and collaborators (see Wilkinson 
1953, Wilkinson and Clegg 1956). Some apparently large impurities 
were found, the largest being about 20°, intensity of 7’=1 impurity in the 
mainly 7'=0 state at 6-88 Mev in B. It is also probable that there is a 
large impurity in the 12-69 Mev state of 44N (Barker and Mann 1957). 
Some large impurities are to be expected in a few cases when two com- 
bining states lie unusually close together, however, in view of what has 
been said above, the explanation of Wilkinson and Clegg (1956) for the 
large mixing in the !°B state is not tenable. 
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REVIEWS OF BOOKS 


Relativity: The Special Theory. By J. L. Syneu. (Amsterdam : | North- 
Holland Publishing Company, 1956.) [Pp. xiv+450.] Price 38 guilders. 


Tur keynote of this admirable treatise on the special theory of relativity is 
sounded in the preface where Professor Synge says “To understand a subject, 
one must tear it apart and reconstruct it in a form intellectually satisfying to 
oneself’. It is in this spirit of constructive criticism that the author gives a 
most complete account of the theory in which he discusses not only the familiar 
topics of the Lorentz transformation and its application to mechanics and 
electromagnetic theory but also investigates a number of interesting and 
important and puzzling questions which must also have come into the minds 
of students but which are not discussed in any of the standard works. Among 
these special topics which Professor Synge has investigated ab initio we may 
mention in particular the relativistic form of a discrete Newtonian system with 
action and interaction between its particles, the elastic collision of particles, 
the conservation of angular momentum, and models of uncharged ‘ particles ’ 
of finite energy constructed out of non-singular solutions of Maxwell’s equations 
in vacuo, together with model photons constructed according to the same plan. 
No difficulty is shirked in this work and detailed proofs are given of all the 
important theorems. In particular we may mention the particularly useful 
analysis of the Michelson—Morley experiment with special reference to the 
concept of the rigidity of the framework supporting the interferometer. The 
book is beautifully printed and will take its place as a standard treatise on the 
subject. Car 


Physical Mathematics. By CuEsTER H. Pacer. (Princeton, New Jersey: 
D. van Nostrand Co. Inc. 1955.) [Pp. x+329.] 


THERE is, among university students of physics and engineering, an increasing 
need for an understanding of advanced mathematical techniques. If they have 
been trained mainly in the experimental side of their subject, they often lack 
practice in mathematical discipline, and consequently find some of the standard 
mathematical text-books difficult to follow. Dr. Page’s book gives a very 
careful treatment of some selected topics, and could usefully be recommended 
to good graduate or advanced students who have such limited mathematical 
experience. The arguments are clearly presented, and illustrated with physical 
examples. 

Several chapters are devoted to the general theory of eigenfunctions and 
eigenvalues of differential equations, treated by variational methods and by 
separation of the variables, and illustrated by the properties of the Helmholtz 
equation. The detailed formulae are given for the transformations between 
the main coordinate systems in three dimensions, and the differential equations 
in the separated variables are classified in terms of their singularities. There 
is, however, very little mention of the special functions which satisfy these 
differential equations, and no mention of asymptotic approximations, the 
W.K.B. method, or the Airy Integral function. Green’s functions and integral 
equations are fully dealt with, and this is perhaps the most useful feature of 
the book. The last two-fifths of the book are devoted mainly to Fourier 
integrals and Laplace transforms. The illustrations are almost entirely from 
electrical circuit theory. The basic theory of the complex variable is treated 
only briefly in an appendix in one of the chapters. There is a full discussion 
of the relations between the resistive and reactive parts of an impedance or 
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admittance function, and between the amplitude and phase, known to com- 
munication engineers as Bode’s theorem. 

The book will be of most use in the fields of circuit theory and wave propaga- 
tion, but students of other branches of physics will find much to interest them, 
especially in the earlier chapters. KoGe8: 


Science and Information Theory. By L. Brittourn. (New York: Academic 
Press Inc. ; London : Academic Books Ltd., 1956.) [Pp. xii+320.] Price $6.80. 


In recent years Information Theory must have been a source of uneasiness for 
many physicists, as one of the subjects ‘‘ on which one hears so much and on 
which one ought to read up’’. Here at last is a book which one can whole- 
heartedly recommend to all physicists and scientists of an advanced level. 

_ The central theme of the book is the relation of the process of collecting 
information on physical systems to the 2d principle of thermodynamics. 

More than half of Brillouin’s book is devoted to the elucidation of the by no 
means quite simple connection between information and entropy, which he has 
widely generalized and brilliantly illustrated. It can now be said that the 
principle of indeterminacy and the ‘ negentropy ’ principle are the two funda- 
mental limitations of all physical measurements, but while the importance of 
the first is long recognized, the second has not so far received all the attention 
which it deserves. 1 regret only that (apart from a few brief remarks on p. 233) 
the author has not emphasized that the exorcism of the Maxwell Demon 
‘incidentally ’ proves also the inconsistency of classical physics with thermo- 
dynamics, that is to say that the 2d Printiple, like the 3d, derives its validity 
from quantum physics. This is a contention first put forward by Max Born, and 
not yet universally accepted. 

It is as yet hard to say what influence information theory will have on the 
future development of physics. Of all scientists the physicists are most reluctant 
to adopt statistical methods. The structural simplicity of their subject enables 
them to consider parameters one by one, and to reject conclusions obtained 
by elaborate statistical methods from complex material which does not show 
a clear result on inspection. But information theory, though it is in essence a 
statistical science, does not suggest giving up this noble prejudice. It suggests 
merely that in many cases a set of parameters can be more accurately 
obtained from complex experiments than from one-by-one tests. Moreover 
the theory gives upper bounds for the legitimate number of questions which 
a given experiment can be expected to answer. of 

Though the book will be most enjoyed by advanced physicists, it can be 
equally recommended to all scientists who had no previous contact with 
information theory, and to electrical engineers who had no intimate knowledge 
of statistical physics or with thermodynamics. The main results of com- 
munication theory are all derived from first principles, with just that level of 
mathematics which theoretical physicists have long accepted as their standard, 
that is to say with careful, even painstaking discussions of consequences, but 
without digging down into the foundations as deep as set theory. D.G 


Vacuum Deposition of Thin Films. By L. Houuanp. (London : Chapman & 
Hall Ltd., 1956.) [Pp. xx+-541.] Price 70s. 
Tur author is the Head of the Vacuum Coating Research Laboratory of 
Edwards High Vacuum Ltd., and his experience of the technology of his subject 
is reflected in the practical approach throughout the book and the range of 
topics chosen. A preliminary chapter on vacuum equipment is rigorously 
selective, and does not attempt to take the place of a text on vacuum practice. 
though some very useful advice is given. Chapters follow on the degassing of 
plastics in vacuo, various methods of cleaning surfaces, three on the evaporation 
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of various materials and one on the growth, structure and properties of vacuum - 
deposited films. Then the preparation of films for particular applications is 
described : Chapter 8 on metal films for electrical purposes such as metallized 
resistors (this section might have been expanded. with advantage), zinc coating 
of tissue for condensers and the deposition of electrodes on to crystal vibrators ; 
9 and 10 on optical interference films; 11 on evaporated aluminium films, 
especially for mirrors; 12 on evaporated coatings on plastics and lacquered 
surfaces: 13 on shadow-casting and the production of surface replicas for 
electron microscopes; 14 on cathodic sputtering ; 15 and 16 on oxide films, 
including electrically-conducting oxides. There is a valuable appendix 
collecting together much information on deposition characteristics and refrac- 
tive indices of inorganic compounds. 

Throughout the book the treatment is selective and critical in the light of the 
author's own experience. It is up-to-date and well documented. The list 
of 557 references at the end is numbered in the order of citation in the text, but 
the author index refers to this list and enables particular papers to be traced. 
A short bibliography of relevant books is included in the author's preface. 
Where needed the fundamental physics of the processes involved is discussed, 
usually with a sound grasp of the main principles. It would have been outside 
the scope of the book to have delved deeply into theory, but it would have 
added considerably to the value of some sections to have been more quantitative 
and precise, for example when discussing adsorption, re-evaporation and life- 
times of surface atoms. As Professor Tolansky rightly remarks in his appre- 
ciative foreword, the subject of vacuum deposition is still as much an art as a 
science and anyone taking up experimental work in it needs an experienced 
guide. Mr. Holland through this book supplies that need. Other workers in 
this and neighbouring fields will be able to add information from their own 
experience, e.g. the effects of chlorine liberated from some glasses when baked, 
but they will be grateful to the author for providing so much practical advice, 
so many data and such a good selection of references to other work. F. A. V. 


The Atomic Nucleus. By R. D. Evans. (London: McGraw-Hill.) [Pp. xv+ 
972.] Price 109s. 


Turs massive book is based on the graduate course in nuclear physics which the 
author has given at the Massachusetts Institute of Technology for the past 
twenty years; there can be indeed little in it that has survived those two 


which is that the book consciously or unconsciously tries to be a text-book 
containing only what is unquestionably right—while pretending to give an account 
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of present-day nuclear physics. The two are incompatible. As a result 
author usually stops short of any novel, MeOOReOdaLA or controversial eee 
and the reader gets little feeling for the way in which nuclear physics is in fact 
developing today. Despite these remarks this is by far the best book on experi- 
mental nuclear physics presently available. 

The references are indexed into the main text—a most useful feature—and 
the index itself is very good (indeed on the Lamb shift it is more informative 
than the text). For a first edition the book is remarkably free from misprints 
and trivial errors and it is a monument to the care of the author. The large 
number of problems set in the text is a useful feature. 

We finally note that a book which figures in a series whose title contains the 
word “ International’? should not date scientific events from local history. 
Mendeleev’s establishment of the periodic table is dated by Dr. Evans as 
“ About the time of the American Civil War”. Enquiries by the present 
reviewer among a statistically significant number of his colleagues have enabled 
him to fix this event in the year 1842+30. However lest the reviewer appear 
ungrateful he must record his thanks for the very precise information gleaned 
from p. 755 that “on June 5, 1951, Dom DiMaggio had a batting average 
of 0-359, had been at bat 189 times in 44 games and had hit safely at least once 
in each of his last 25 consecutive games’. Whatever that means. D. H.W. 


Recent Advances in Science. By M. H. SHamos and G. M. Murpuy. (New 
York: Interscience Publishers, Inc.; London: Interscience Publishers 
Ltd.) (Published jointly with the New York University Press.) [Pp. xi+ 
384.] Price $7.50. 


THis is an interesting volume containing articles on recent advances in physics 
and applied mathematics written at a level which is meant to be acceptable to 
all scientists and engineers in whatever field they work. Most of the authors 
are extremely eminent and leaders of their subjects ; for instance Bethe writes 
on nuclear structure, Morse on operations research and Shockley on transistors. 
On the whole the authors seem to the present reviewer to have been successful 
in writing at the right level though of course there are the usual anomalies, 
such as a discussion in one article of what is meant by a wave and the use of 
Maxwell’s equations in another. None the less the book can be recommended 
to those for whom it is meant, namely anyone with a general knowledge of 
science who wishes to know what is happening in certain branches of pure and 


applied physics. N. F. M. 


Introduction to Solid State Physics. Second Edition. ByC.Krrren. (London: 
Chapman & Hall.) [Pp. xvii+617.] Price 96s. 

A sEcoND edition of Kittel’s well known “ Introduction to Solid State Physics ” 

is welcome. It is about two hundred pages longer than the first edition and 

contains a fuller account of the theory of energy bands and new material on 

alloys, semiconductors and imperfections in solids. NL. 


Photoconductivity Conference held at Atlantic City, November 4-6, 1954. By 
R. G. Breckenripcr, B. R. Russern and E. E. Hann. (New York : 
John Wiley & Sons Inc. ; London: Chapman & Hall Ltd.) [Pp. xiii+653.] 
Price 108s. 

Tus volume contains the papers presented at a conference in Atlantic City 

in November 1954 sponsored by the University of Pennsylvania, the Radio 

Corporation of America and the U.S. Office of Naval Research. It contains 

papers on a very wide range of subjects, theoretical and applied, and will be of 

considerable value to workers in this field. N.F.M. 
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CORRIGENDUM 


The Effect of Free Electrons on Lattice Conduction 


By J. M. Ziman 
Cavendish Laboratory, Cambridge 


In a recent paper (Ziman 1956), the mean free path of phonons scattered 
by electrons was calculated, but the phonon creation processes, which 
tend to balance out the annihilation processes to some extent, were 
neglected. Their effect is to increase the relaxation time by an amount 
which, as may be shown by standard methods (Wilson 1953), corresponds 
to the deletion of the factor (1—e-*)~! in eqns. (9) and (19) of my paper. 
The curves for 4 have been recalculated, but turn out to be very similar 
to those drawn in fig. 1. The low temperature asymptote is almost 
unchanged, having the form 


(a) ToT 3) as T 0, d= 02088 tas are od 


The minimum values of ¢ are a few per cent larger, but come at the same 
values of T/T’. Only at high temperatures is there a significant difference, 
with 


= ow A/a [UN 32 (oP T's 
(c) as T'+ o, oS ay (7) (7-)exn (7) 


( 475 1 Ty 3/2 4 T » 3/27" 5 ce 
a2 eS ae See RTT at | By Bards) 
aia ade a(F) +=, (7) ‘a HAE! 


The mean free path thus rises as 7°’? which is much steeper than any 
observed effect. 


I am indebted to Dr. R. Stratton for pointing out this error. 
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